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ABSTRACT 

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space 
telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Markarian 266. The SW 
(Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and SO/a (pec), 
respectively. Both companions are more luminous than L* galaxies and they are inferred to each 
contain a w 2.5 x 10 8 Mq black hole. Although the nuclei have an observed hard X-ray flux ratio of 
,fx(NE)/fx(SW) = 6.4, Mrk 266 SW is likely the primary source of a bright Fe Ka line detected 
from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. 
Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous 
H 2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding 
the NE nucleus. A soft X-ray emission feature modeled as shock-heated plasma with T <~ 10 7 K 
is co-spatial with radio continuum emission between the galaxies. Mid-infrared diagnostics provide 
mixed results, but overall suggest a composite system with roughly equal contributions of AGN and 
starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been 
detected, with most having estimated ages less than 50 Myr. Detection of 24 /im emission aligned with 
soft X-rays, radio continuum and ionized gas emission extending ~34" (20 kpc) north of the galaxies 
is interpreted as ~ 2 x 10 7 M Q of dust entrained in an outflowing superwind. At optical wavelengths 
this Northern Loop region is resolved into a fragmented morphology indicative of Rayleigh- Taylor 
instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" 
phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, 
and rapid gas consumption in luminous dual AGNs with kpc-scale separations early in the merger 
process may explain the paucity of detected binary QSOs (with pc-scale orbital separations) in 
spectroscopic surveys. An evolutionary sequence is proposed representing a progression from dual to 
binary AGNs, accompanied by an increase in observed L x /Lj r ratios by over four orders of magnitude. 

Subject headings: galaxies: interactions — galaxies: nuclei — galaxies: active — galaxies: Seyfert — 
galaxies: starburst — galaxies: star clusters: general 



1 Infrared Processing & Analysis Center, MS 100-22, 
California Institute of Technology, Pasadena, CA 91125; 
mazz , b chan@ ipac . caltcch.edu 

2 ICREA and Institut del Ciencies del Cosmos (ICC), Uni- 
versitat de Barcelona (IEEC-UB), Marti i Franques 1, 08028 
Barcelona, Spain; kazushi.iwasawa@icc.ub.edu 

3 Department of Physics and Astronomy, State University 
of New York at Stony Brook, Stony Brook, NY 11794-3800; 
vavilkin@grad.physics.sunysb.edu 

4 Spitzer Science Center, MS 314-6, Califor- 
nia Institute of Technology, Pasadena, CA 91125; 
lee jhhowclljason, inami, haan@ipac.caltech.edu 

5 Department of Astronomy, University of Virginia, Char- 
lottesville, VA 22904-4325; aevans, gcp8y@virginia.edu 

6 National Radio Astronomy Observatory, 520 Edgemont 
Road, Charlottesville, VA 22903-2475; dkim,dfrayer@nrao.edu 

7 University of Oregon, Physics Department, Eugene OR, 
97402; nuts@bigmoo.uoregon.edu 

8 Department of Astronomy, Cornell University, Ithaca, NY 
14853; spoon@isc.astro.cornell.edu 

9 CSIRO Astronomy & Space Science, ATNF, Marsfield NSW 
2122, Australia; scbastian.haan@csiro.au 

10 NASA Herschel Science Center, MS 100-22, Cal- 
ifornia Institute of Technology, Pasadena, CA 91125; 



lord,apple,cxu@ipac.caltech.edu 

11 Astronomy Department, California Institute of 
Technology, MS 249-17, Pasadena, CA 91125; ja- 
son.mashall@caltech.edu, ap@astro.caltech.edu 

1 Subaru Telescope, National Astronomical Observatory of 
Japan, Hilo, HI 96720; cmishida@mac.com 

3 Institute for Astronomy, University of Hawaii, 2680 Wood- 
lawn Drive, Honolulu, HI 96822; sanders, vivian@ifa.hawaii.edu 

14 Byurakan Astrophysical Observatory, Academician of 
National Academy of Sciences of Armenia, Byurakan, 378433 
Aragatsodn Province, Republic of Armenia; khache@bao.sci. am 

15 The Observatories, Carnegie Institution of Wash- 
ington, 813 Santa Barbara Street, Pasadena, CA 91101; 
barry@obs.carncgiescience.edu 

16 Max-Planck-Institut fr extratcrrestrischc Physik, Postfach 
1312, 85741 Garching, Germany; sturm@mpc.mpg.de 

17 Department of Astronomy, University of Maryland, College 
Park, MD 20742; veilleux@astro.umd.edu 



2 



Mazzarella et al. 



1. INTRODUCTION 

1.1. LIRGs, ULIRGs and GOALS 

Luminous Infrared Galaxies (LIRGs; 10 11 L < Li r < 
10 12 Lq) and Ultraluminous Infrared Galaxies (ULIRGs; 
Lir > 10 12 £©) are intriguing objects with widespread 
implications for galaxy evolution. They contain the high- 
est known rates of star formation, they exhibit a high 
frequency of active galactic nuclei (AGNs) and large- 
scale outflows (superwinds), and mounting evidence in- 
dicates that the majority of local (U)LIRGs represent 
an evolutionary process involving the transformation of 
major mergers between dusty, gas-rich disk galaxies into 
elliptical galaxies hosting cla ssical UV-excess QSOs or 



powerful radio galaxies (e.g., Sanders & Mirabel 1996 
Veilleux|[2006l ). At redshifts of z ~ 1, LIRGs have a 
higher space density than ULIRGs and d ominate the to 



tal st ar- format ion density at that epoch ( |Le Floc'h et al. 
20051). At z - 2 the contributions of LIRGs and ULIRGs 



to the total IR lum inosity density are about equal (e.g., 
Caputi et al. 2007 1 . Since LIRGs and ULIRGs were much 



more common in the early universe, these populations are 
fundamental in understanding both star formation and 
galaxy evolution. 

The Great Observatories All-Sky LIRG Survey 
(GOALS Y°\ is utilizing imaging and spectroscopic data 
from NASA's Spitzer, Hubble, Chandra and GALEX 
space observatories, in combination with ground-based 
observations, in a comprehensive study of more than 200 
of the most luminous in frared galaxies in the local uni- 
verse (Armus et al. 2009). The sample consists of 181 
LIRGs and 23 ULlRGs that form a statistically complete 
subset of the flux- limited IRAS Revised Bright Galaxy 
Sample (RBGS), which itself is comprised of 629 extra- 
galactic objects with 60 /im flu x densities above 5.2 4 Jy 
and Galactic latitudes |6| > 5° (Sa nders et al.|[2003 ). 



1.2. Dual AGNs in Major Mergers 

A key scientific driver for GOALS is the investigation 
of the relative time-scales and energetics of active star 
formation and AGN phenomena during different phases 
of the merger sequence. It is now widely accepted that 



(SMBH) in their centers (e.g., 
with a SMBH mass proportiona 
bulge (M S MB H /M bulge ps 0.2%; 



Ferrarese & Ford 


2005) 


to the mass of the stella 


Marconi k. Hunt 


2003 



ing of paired SMBHs during major mergers. While the 
nuclei are still separated by kpc scales, how common is it 
for both SMBHs to have accretion rates high enough to 
produce energetically significant dual AGNs, as opposed 
to one or both nuclei being powered predominantly by 
a burst of star formation? How do AGN characteristics 
depend on properties of the host galaxies and dynamics 
of the encounter? Are the fuel supplies and accretion 
rates sufficient to sustain two luminous AGNs well into a 
true binary stage (e.g., binary QSOs), when the SMBHs 
are closely bound (r < 100 pc) in Keplerian orbits inside 
a dynamically relaxe d (elliptical) merger remnant (e.g., 



Colpi & Dotti 2009)? Is there something special about 
the physical conditions in systems that host dual AGNs, 
or is their lack of detection in some (U)LIRGs merely 

18 http : // goals . ipac .caltech.edu / 



a matter of observing them during the right stage prior 
to merging (i.e., an evolutionary timing coincidence), or 
accounting for differences in nuclear dust obscuration? 
What is the relative importance of AGNs and star for- 
mation in the energetics of the "supcrwind" phenomena 
that appe ars to be ubiqu itous in (U)LIRGs (e.g., see the 



review by Veilleux 



2006)? 



Although (U)LlRGs are predominantly systems in- 
volving major mergers, extensive observations utilizing 
optical and infrared spectroscopy have so far turned 
up surprisingly few systems containing confirmed dual 
AGNs. This has a number of possible explanations. 
First, (U) LIRGs contain large quantities of centrally con- 
centrated dust that can effectively obscure one or both 
AGNs at optical , near-infrared, and even mid-infrare d 
wavelengths (e.g. Laurent et al.|2000 Spoon et al.| 2004). 
Therefore, circumnuclear star formation can dominate 
the observed spectra while in reality a powerful AGN 
may be buri ed inside. Second, about 30% of observed 
LIRGsJe^ JVeilleux et al.|1995[ ), and -40% of ULIRGs 
(e.g., |Veilleuxet al.|1999| ) are classified as LINERs based 
on optical spectroscopy. It has proven difficult to dis- 
tinguish between low-luminosity AGNs (photoionization 
from radiation due to accretion onto a SMBH), photoion- 
ization from very hot stars, and shock heating (induced 
by SNe and starburst-driven super winds) as the primary 



source of ionization in LINERs (e.g., Veilleux et al.|1999 ) 
Recent results indicate that most nearby LllNERs are 
dominated by photo-ionization rather than shock heat- 
ing, and that they are an important class of AGNs dis- 
tinguished pri marily by a lower accretion rate than in 
Seyfert nuclei ( [Kewley ■ et al |2006[|H~o|2008| ). This implies 
that the frequency of dual AGIN s may be much higher 
than inferred to date. Only recently, with the capabil- 
ity of the Chandra X-ray Observatory to penetrate their 
extensive dust cocoons with high-resolution hard X-ray 
observations, have suspected dual AGNs been confirmed 
in a small number of (U)LIRGs. Perhaps the best known 
example is NGC 6240, in which hard X-rays and strong 
Fe Ka emission lines indicate the presence of two AGNs 
with a p rojected separation of l'/4 (1 kpc) ( |Komossa et| 
al.|2003j >. A second case is the ULIRG Mrk 463, a system 
first pointed out as an active double-nuc leus galaxy by 



Petrosian, Saakian & Khachikian (19781. The eastern 
component undoubtedly hosts a dust-enshrouded Type 
1 AGN that is powering apparently young radio jets 



(Mazzarella et al. 
tral classifications 



1991a) 



but conflicting optical spec 
rom various studies left the nature of 
the western nucleus in doubt (LINER, Seyfert 2, or star- 
burst/ AGN composite). Dual AGNs in Mrk 463 have 
been confirmed recently via detection of two c ompact 
luminous so urces of hard X-rays using Chandra ( jBianchi 
et al.||2008|). To date, very few systems in the GOAL 



sample have been confirmed to contain dual AGNs via 
X-ray observations. The first system is NGC 6 240, and 



th e seco nd is Arp 299 (NGC 3690 + IC 694) flBallo et 
|al.||2004[ ) . This article presents an in-depth investigation 
of a third system in this rare class of objects that has 
considerable importance for understanding the evolution 
of galaxy mergers and their remnants. 

1.3. Complex Phenomena in Mrk 266 

One of the most remarkable LIRGs in the local Uni- 
verse is Mrk 266 (NGC 5256). The system was first 
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called out as an extraordinary object in the Markarian 
Survey (First Byurakan Sky Survey) of ultraviolet-excess 
galaxies due to the presence of t wo nuclei within what 



appeared to be a singl e galaxy (Petrosian, Saakian & 



|Khachikian||1978 



vestigation at op 



1979). In a detailed spectroscopic m- 



lcal wavelengths, |Petrosian, Saakian fe| 



Khachikian ( 1980 ) showed that the nuclei have a radial 
velocity difference of 280 km s _1 , and based on the as- 
sumption that the nuclei revolve about a common cen- 
ter of gravity, masses of 7 x 10 9 and 10 x 1O 9 M were 
estimated for the SW and NE nuclei, respectively. As 
for all objects in the Markarian Survey, Mrk 266 stood 
out because its UV-to-optical flux ratio is higher than 
in normal galaxies. However, Mrk 266 is also a LIRG, 
with an infrared luminosity of L ir = 3.6 x 10 11 Lq (8 
- 1000 //m, a s estimated from f lux densities in all four 
IRAS bands: Armus et al.| [2009). As shown in the cur- 
rent study (j]3.5p, the global spectral energy distribution 



indicates the bulk of its energy is emitted in the mid- 
to far-infrared; much more energy is being emitted at 
the peak of the SED than is escaping in UV radiation 
(i/L v (70/im)/i/L,,(0.2/im) « 20). 
Mrk 266 contains a filament ary nebula of ionize d hy- 



drogen ~30 kpc in diameter (Armus et al. |1990[ ) sur- 
rounding the two nuclei. A luminous X-ray nebula ^100 
kpc in extent detected by ROSAT (also surrounding the 
two nuclei) and complex kinematics derived from opti- 
cal spectroscopy revealed one of the most spectacular 
examples of an o utflowing, starburst-driven superwind 
( Wang et al.|1997 ). In addition to radio continuum emis- 
sion from the two nuclei, a bright radio source located 
between the nuclei was detected and interpreted as en- 
hanced synchrotron emission induced by extensive shocks 



at the interface of the mergin g galaxies ([M azzarcll a et al. 

Imaging in Ha + [Nil] jMazzarella et al.|1988| |Ar-, 
mus et al.|l990j, [O III] A5007 ([Hutchmgget al.|1955f and 



soft X-rays ( |Kollatschny fe Kowatsch|1998[ ) also revealed 
a bright loop of emission extending ~24"(17 kpc) to the 
north-east and con nected to the SW (Seyfert 2) nucleus 
(Wang et al. 1997). Recent Chandra observations have 



resolved X-ray emission from both nuclei, and detected 
diffuse emission asso ciated with the northern f eature and 
between the nuclei (Brassington et al. 20071. The na- 



ture of the northern emission region (hereafter called the 
Northern Loop) has remained controversial, wit h sugges- 
tions inc luding a component of the superwind (|Wang et| 
al.|1997| ), a "jet" from an AGN (|Kollatschiry fc Rowatsch| 
1998J7and a tidal tail (IBrassington et al.||2007p. 



Mrk 266 is worthy of further detailed study because 
it manifests complex processes that are rarely detected 
simultaneously, presumably due to their relatively short 
time scales: vigorous star formation, a candidate dual 
AGN (depending on the ionizing source of the LINER) , 
a powerful large-scale superwind, and substantial radio 
continuum and X-ray emission between the nuclei. The 
relative proximity of the system permits close-up inves- 
tigation of an important stage in the evolution of major 
mergers that may involve the transformation of LIRGs 
into ULIRGs, AGN feedback with the ISM, and impli- 
cations for fueling (or quenching) of accreting SMBHs 
which galaxy evolution models predict to have occurred 
in large numbers during earlier epochs. 
In this study, new observations and re-processed 



archival data from the Spitzer, Hubble, Chandra, XMM- 
Newton and GALEX space telescopes are analyzed in 
combination with multi-wavelength ground-based data. 
Imaging, photometric, and spectroscopic data are pre- 
sented in f|2] In glthc new data are interpreted to obtain 
new i nsigh ts into the nature of the colliding/mergin g sys - 
tem (j |3.1[ ), the nuclei and circumnuclear regions (S3.2), 
the reg ion b etween the nuclei ({ 3.3 ), the extensive s uper - 
wind (j |3.4[ ), the SEDs of the major components (?3.5), 
newly detected star clust ers ( j3.6), and the molecular 
gas properties Q3.7). In { 3.81 Mrk 266 is examined in 



comparison with other interacting systems with strong 
radiation sources between the galaxies, and a sequence 
is proposed that may represent an evolutionary progres- 
sion from dual AGNs (with kpc-scale separations) to bi- 
nary AGNs (with parsec-scale orbits). Due to the large 
number of new measurements spanning many spectral re- 
gions, the basic observations and data reductions are de- 
scribed in fj2j but most of the corresponding figures and 
tables are displayed alongside their interpretation and 
analyses in <|3j Conclusions are summarized in ^J4j A sys- 
temic heliocentric recession ve locity of 8353 ± 13 km s _1 
(Ide Vaucouleurs et all 119911) is adopted, corrected to 



8825 ± 2 2 km s" 1 via the flow model of |Mould et aL. 
( 2000a|b ) that accounts for three major attractors in 



the local Universe. We adapt cosmological parameters 



H n = 70 km s~ + Mpc 



fl M = 0.28, and ft v = 0.72 
( jHinshaw et al.| 20091). The corresponding luminosity 
distance to Mrk 266 is 129 Mpc (distance modulus 35.55 
mag), and the spatial scale is 0.59 kpe/arcseep 9 ] 

2. OBSERVATIONS 
2.1. Optical and Near-Infrared Imaging with HST 

Mrk 266 was observed 2005 November 17 with the Ad- 
vanced Camera for Surveys (ACS) Wide Field Channel 
(WFC) as part of a HST Cycle 14 GOALS program to 
investigate 88 LIRGs with L ir > 10 114 L o (PID 10592, PI 
A. Evans). The system was imaged with an integration 
time of 22.5 minutes in the F435W filter (B band, 0.435 
/im, image #J9CV48010) and 13 minutes in the F814W 
filter (I band, 0.814 fim, image #J9CV48020). The WFC 
has a field of view of 202" x 202" with pixels subtending 
0705. The images were processed through the standard 
STScI pipeline to remove instrumental signatures, sub- 
tract a dark frame, apply flat fielding, remove geometric 
distortion, and apply flux calibration. Residual cosmic 
ray removal, background subtraction, and accur ate as- 
trometric calibra tion using sources from 2MASS ( |Skrut-| 
skie et al. 2006) were also performed using procedures 
detailed by Evans et al. (2012, in preparation). Figure [I] 
presents a composite color image constructed using the 
B and I band data, and labels identify major features 
that are analyzed throughout this article. 

Near-infrared observations at 1.6 /im (F160W, H band) 
with the NIC2 aperture of NICMOS were acquired on 
1997 September 13 in Cycle 7 (PID 7328, PI M. Malkan). 
The image (#N44B40010) was retrieved from the HST 
archive, corrected for bad pixels, cosmic rays were re- 
moved, and astrometric calibration was applied by align- 
ing the galaxies with the ACS images. The central re- 
gions of the SW and NE galaxies in these three HST 



19 Provided by NED at http://ned.ipac.caltech.edu/ 



4 



Mazzarella et al. 
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(LINER) ' 



Filaments 
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Northern Loop 
(Superwind) 



Star Clusters & 
Associations 




Mrk 266 SW 
(Seyfert 2) 




Background 
galaxy pair ? 



5.0" - 3.1 kpt 



Figure 1. HST ACS I band and B band images combined to form a color composite image of the Mrk 266 system and its tidal features. 
North is up and east is left. The image field of view is lf4 x 1.'6, with a scale bar in the lower left corner. The color table has been chosen 
to maximize contrast for faint extended structures. The limiting surface brightness is 25.0 and 24.5 mag arcsec - 2 in the B and I band, 
respectively. Labels identify major structures that are studied throughout this article. 
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the data r eduction procedures can be found in |Iwasawa et| 



J 



I. j^' 



Z0" (12 kpc) 



Figure 2. The result of stacking B, V, and I band images ob- 
tained with the University of Hawaii 2.2 m telescope (effective see- 
ing ~T.'5 FWHM). The field of view is 150 kpc x 150 kpc (4f2 
x 4'2) , and a scale bar is provided in the lower left. Faint, low 
surface-brightness emission spanning 2.'9 (~ 103 kpc) is detected. 

band are displayed in context with their analysis in j]3.1| 




2.2. Deep B+V+I Band Ground-Based Imaging 

To complement the HST high-resolution imagery and 
reveal the faintest emission possible, Figure [2] shows the 
result of spatial registration and stacking of B, V, and 
I band images obtained on 1998 December 22 and 2000 
May 04 with the University of Hawaii 2.2 m telescope 
(effective seeing ~l'/5 FWHM). Det ails of t he da ta ac- 
quisition and reduction are given by Ishida (2004). The 
data are presented using Gaussian equalization to con- 
trast the low and high surface-brightness features in the 
same view. Highly asymmetric, low surface-brightness 
emission can be seen spanning «103 kpc (2(9) from the 
SE to NW extremities, which is approximately twice the 
extent of the faintest emission visible in the HST ACS 
data, SDSS, and previously published imagery. 

2.3. X-Ray Observations with Chandra and 
XMM-Newton 

2.3.1. Chandra Observations 

Chandra observations were obtained with the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S) on 2001 
November 02 (PI A. Read) with a total integration time 
of 19.95 ks. The spatial resolution is ~ 1" FWHM. These 
observ ations were first presented by |Bras sington et al. 
(2007). In order to make refined measurements, high- 
light additional features, and analyze the X-ray prop- 
erties alongside observations at other wavelen gths , the 
data were taken from the Chandra Data Archival (OB- 
SID #2044) and re-processed using standard procedures 
in CIAO 4.0. The data were corrected for detector re- 
sponse and converted to flux density (f„) units. Details of 



al.| ( 2011[ ). Figure |3] shows the resulting full band (0.4 - 7 
KeV) image, the hard band (2-7 keV) image, a smoothed 
version of the full b and data, and the 20 cm radio con- 



tinuum image from Mazzarella et al. fi988). Figure 



illustrates apertures with parameters listed in Table T 
that were used to define regions of interest, and Figure 
shows the X-ray spectra constructed in these regions. 



(a) 0.4 - 7 keV 

■ ■■■■■ ■ ■. ■ ■ v _5k . . 


(b) 2 - 7 keV 




■ 






(c) 0.4 - 7 keV 

V 


(d) 20 cm 

• 


10" 





Figure 3. The X-ray and radio continuum emission in Mrk 266: 
(a) Chandra 0.4 - 7 keV (full band) X-ray image; (b) 2 - 7 keV 
(hard band) X-ray image; (c) 0.4 - 7 keV X-ray data (same as 
panel a) smoothed with a 2-pixel G aussian kernel; and (d) 20 cm 
radio continuum image (VLA) from Mazzarella et al. (1988). The 
scale bar indicates 10". 



2.3.2. XMM-Newton Observations 

X-ray observations with XMM-Newton, reaching sig- 
nificantly higher energies (10 keV) than Chandra (7 keV), 
were acquired with the European Imaging Photon Cam- 
era (EPIC) on 2002 May 15 (PI A. Read). The data 
obtained with the EPIC pn-CCD camera were utilized, 
providing an angular resolution of ^6" FWHM and a 
spectral resolution of (E/SE) ~ 40. The observation 
was performed with the Thin Filter in the extended full- 
window imaging mode. The data set (#0055990501) was 
taken from the XMM-Newton Science Archivep"| and pro- 
cessed with standard procedures in the software packages 
SAS 8.1 and HEASoft 6.6 utilizing the latest calibration 
files (as of 2009 January). Time intervals with high parti- 
cle background levels were discarded, leaving useful data 
with an effective exposure time of 13.2 ks. The count 
rate in the 0.4-10 keV band is 0.19 count s _1 . The X-ray 
emission is only marginally resolved by XMM, therefore 
the spectrum for the total system emission within a large 
aperture of 45" radius is produced here. Figure [6] shows 
the resulting XMM spectrum of the integrated emission 
from the system. Key features include Si XIII Ka (1.853 
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Table 2 

X-Ray Measurements 



Figure 4. Elliptical apertures used to define regions of inter- 
est for extraction of X-ray spectra from the Chandra data: Mrk 
266 NE ("NE"), Mrk 266 SW ("SW"), the central region between 
the nuclei ("Cen"), the Northern Loop ("Loop"), and surrounding 
diffuse emission ("Dif"). The two ellipses marked with red lines 
indicate regions subtracted from the ellipse labeled "Dif" to esti- 
mate the diffuse component. Parameters for the regions are listed 
in Table [I] A circular aperture of radius 33" used to measure the 
integratea system flux is not illustrated. 



Table 1 

Elliptical Apertures for X-ray Measurements 



Region 


RA, DEC (J2000) 


a" b/a 


PA° 


(1) 


(2) 


(3) (4) 


(5) 



SW 


(13 


38 


17 


31, +48 


16 


31 


9) 


1 


7 


1 


00 




NE 


(13 


38 


17 


82, +48 


16 


41 


2) 


1 


7 


1 


00 




Center 


(13 


38 


17 


66, +48 


16 


35 


8) 


6 


9 





47 


155 


Loop A 


(13 


38 


17 


15, +48 


16 


57 


2) 


4 


5 





32 


10 


Loop B 


(13 


38 


17 


83, +48 


17 


01 


4) 


4 


7 





68 


75 


Diffuse a 


(13 


38 


18 


71, +48 


16 


34 


2) 


22 


5 





84 


80 


Diffuse -1 


-(13 


38 


17 


49 , +48 


16 


36 


3) 


4 


6 





61 


25 


Diffuse -2 


-(13 


38 


18 


07, +48 


16 


33 


6) 


2 


5 


1 


00 




Total 


(13 


38 


17 


66, +48 


16 


35 


8) 


33 





1 


00 





a The diffuse emission is approximated by the flux within this large 
ellipse after subtraction of the flux in the two small ellipses in the 
following two rows, which arc proceeded by a minus sign in this table 
and marked with red lines in Figure |4| 

keV) emission, as also observed in the Chandra spectra 
of the two nuclei (Fig. [5| , and strong Fe Ka line emission 
at 6.4 keV. 

2.3.3. X-ray Measurements 



Spectral measurements and results from continuum 
model fitting of the Chandra and XMM-Newton obser- 
vations are presented in Table [2] The sum of the emis- 
sion in the elliptical regions misses some low surface- 
brightness emission that is recovered in measurements 
of the total system in a circular aperture with a radius 
of 33". The total system fluxes measured by Chandra 



Region 


F(Fe Ka) 


kT (soft) 


Fx (soft) 


F x (hard) 


(1) 


(2) 


(3) 


(4) 


(5) 


Mrk 266 SW 


1 2+ - 7 
l - Z -0.5 


r, 77 +0.07 
u -' '-0.08 


2.2 


2.2 


Mrk 266 NE 


< 1.6 


n 72+ ' 11 
u -'^-o.io 


3.0 


11 


Between Nuclei 




u - DD -0.07 


2.6 


1.0 


Northern Loop 




n oer + 0.05 
u - 5o -0.05 


1.4 




Diffuse Region 






4.1 




Total {Chandra) 


1 =+0.8 




18 


20 


Total (XMM) 




18 


23 



Note. — Measurements of the Chandra spectra in regions of 
interest; they are not corrected for absorption. Column (1): Region 
of interest, as illustrated in Fig. H] Column (2): Flux of the Fe 

Ka line in UnitS Of 10 _ 6 ph CHI™ = _ 1 Tll » r«t.framo onor-mr 



The rest-frame energy 



measured from the XMM spectrum for this feature is 6.50^g'Jg 
keV. The value for the NE nucleus is a 2<j upper limit. Column 
(3): Temperature (kT in units of keV) of the soft (0.5-2 keV) X-ray 
emission with la uncertainties. Column (4): Flux in the soft (0.5-2 
keV) X-ray band in units of 10 — 14 erg cm -2 s — x ; the estimated 
relative uncertainty is 20%. Column (5): Flux in the hard (2-7 
keV) X-ray band in units of 10 — 14 erg cm -2 s" 1 ; the estimated 
relative uncertainty is 30%. No flux is detected above 2 keV in 
the Northern Loop or in the Diffuse Region. Mrk 266 SW shows 
some emission above 2 keV, but most of the detected flux appears 
to originate in the Fe Ka line. The NE and SW galaxies also both 
show strong Si XIII emission at 1.8-1.9 keV. 

are, within the uncertainties (i.e., 20% in F x (soft) and 
30% in F x (hard)), in agreement with the XMM total flux 
measurements listed in Table [2] There may be some con- 
tribution of photoionized gas (by the AGNs) in the soft 
X-ray emission, but the S/N is insufficient to fit this com- 
ponent; therefore, the temperatures given in Tableware 
obtained from a model fit consisting of thermal emission 
due to star formation. Mrk 266 NE has an AGN compo- 
nent, which is modeled by an absorbed power-law with 
photon index 1.8 (fixed) and N H = 8.2±^| x 10 22 cm" 2 ; 
correcting for this absorption suggests an intrinsic hard- 
band flux of 2.7 x 10 -13 erg cm -2 s^ 1 . There is a tail 
of hard X-ray emission in the source between the nuclei, 
but its spectral shape is not well constrained; the 2-7 keV 
flux given here is mostly due to this hard tail. 

2.4. Ultraviolet Imaging with GALEX and 
XMM-Newton OM 

Ultraviolet images of Mrk 266 were obtained on 2004 
June 04 with an integration time of 111 sec in the NUV 
(2271 A) and FUV (1528 A) bands during the All Sky 
Imaging Survey (AIS) conducted by the Galaxy Evolu- 
tion Explorer (GALEX). The data were extracted from 
MAST (tile AISCHV3_104_01135) as provided by the 
standard GALEX data processing pipeline, and photom- 
etry was performed directly from t hese i mages using pro- 
cedures outlined by Howell et al. (20101. 

Contemporaneously with the X-ray observations by 
XMM-Newton summarized above, the Optical/UV Moni- 
tor Telescope (XMM-OM) obtained observations of Mrk 
266 in the U (3000 - 3900 A), UVW1 (2450-3200 A), 
and UVW2 (1800-2250 A) filters with exposure times 
of 800, 1000 and 1400 sec, respectively. The XMM-OM 
instrument and data characteristics are described by |Ma-| 
son et al. (2001). These images have significantly better 
spatial resolution (~ 2" FWHM) and higher sensitiv- 
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Energy (keV) 
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Figure 5. Chandra X-ray spectra covering 0.35 - 7 keV for the regions illustrated in Figure El (a) Mrk 266 NE ("NE"); (b) Mrk 266 SW 
("SW"); (c) the central region between the nuclei ("Cen"); (d) the Northern Loop ("Loop"); (e) the diffuse emission ("Dif"); (f) and the 
total system emission ("Total"). The horizontal lines indicate the variable- width energy bands, and the vertical lines are la error bars. 
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zero-points in Table 19 of the XMM-Newton Users Hand- 
bookQwhich are 19.1890, 18.5662 and 16.5719 AB mag 
in the U, UVW1 and UVW2 filters, respectively. Figure 
[7] displays the UV imagery of Mrk 266 in the FUV band 
of GALEX and in three filters of XMM-OM. 



2.5. Mid- and Far-Infrared Imaging with Spitzer 

Imaging observations of Mrk 266 were obtained with 
the Spitzer Space Telescope as part of the Cycle 1 Gen- 
eral Observing program for GOALS (PID 3672; PI J. 
Mazzarella ) . The Infrared Array Camera (IRAC, Fazio 



Energy (keV) 



Figure 6. The 0.5 - 10 keV X-ray spectrum of Mrk 266 (in- 
tegrated system) measured with the EPIC instrument onboard 
XMM-Newton. Primary features include the Fe Ka emission line 
at 6.4 keV and an Fe K absorption edge at 7 keV. 



ity than the AIS observations from GALEX, and they 
also considerably extend the UV coverage. The reduced 
and calibrated XMM-OM images, as produced using the 



methods documented by iKuntz et al.l ([2008]) , were ob 



tained from MAST(dataset #005599050T)7Photometric 



measurements (presented in p. 6.2 ) were calibrated using 



et al."|2004[ ) onboard Spitzer was used to image the source 
at 3.6, 4.5, 5.8 and 8.0 fjm (AOR #12305408). Due to 
the presence of bright galactic nuclei, the high dynamic 
range (HDR) mode was utilized in all four bands. The 
HDR frames were comprised of 1-2 second exposures to 
calibrate pixels that may have saturated in the 30-second 
exposures. The IRAC observations involved dithering to 
achieve redundancy to correct for cosmic ray hits and bad 
pixels. The Multiband Imaging Photometer for Spitzer 
(MIPS jRieke et al.|2004| was utilized to image Mrk 266 
at 24, 70 and 160 /im using the Photometry and Super 
Resolution mode (AOR #12349184). Each observation 
consisted of multiple 3-second integrations to insure ad- 
equate redundancy. Mapping cycles were used at each 
wavelength to correct the data for transients in the de- 
tectors. Details of the post-pipeline data processing and 
calibration of these data are given in Mazzarella et al. 

22 http://heasarc.gsfc.nasa.gov/docs/xmm/uhb/XMM_UHB.html 
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Figure 7. UV imagery of Mrk 266. (a) FUV (1528 A) band of 
GALEX; (b) UVW2 (1800-2250 A) band of XMM-OM; (c) UVW1 
(2450-3200 A) band of XMM-OM; (d) U (3000 - 3900 A) band of 
XMM-OM. 

2012 (in preparation). Various Spitzer images and pho- 
tometric measurements for Mrk 266 are incorporated into 
the following sections in conjunction with the analyses of 
data in other spectral regions. 

2.6. Mrk 266 from Radio Through X-Rays 
2.6.1. Multi-Waveband Imagery 

Figure [8] shows the structure of Mrk 266 in nine spec- 
tral regions spanning radio through X-rays using the ob- 
servations discussed above. All images are displayed with 
the same orientation and scale to facilitate comparison. 
In each panel, contou rs of the 6 cm (4.885 G Hz, 0?!3 x 
Of! 4 beam) image from Mazzarella et al. ( 19881) are super- 
posed, primarily for astrometric reference. It is notewor- 
thy that the elongated radio continuum structure with 
double peaks observed between the galaxy nuclei, as well 
as the elongated morphology of the SW galaxy, was con- 
firmed in more se nsitive 8.4 GHz obse rvations with 0?!27 
x 0722 resolution ( |Leipski et aL|2006| ). 



2.6.2. Photometry and Spectral Energy Distributions 

Using the data displayed in Figure [8] and other pass- 
bands discussed above, aperture photometry was per- 
formed in elliptical regions 
pcndix (Table 



aarameterized 
31). However 



in the Ap- 
an alterna- 



11 and Figure 

tive measurement procedure was necessary for the lower 
resolution far-infrared data. The spatial resolution of the 
MIPS images, determined by the size of the Spitzer PSF 
(A/D) are 6", 18", and 40" FWHM at 24, 70 and 160 /im 
respectively. The 24 /im IRAC image (Fig. |8b) clearly 
resolves the two galaxies, and it also reveals oust emis- 
sion from the Northern Loop which is analyzed in { 3.4.3| 
The 10" angular separation of the two galaxies is too 
small to be resolved by MIPS at 70 and 160 /im. How- 
ever, as shown in Figure [8b, the centroid of the 70 /im 



emission clearly indicates that the SW galaxy dominates 
the far-infrared emission from the system. The centroid 
of the 160 /im emission is also consistent with the SW 
galaxy providing the majority of the flux. 

The 70 and 160 /im flux densities for the two g alaxies 
wer e estim ated as follows: The relation derived by |Xu et" 

aDpnol), 



log(Li r ) =log(L 24 )+0.87(±0.03)+0.56(±0.09)xlog(L 8 /L24), (1) 

was used with measured monochromatic luminosities at 
8 and 24 /jm to derive L ir (1 - 1000 /jm) estimates of 
2.3 x 10 u , 6.6 x 10 10 , and 3.4 x 10 11 L Q for Mrk 266 
SW, NE, and the total system respectively. Then, the 
ratios of L7o/L; r = 0.58 and Li6o/L; r = 0.19 measured 
for the total system were applied to estimate the 70 and 
160 /im flux densities for the two galaxies (Table 13). 
The estimated uncertainty in the fluxes derived using this 
procedure is «20%. Comparison of the summation of the 
estimated component flux densities with the measured 
total 70 and 160 /jm flux densities for the total system 
(Table 12 ) indicate s co nsis ten cy within the uncertainties. 



Appendix Tables 12 and 13 present photometry for the 
total system and forthe components (SW, NE, the cen- 
tral region between the nuclei, and the Northern Loop). 
Figure ©[displays the resulting Spectral Energy Distribu- 
tions (SEDs). 

2.7. Mid-Infrared Spectroscopy with Spitzer 

Mid-infrared spectroscopy of Mrk 266 obtain ed in star- 
ing mod e with the Infrared Spectrograph (LRS,| Houck et 
al.|2004 ) onb oard Spitzer h ave been pu blished byjBrandl 
^aUj2006j) , |Dudik et aL] (|2007[ ), and |Bernard-Salas et 
aL~ J 2~t)09~ [TUnfortunately, these observations suffer from 



aperture affects that resulted in missing emission from 
some components in different spectral regions. Appendix 
|5.2| provides clarification of which emission features have 
been measured for each component of Mrk 266. 

In order to overcome the limited spatial extent of the 
staring-mode observations, the IRS spectral map ob- 
tained on 2005 January 10 available in the Spitzer archive 
(program P03269; PI J. Gallimore; AOR #12459264) 
was used to construct spectra for each nucleus and for 
the total system. Spectral mapping with the SL module 
consisted of 13 integrations, each 6 seconds in duration 
and stepped perpendicular to the slit with an interval of 
1"8 (half the slit width). Mapping with the LL mod- 
ule was performed with 5 integrations, each 6 seconds 
in duration and stepped perpendicular to the slit with 
an interval of 4'.'85 (half the slit width). The individ- 
ual calibrated spectra produced by version S15.3 of the 
Spitzer Science Center pipeline were used t o construct a 



spect ral cube with the CUBISM package (Smith et al. 
2007|). In all cases, background observations were assem- 



bled from the non-primary slit and subtracted from the 
data. CUBISM was used to construct images over spe- 
cific wavelength ranges and spectral features, including 
5.5, 10, and 14 /im continuum and continuum-subtracted 
6.2, 7.7, and 11.3 /im PAH, 9.7 /mi H 2 S(3), and [S IV] 
10.5 /i m e mission lines. These images are presented in 
Figure |10[ where the first panel illustrates rectangular 
apertures used in CUBISM to extract 1-d spectra. 

With a pixel size of 5'.'lx5'.'l and a spatial resolution of 
10" on the red end, we cannot use the LL observations 
(14-38 /im) to separate the two galaxies as we can with 
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Figure 8. Imagery of Mrk 266 spanning radio through X-rays. The images are displayed at their native resolution (no PSF matching) 
and registered co-spatially within the uncertainty of the astrometry in each bandpass. Each panel shows the same 32'.'5 X 43'.'7 field of 
view. Labeled in each panel are the bandpass and telescope (upper left), and a scal e bar indicating 5" (3.1 kpc). From the upper left 
to the lower right, each panel shows contours of the 6 cm radio continuum emission jMazzarella et aI7||1988[ l superposed on a grayscale 
representation of the following: (a) 20 cm radio continuum from the VLA; (b) 70 /im Spitzer MIJ-'S; (c) ^4 (im Spitzer MIPS; (d) 8 /im 
Spitzer IRAC; (e) 1.6 /im HST NICMOS; (f) I band (0.81 /im, F814W) HST ACS; (g) B band (0.44 /mi, F435W) HST ACS; (h) U band 
(0.35 Mm) XMM-OM; (i) X-rays (0.4 - 2 keV) Chandra. 



the SL data (5.2-14.5 /xm). Further, the LL map did 
not completely cover the NE nucleus. A "total" SL+LL 
1-d spectrum was extracted that contains most of the 
integrated emission of the system over 5-38 /zm in a re- 
gion 4x4 LL pixels in size (outlined in brown in Figure 
10 1). Another 1-d spectrum was extracted that fully 
covers both galaxies in the 5-14 /im SL spectrum in a re- 
gion 13x13 SL pixels in size, where each pixel is l'.'8xl'.'8. 
The resulting 1-d spectra, representing the sum of the 



and SW nuclei; these regions are shown as red and blue 
rectangles (each covering 5x5 SL pixels) in Figure 10 1 



and the resulting 1-d spectra are presented in Figures 
12b and [121: 



emission from both gala xies , are presented in Figure 11 
(SL+LL) and in Figure 12 1 (SL only), respectively. In 
addition, 1-d extractions from the SL map were con- 
structed that provide a clean separation between the NE 



Due to the very high mid-infrared surface 
Enghtness of the galaxies, with a spatial resolution of 
3'.'6 (two 1'.'8 pixels), even the SL spectral map does not 
enable a reliable separation of the nuclear emission from 
the diffuse emission surrounding the nuclei. 

This is the first time mid-infrared spectra are avail- 
able for each nucleus in Mrk 266 separately, as well as 
the integrated emission from both galaxies. Table [3] lists 
measured line fluxes and equivalent widths for the to- 
tal system and, where possible, the individual galaxies. 
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Figure 9. SEDs for Mrk 266 SW (red), Mrk 266 NE (blue), the Central Region (green) between the nuclei, the Northern Loop 
(purple), and the integrated emission from the Total System (brown). The legend shows the color coding used to identify the various 
components. Most of the photometric measurements are presented here for the first time and some data were selected from NED, as 
listed in Tables [l2"| and [l3| The monochromatic fluxes plotted are computed via log(ff„[JyHz]). Conversions to other units are as follows: 
log(i/f„[Wm- 2 ]) = log(i/f„[JyHz]) - 26.00; logi>L„[W]) = log(i/f„[Wm- 2 ]) + 50.30; log(^L„[L Q ]) = log(i/f„[W]) - 26.58. 



The fluxes were converted from surface brightness units 
(W m -2 sr" 1 ) to W cm~ 2 by multiplying by the ex- 
traction aperture. Most of the spectral features in the 
SL+LL spectru m of the integrated emission from the 
galaxy p air (Fig. [IT]) were m easured using the PAHFIT 
package (Smith et al. 2007), which includes decompo- 
sition of various spectral components and an extinction 
correction based on a fit with a Tg.7p m = 0.30 dust screen, 
corresponding to Ay ~ 5.7 mag. However, since the nu- 



clei were resolved only in the SL data (Figure 12 d and 



12 :), and PAHFIT produces most reliable continuum and 



feature fits when SL and L L are combined, we have used 



a spline- fitting technique (Hony et al. 2001 Peeters et 
aL]|2002 ) to measure the SL features listed in Table 3] for 
the SW and NE galaxies. This further allows for direct 
comparison to other sources in the literature where this 
technique has been employed. Spline-fitting was also uti- 
lized for the PAH EQW measurements to estimate the 
local con tinuum and the sili cate absorption strength. As 
noted by Smith et al. (|2007 1 , spline fitting systematically 
underestimates PAH fluxes compared to full spectral de- 
composition via PAHFIT. The various measurements are 



identified in Table [3] with superscripts p (PAHFIT) and 
s (spline fitting), respectively. 

3. ANALYSIS & DISCUSSION 

Mrk 266 is a complex collection of many different phys- 
ical subsystems, each of which are identified in Fig. [TJ 
and discussed in the following subsections. At the macro 
level, Mrk 266 is clearly a merging system with two dis- 
tinct AGNs coupled with strong spatial differences in cur- 
rent star formation, outflow, merger dynamics and AGN 
heating. But within that framework there are several 
specific questions that we can directly address using our 
multiwavclcngth data set: What are the structural pa- 
rameters and underlying stellar populations of the two 
host galaxies? Is there evidence for dynamic outflow 
from each nucleus? Is the NE LINER a genuine AGN 
or is it powered by star formation? Is the origin of the 
observed diffuse emission tidal or due to galactic super- 
winds? What are the properties of the suspected star 
clusters that are forming in the merging system? Is Mrk 
266 a unique dual AGN system, or are we observing a 
short-lived phase of an evolutionary process that com- 
monly occurs when two massive, dusty disks merge? 
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(a) 5.5 um Continuum (b) 6.2 um PAH (c) 7.7 um PAH 

SL+LL Total 




(g) 11.3 um PAH (h) 14 um Continuum 




Figure 10. Image slices through the 2-d spectral IRS map. Panels (a) - (h) show images in the indicated spectral features (pure line 
or continuum emission). All images have north up and east left, and have the same scale indicated by the bar in panel (a). The colored 
rectangles overlaid on the 5.5 (im image in panel (a) indic ate the regions of the 2-d spectral cube used to extract 1-d spectra. Brown box: 
"total" emission in the SL+LL modules plo tted in Figure [TT] (Mrk 266NE was not fully covered by the LL aperture); oran ge b ox: "total" 
emission in the SL module plotted in Figure |l2^ ,; red box: emissio n fro m the NE nucleus in the SL module plotted in Figure [12p ; blue box: 
emission from the SW nucleus in the SL module plotted in Figure [12} :. 
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Figure 11. 1-d spectrum of Mrk 266 extracted from a low- 
resolution spectral map constructed from the SL and LL modules 
of the IRS. The region use d to produce this spectrum is shown as 
a brown box in Figure [To] This is our best attempt to represent 
the integrated emission tor both galaxies across the full 5-38 (im 
range covered by the low-resolution modules of the IRS, with the 
caveats noted in the text. The colors in this spectrum represent 
the various or der s in the IRS instrument and are not related to the 
colors in Fig. |10| 

3.1. Interacting/ Merging Galaxies Revealed 
3.1.1. Morphology and Bulge/Disk/Bar Decomposition 

High resolution HST imagery (Fig. [I]) has revealed 
a wealth of morphological structures within Mrk 266. 
At optical wavelengths (B and I) there are two distinct 
galaxies within a diffuse and highly asymmetric envelope. 
Patchy dust obscuration and clumpy star formation re- 
gions are evident on a variety of scales. The underlying 
structure of the two galaxies is best revealed in the NIC- 
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Figure 12. The spatial resolution of the SL observations permit 
a clean separation of the two galaxies that is not possible with 
the LL observations. Each panel is a 1-d extraction from the SL 
spectral map presented in Figure |10| (a) orange region covering 
the entire region; (b) red region containing Mrk 266 NE; (c) blue 
region containing Mrk 266 SW. 
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Table 3 Table 4 

Measurements from Low-Res Spectral Mapping GALFIT Results 



Feature Flux (1 X 1(T 20 W cm' 2 ) EQW (fim) Band C(n) § r c | PA° mag AB 

(1) (2) (3) (1) ^2) (3) (4) (5) (6) (7) (8) 



Mrk 266 Total 3 



Mrk 266 SW a 



6.2^im PAH 


40.2 


(±1.8)p 


H 2 S(5) 6.91 Mm 


0.66 


(±0.50) p 


[Ar II] 6.99Atm 


2.34 


(±0.72)P 


7.7/jm PAH Complex 


163.0 


(±6.1)p 


8.3/im PAH 


13.7 


(±1.8) p 


8.6Mm PAH 


26.8 


(±1.5) p 


[Ar III] 9.02^m 


0.80 


(±0.49)p 


H 2 S(3) 9.65Atm 


1.42 


(±0.38)P 


[S IV] 10.51^m 


1.59 


(±0.29)p 


11.3/im PAH Complex 


31.9 


(±1.2)p 


12.0/im PAH 


11.85 


(±0.83)P 


H2 S{2) 12.27Atm 


0.89 


(±0.15)p 
(±0.73)p 


12.6/im PAH Complex 


23.78 


[Ne II] 12.81/Ltm 


6.78 


(±0.18)p 


13.6/iin PAH 


8.86 


(±0.70)p 


14.2/xm PAH 


2.80 


(±0.52)p 


[Ne III] 15.56A*m 


3.57 


(±0.24)p 


16.4/mi PAH 


2.08 


(±0.38)p 


H 2 S(l) 17.04^m 


1.61 


(±0.19)p 


17^m PAH Complex 


19.4 


(±1.7)p 


17.4/Ltm PAH 


1.01 


(±0.36)p 


[S III] 18.71/im 


4.16 


(±0.28)p 


[O IV] 25.9Atm 


4.52 


(±0.14)p 


[S III] 33.51/im 


3.71 


(±0.28)p 


[Si II] 34.86Aim 


5.69 


(±0.39)p 



0.47(±0.03) s 



46(±0.02) s 
15(±0.01) s 

51(±0.01) s 



Mrk 266 SW b 



6.2^tm PAH 


17.78 


(±0.35) s 


0.67(±0.03) s 


7.7/rni PAH Complex 


30.46 


(±0.35) s 
(±0.14) s 


0.54(±0.01) s 
0.17(±0.01) a 


8.6fj,m PAH 


5.45 


[S IV] 10.5lAtm 


0.95 


(±0.13) s 




11.3/im PAH Complex 


10.376 (±0.088) s 


0.73(±0.01) s 


[Ne II] 12.81/mi 


6.15 


(±0.99) s 





Mrk 266 NE b 



6.2^m PAH 


2.418 (±0.019) s 


0.4 (±0.1) s 


7.7/im PAH Complex 


2.86 (±0.28) s 
0.69 (±0.11) s 


0.25(±0.03) s 
0.10(±0.02) s 


8.6fj,m PAH 


H 2 S(3) 9.65/im 


0.52 (±0.14) s 




11.3Atm PAH Complex 


2.170 (±0.078) s 


0.34(±0.02) s 


H2 S(2) \2.21y,m 


0.308 (±0.063) s 




[Ne II] 12.81/mi 


1.84 (±0.29) s 





a Spectral measurements for the integrated emission from both galax- 
ies in Mrk 266. Features in the 5.2-14.5 fim (SL) and 14.5-38 /im (LL) 
regions were measured from the 1-d spectrum plotted in Figure |11| 
this region is 4x4 LL pixels in size, where each pixel is 5'. 1 x 5. 1, as 
outlined in brown in Figure [lO[ a) . 

k Spectral features measured from the SL spectra of Mrk 266 SW and 
NE as plotted in Figure [l2^ and |12| j. These 1-d spectra were extracted 
from the blue and red regions of the spectral map shown in Figure [l0| 
which have dimensions of 5x5 SL pixels, where each pixel is 1. 8 x 1. 8. 
p These line fluxes were measured using PAHFIT and include an ex- 
tinction correction based on a fit with a tq _j^ m — 0.30 dust screen, 
which corresponds to Ay m 5.7 mag. 

s EQW measurements for the PAH features were measured using a 
spline-fitting technique for estimating the local continuum. Spline fit- 
ting was also used to estimate fluxes of features in the SL spectra of 
the individual galaxies which could not be measured with PAHFIT. Sec 
text for details. 



MOS H band image. The SW galaxy (Fig. n3j) is a barred 
spiral with two arms and predominant dust teatures, and 
the estimated Hubble type is SBb (pec) . The NE galaxy 
(Fig. 14) has no discernible spiral arms and appears to 
be Hubble type SO or SO/a (pec). 
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Note. — Fitting parameters and results of modeling the HST data 
using GALFIT. Column (1): HST Filter: B = F435W (0.44/im); 
I = F814W (0.81/W); H = F160W (1.6A*m). Column (2): Per- 
centage of the total flux of the GALFIT model of the galaxy pair 
accounted for by the indicated component of the model image 
at 0.44, 0.81 or 1.6 fim. Column (3): Model components coded 
as follows: B=Bulge (Sersic index n=4); D=Di sk (Sersic in dex 
n=l); b=bar; P=Point source based on Tiny Tim |Krist| ( |1993|l; T 
= Total flux from all GALFIT components for the galaxy. Col- 
umn (4): Bulge-to-Disk ratio. Column (5): Effective radius of 
the Sersic fit (kpc units). Column (6): Axis ratio of the compo- 
nent fit. Column (7): Position angle of the component fit (de- 
grees east of north). Column (8): Apparent AB magnitude of 
the sum of the GALFIT model components for each galaxy, where 
f v (Jy) = iQQ-4(8.9-mAB). Conversions between the AB and V ega 
systems (http://www.stsci.edu/hst/acs/analysis/zeropoints/) are 
as follows: BAB - tiVega = -0.11;1aB - IVega = U.42;Mab — 
Hvcga = 1.32. 

a The 1.6 Aim imago of Mrk 266 S W is modeled as the sum of Pointsource + Bar + 
Disk (n = 1) -|- Bulge (n = 4) components. The 0.44 and 0.81 rim images are too 
confused by extinction and H II regions to derive reliable models of the underlying 
stars; simple Bulge (n = 4) models were fit at these wavelengths only to examine 

the residual images for fine structure in the star-forming regions (Fig. [ToV 

b 1 1 

The 1.6 /rrn i ma ge of Mrk 266 NE is modeled as the sum of Pointsource + 

Bulge (n = 4) components; there arc no significant Bar or Disk components de- 
galaxy. 

x (apparent magnitude) accounted for in the GALFIT model for the 
each HST passband. 
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C The total fl 
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In the extremities of the system are three notable fea- 
tures at optical wavelengths (Fig. [T]): (1) Extending 
~25" (15 kpc) to the north is a fragmented, filamen- 
tary feature known as th e No rthern Loop; this structure 
is examined in detail in { 3.4 (2) Faint emission extend- 
ing -60" (36 kpc) to the south and -40" (24 kpc) to 
the south-east from the center of the system appears to 
be tidal debris. The deep B+V+I image (Fig. [2]) re- 
veals much fainter, asymmetric emission spanning sal03 
kpc (2.'9). The peculiar morphology and vast extent of 
this emission are consistent with numerical simulations 
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Figure 13. A close-up of HST imagery of Mrk 266 SW showing 
a color composite of ACS B+I bands (a); NICMOS H band (b); 
ACS B band (c); and ACS I band (d). The images are displayed 
with an emphasis on the high surface brightness emission near the 
nucleus. The spatial resolution is very similar in all 4 bands, but 
there has been no PSF matching. The scale bar indicates 2". 




Figure 14. HS T im agery of Mrk 266 NE displaying the same 
bands as in Figure [13] 



Mi- 



of tidal debris created du ring a major merger (e.g. 
hos|1999||C~ox et al.|2006| . (3) Approximately 25" to'tne 
south-west are two objects that could be either knots 
near the end of a tidal tail or background galaxies. The 
latter interpretation is favored due to what appear to 
be compact nuclei and spiral arms. If these are back- 
ground galaxies with similar physical size and separation 
as Mrk 266, their projected diameters («3") and sepa- 
ration (w2") imply they are ^5 times further away than 
Mrk 266 (distance ^650 Mpc, z ~ 0.15); redshift mea- 
surements are nee ded for confirmat ion. 

Using GALFIT ( |Peng et al.|2002[) , quantitative isopho- 
tal analysis was performed on both component galaxies 
in the three HST bands. At 1.6 /zm (H band) where 
dust obscuration is minimal and thus permits the most 
reliable decomposition of the underlying galaxies, vari- 
ous trial fits involving possible point source, bulge, disk 
and/or bar components were made; potential nuclei were 
modeled using a point spread func tion gener ated from 
the Tiny Tim package ( |Krist|[l99 3 p) Figure [15] shows 



the components of the best-fitting models in the H band. 
Due to extensive, patchy dust obscuration and complex 
star-forming regions, only models with possible bulge and 
disk components were considered in the B and I bands. 
Subtraction of the best-fitting models from the direct 
images are useful to elucidate further fine structures, as 
shown visually for all three bands in Figure [16] The 
results of this structural analysis at H band are the fol- 
lowing: a) The SW galaxy has no detected point source 
emission; its nucleus is therefore obscured by dust even 
at 1.6 fim. b) The NE galaxy exhibits no significant disk 
component and therefore can be modeled using only a 
point source plus bulge profile, c) Approximately 10% of 
the total H-band light is in the form of diffuse emission 
which is likely tidal in origin (Fig. 15 1). In the optical 
bands the residual images reveal the following: a) The 
SW galaxy has many knots of emission which are likely 
star clusters or a ssociations; these features are analyzed 
in detail in |3.6| b) The largest features in the B-band 
residual image appear to be giant star-forming regions 
(superassociations) situated on both ends of the stellar 
bar. This phenomenon, in which superassociations occur 
on both sides of stellar bars much more frequently than 
they occur on a single side, was pointed out by Magtcsian 



& Khachikian (1991). c) The NE galaxy has numerous 
radial, asymmetric filaments (most pronounced in the B 
band) indicative of bi-conical out flow fr om the AGN; this 
region is analyzed more fully in §3.2.4| Parameters from 
the GALFIT modeling are listed in Table Ej Adopting a 
distance modulus of 35.55 mag (see j ]1.3[ ), the apparent H 
magnitudes (Vega system) of Mrk 266 SW and NE (Table 
4 ) correspond to remarkably similar absolute magnitudes 
of M H (SW) = -24.32 and M H (NE) = -24.16. 



3.1.2. Galaxy Luminosities and Derived Stellar Masses 

The H-band (1.6 fim) luminosity function of galaxies in 
the Coma Cluster has been appr oximated as a Schechte r 
function with = —23.9 mag ( |de Propris et al.]|1998[ ). 
Therefore, Mrk 266 SW and NE are «50% and 30% more 



luminous than the typical L^j galaxy in the local universe. 
Applying the mean mass-to-light ratio at H-band of 4.6 
(M & /Lq) for spiral galaxies (i ncluding early types SO 
and SO/a; Gavazzi et al. (1996)) results in total galaxy 
mass estimates of 6.3 x 10*M© (SW) and 5.3 x 10 10 M Q 
(NE). Estimates of the stellar mass in old stars derived 
from the mean ratio of M» /i^ L t/ (2.2jnm) = 5.6( ±1.5) 
based on model predictions by Lacey et al. ( 2008[ ) pro- 
duces similar results of M* = 6.1 x 10 iu M Q (SW) and 
M* = 4.4 x 10 10 M (NE)QSED model fitting can im- 
prove on st ellar masses est imated from monochromatic 
luminosities U et al.| (2011 1. Next we turn attention to 
the nuclei and their circumnuclear regions. 



The model curve for z : 



http://www.stsci.edu/software/tinytim/tinytim.html 



: galaxies in Figure 13c of |Lacey et| 
|al,| l |2008|l , adjusted to h = 0.70, corresponds to a range of 4.5 - 7.0 
m lVU/^L 1/ (2. 2/^m). A significantly lower ratio of M«/L2.2jj m ~ 1 
was found by |Arnouts et al.| | |2007| and |Cole et aT| ( |2001^ 
nature of this discrepancy is beyond the scope of this article 
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Figure 15. Results of GALFIT modeling of the 1.6 NICMOS 
image of Mrk 266. The top row shows (a) the direct image, (b) 
the GALFIT model fit, and (c) the residual image (direct image 
minus model fit). The center row shows the components of the 
model for the SW galaxy: (d) the disk, bulge (e), and stellar bar 
(f). The bottom row shows the components of the model for the 
NE galaxy: (g) the bulge and (h) the point source modeled with 
Tiny Tim. The scale bar indicates 10" (6 kpc). 

3.2. The Nuclei and Circumnuclear Regions 
3.2.1. Optical and Infrared Spectral Diagnostics 

Using optical emission-line diagnostics, the NE nucleus 
has been classified as a LINER and the SW nucleus as a 



Seyfert 2 (e.g., Hutchings et al.|1988 Mazzarella & Boro- 



son 



1993 Ishig aki et al.||2000| )l At optical wavelengths, 
the N E nucleus has brighter continu um emission and less 
extinction than the SW nucleus. Yuan et al. (20101 



classified Mrk 266 NE as a composite AGN/starburst 
nucleus, and they concluded that the NE and SW nu- 
clei have similar relative AGN and starburst contribu- 
tions within the LINER and Seyfert 2 branches of the 
[O III]5007/iJ^ vs. [O I]6300/Ha diagnostic diagramF 
Since the I RS a pertures encompass an area of 5.4 x o~4~ 



Figure 16. Residuals resulting from subtraction of GALFIT 
models parameterized in Table[4]from the direct HST images. The 
left column shows results for the NE galaxy at 0.44 fim (a), 0.81 
fim (b) and 1.6 (im (c). The right column shows results for the SW 
galaxy at 0.44 fim (d), 0.81 fim (e) and 1.6 fim (f). Black (white) 
regions represent positive (negative) flux with respect to the image 
model fits. The tick marks are separated by 0"5. 

(Table[3|) indicate the prominence of star- formation (e.g., 
Armus et al.|2007[ ). The NE source has lower PAH fluxes 
and equivalent widths, indicating excess warm dust emis- 
sion compared to the SW source. 

The 6.2 urn PAH EQWs and 9.7 /im silicate absorption 
strength^" 
vertex oft 



kpc (Fig. 10 1, whereas published optical diagnostics 
were obtained through smaller apertures of ~1 x 1 kpc 2 , 
nuclear emission is more diluted by extranuclear emis- 
sion in the mid-IR spectra than in the optical spectra. 
However, since the mid-IR spectra penetrate much more 
dust than the optical spectra, if there is sufficient cen- 
trally concentrated dust heated by an embedded AGN, 
this emission can potentially overpower the extended star 
formation (as in many ULIRGs). The IRS spectra of the 
two nuclei (Fig. 12 ) show them to be similarly dominated 
by PAH emission, and the large PAH equivalent widths 

25 The identifications of the NE and SW nucl ei of Mrk 266 were 
mis takenly swapped in th e tables and figures o f|Kim et al.||l995| 
and |Veilleux et al.| l |1995| l, and thus also in Table 2 ot |Yuan et al. | 
(201U). Tne identifications are corrected here. 



place the SW galaxy on the bottom right 
re silicate depth versus PAH EQW plot of 
Spoon et al. (2007), a region dominated by galaxies with 
high star formation rates; the NE galaxy (with similar Si 
strength but a lower 6.2 ^m EQW of 0.4 /im compared 
to 0.67 /mi for the SW galaxy) is in a region populated 
with AGN and ULIRGs having substantial dust heating 
by non-stellar radiation. The observed silicate strength 
in the SW galaxy (—0.60 ±0.06) is intermediate between 
the values observed for the face-on galaxy NGC 7714 
and the edge-on galaxy M82, which is consistent with 
the derived disk axial ratio of 42° for Mrk 266 SW 



26 The empirical method of |Spoon~e t a l. (2 007|l was applied to 
the SL spectrum (5-14 fim region; Fig. |T2| | to measure silicate 
strengths of -0.62 (±0.06), -0.69 (±0.12), and - 0.66 (±0.07) for 
the SW nucleus, the NE nucleus, and the total system, respectively. 
A similar absorption strength, —0.60 (±0.06), was meas ured using 
the 5-40 fim SL±LL data for the total system (Fig. 111}. 
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Table 5 

Estimated AGN Contributions to the Mid-Infrared and Bolometric Luminosity 



Diagnostic SW Galaxy NE Galaxy Total System 

Value Sgn i±agn Value Sgn Value Sgn 

L Diag ^Bol ^Diag L Bol L Diag 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


w 


(10) 


(1) [O IV] 25.89Atm/[Afe II] 12.81/mi a 


0.93 


23% 


73% 








0.67 


17% 


64% 


(2) [Ne V] 14.32^m/[JVe II] 12.8lAim a 


0.14 


11% 


53% 














(3) PAH 6.2/tm EQW [fim] h 


0.67 


0% 




0.4 


30% 




0.47 


10% 




(4) PAH 7.7^m EQW H b 


0.54 


69% 


43% 


0.25 


79% 


56% 


0.46 


71% 


46% 


(5) f(PAH 6.2 M m)// 5 . 5At m VS. /(15Aun)// 5 . 5M m b 


1.6,8.6 


47% 


9% 


0.56,6.4 


81% 


33% 


2.0,2.9 


52% 


11% 


(6) /430)//,(15) c 


3.9 


89% 


82% 


3.3 


92% 


87% 


4.8 


85% 


76% 


Mean d 




40% 


52% 




71% 


59% 




47% 


49% 


Std. Dev. 




35% 


29% 




28% 


27% 




33% 


28% 



Note. Column (1): Mid-infrared spectral diagnostic used to estimate the fractional contribution of radiation from an AGN. 

Columns (2)-(4): For the SW galaxy, the value of the diagnostic followed by the percentage of AGN contribution to the diagnostic 
luminosity and the e stima ted pe rcentage of AGN contribution to L^oi- The latter uses the zero-points and bolometric corrections from 
|Veilleux et aDj2009} ; see g3.2.1| Columns (5)-(7): The same as columns (2)-(4), but for the NE galaxy. Columns (8)-(10): The same as 
columns ('A)-(i), but for the total system. 

a Based on published high-resolution, staring-mode data (Appendix Table [14} . The NE nucleus was not fully covered in these observations. 
b Based on new low-resolution measurements from the IRS spectral map (Table [31. 

c The continuum flux density ratio estimated via interpolation in the SED of the global system, /„(30/im,)//„(15/jm) = 4.8, is within 8% 
of the value /„(30/im)/ /„(15/im) = 5.4 computed from the IRS (SL+LL) spectrum. Given this close agreement, since the 15 and 30 (im 
fluxes could not be measured for the individual galaxies with the IRS (because they were unresolved by t he L L module), the flux ratios 
estimated by interpolation from the broad-band photometry in the SEDs for the individual galaxies (Table [l3| were substituted. 
cl The mean Lagn/LboI values omit the 6.2 /im PAH EQW because there is no bolometric correction available for this diagnostic. 



(Table [4]). The observed silicate absorption strengths 
correspond to Nh = 2.2 x 10 22 and 2.4 x 10 22 cm" 2 for the 
SW and NE nuclei, respectively; these values result from 
applying the relationsh ips Nh = T9.7 um x 3.5 x 10 22 cm~ 2 
and A 9 . 7/im /A v = 0.06 (|Roche & Aitken|1984|). However, 



for clumpy dust configurations expected m AGN tori, 
emission and absorption cannot be disentangled along 
the line of sight and therefore these are lower limits to 
the true silicate column to the central sources. 

Various empirical diagnostics have been published to 
quantify the relative contributio ns to dust heating from 
AGN and starburst sources (e.g., Armus et al.|2007 ). We 



focus primarily on the suite of methods appl ied to local 
ULIRGs and QSOs by |Veilleux et al.| ( |2009| ), including 



application of their bolometric corrections to express the 
estimated AGN dust heating contributions as a fraction 
of Lboi- All methods assume that the line emitting re- 
gions are subject to the same amount of extinction within 
the aperture, which may be questionable, especially for 
Mrk 266 SW. Since the NE galaxy was not observed in 
the high-resolution IRS mode, methods that involve the 
fine structure lines can only be applied to the SW galaxy. 
The results are summarized in Table [U For the SW 
galaxy, the [O IV] /[Ne II] ratio indicates an overall AGN 
contribution to L|-, i of 73%, whereas the [Ne V]/[Ne II] 
ratio suggests 53%. The 6.2 /zm EQW indicates that the 
mid-IR luminosity of the SW component is due entirely 
to a starburst, and a 30% AGN contribution is inferred 
for the NE galaxy. This diagnostic suggests a 10% AGN 
contribution to the mid-IR luminosity of the total sys- 
tem^] The AGN fractions estimated from the 6.2 /im 
EQW are substantially lower than those obtained using 
the EQW of the 7.7 /mi PAH feature (see Table |]|. This 
discrepancy may be due to greater difficulty measuring 

27 The CAFE SED fitting package provides an independent, 
model-dependent estima te of the AGN contribution to the total 
infrared luminosity. See §3.51 



the 7.7 /im PAH feature because of confusion from sur- 
rounding features (8.6 /mi PAH and silicate absorption) 
and complications in defining the proper continuum level. 

The continuum shape over 3-16 /mi provides another 
diagnostic. The physical basis is that dust grains lo- 
cated near an AGN can reach temperatures as high as 
500-1000 K (before they sublimate), and such high dust 
temperatures would completely dominate the continuum 
emission compared to du st heated by a starbu rst. This 
m ethod was developed by Laurent et al. ( 2000[ ), adapted 
by Peeters et al. ( 2004 ) for interpretation of mid-IR spec- 
tra from the Infrared Space Observa tory (ISO), and mod - 
ified further for Spitzer IRS data by Armus et al. (2007). 
Continuum fluxes in the relevant pass bands were com- 
bined with the 6.2 /im PAH fluxes, and the formulation of 
VeilleuxetaLl (2009) was used to compute AGN fractions 



before and after bolometric corrections. The results sug- 
gest Lagn /Lboi is 33% for the NE galaxy, but only about 
10% for the SW galaxy and the total system. The contin- 
uum flux density ratio f 1/ (3 0/um)/f t/ (15/zm) provi des our 



2009[ ) found 
All-tree and 



final empirical diagnostic; |Veilleux et al. 
this to be an effective substitute tor the 
silicate-free MIR/FIR ratio. 

On average, with a standard deviation of nearly 30%, 
these diagnostics indicate approximately equal contribu- 
tions of AGN and starburst heating of the dust in both 
galaxies (Table [5]) . This is similar to the average AGN 
contribution to Tboi found for ULIRGs in the QUEST 



sample wit h Seyfcrt 2 spectral classifications (Veilleux 



et al. 2009). The variation in these diagnostics may be 
due, in part, to the different optical depths probed by 
different wavelengths. 

Finally, there is a significant difference in the PAH 
pro pertie s of the galaxies. The modeling of Draine <£] 
Li (12001 ) combined with the ratios of the PAH features 



(Tii.3/Wf7.7Mm and i 6 . 2^/^7.7^) indicate a dominance 
of relatively warm, ionized PAHs in the SW component 
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and cooler, neutral PAHs in the NE component. The 
PAH line flux ratios for the total system reflect the fact 
that the SW galaxy dominates the total mid- and far- 
infrared emission from Mrk 266. The iii.^^m/h.itira PAH 
ratio in Mrk 266 NE (0.76) is substantially larger than 
the ratio in the SW galaxy (0.34). The former is similar 
to the ratios observed in other galaxies containing AGNs, 
and the latter is near the mean value (0.25) observed in 



galax ies with pure starburst (H II) spectra (|Smith et al 
20071. Following the inferences proposed by Smith et 



al., this suggests that the hard radiation field surround- 
ing the partially obscured AGN in the NE galaxy may 
be preferentially destroying PAH molecules small enough 
to emit at 7.7 //m, thus increasing the 11.3/xm/7.7/xm 
PAH ratio compared to the SW nucleus which has a more 
heavily obscured AGN and a higher star formation rate. 

3.2.2. Warm Molecular Gas 



In the NE galaxy, the IRS spectral maps (Figs. 10 
and 12 ) reveal bright rotational lines of molecular hy- 
drogen, specifically H2 S(3) 9.7/im and H2 5(2) 12.3/im. 
In the SW galaxy, the H 2 5(3) 9.7/xm emission is much 
weaker, and any H2 5(2) 12.3 /xm emission that may be 
present intrinsically is overshadowed by the blue wing of 
the strong 12.6 /im PAH complex. In contrast, the SW 
galaxy has bright [S IV] 10.51 /j,m emission, whereas the 
NE nucleus is undetected in this feature. Give that the 
NE galaxy has 8 and 24 /xm flux densities that are 1.8 
and 5.7 times fainter than the SW galaxy (see Table 13 
and Fig. [9]), it is surprising that the H2 line emission 
in Mrk 26o NE is substantially more luminous than in 
Mrk 266 SW. This suggests that the majority of the de- 
tected H 2 line emission in the NE galaxy has a different 
origin than in the SW galaxy. A prime candidate for the 
enhanced H2 emission in Mrk 266 NE is shock excitation. 



et al 



Ana l ysis of near-infrared spectra of Mrk 266 led Davies 
] ( |2000p to conclude that fast C 
km s _i ) are responsible for «70% 



shocks (v s 
of the H 2 



40 km s i ) are responsible for «70% of the H 2 1- 
S(l) emission in the two nuclei, with fast J shocks 
(v s ~ 300 km s -1 ) and UV fluorescence responsible for 
the remaining «30%. Davies et al. demonstrated that 
the powerful radio continuum and 1-0 S(l) emissions ob- 
served in Mrk 266 NE cannot be explained with pure 
starburst models, and they are likely generated predom- 
inantly by shocks. Furthermore, the presence of 1-0 S(l) 
emission extended out to ~1 kpc from the NE nucleus, 
with a large equivalent width and low S(l) /B7 ratio that 
rule out star formation, provides additional evidence that 
the H 2 emission is excited by spatially extended shocks. 
Davies et al. speculated that Mrk 266 NE may itself be 
a very close pair of galactic nuclei with shocked molecu- 
lar gas settled between them, as observed in NGC 6240. 
Our interpretation of the latest data from Spitzer, HST 
and Chandra favors an AGN-driven o utflow as the origin 
of strong shocks in Mrk 266 NE (see j|3.2.4|. 



3.2.3. AGN Black Hole Masses 



Marconi & Hunt (2003) established an empirical re- 
lation between H-band bulge luminosity and nuclear 
black hole mass. Combining measurements of the to- 
tal H-band absolute magnitudes with the GALFIT im- 
age decomposition results showing that 78% and 99% 
of the total H-band light is in the bulge components 



of the SW and NE galaxies ((3.1.1), the calibration of 
Marconi & Hunt (Group 1 Galaxies) yields black hole 
masses of 2.3 x 10 8 (SW) and 2.6 x 10 8 M (NE). 
An independent diagnostic is available from the high- 
ionization, mid-IR emission lines detected i n Mr k 266 
SW. 



Using the approach of Dasyra et al 
[Ne V] 14.32iim ([O IV] 25.89/mi) emission 



Q2008D , the 
ine luminos- 

ity of 1.58 (10.5) x 10 41 erg s _1 computed from published 



high-resolution IRS measurements (Appendix Table 14 1 
corresponds to M BH = 1.7 (3.2) x 10 8 M Q for the SW nu- 
cleus, which is consistent with the IR photometry. The 
inferred black hole masses are therefore in the typical 
range known to power other luminous AGNs. 

3.2.4. Evidence for Nuclear Outflows and Shocks 

In the previous section the radiation field heating the 
dust in the two nuclear environments was established. In 
this section evidence is presented for dynamical outflows 
from the nuclei. Within a radius of 2" (1.2 kpc) around 



the NE nucleus (Figs. 14:,d and 16 1), the HST optical 



imaging reveals structure that is strongly suggestive of a 
bi-conical outflow with patchy dust extinction. Evidence 
for this interpretation rather than spiral arms is as fol- 
lows: (a) Multiple filamentary structures extend radially 
and asymmetrically from the nucleus rather than in a 
symmetric spiral pattern, (b) There is the lack of im- 
proved definition in any spiral arm pattern in the I-band 
image relat ive to B-band. (c) Photometric deconvolu- 
tion (j |3.1.1 1 indicates the NE galaxy is bulge-dominated. 
There are also conspicuous knots of concentrated emis- 
sion along a concave arc located 0.4 W of the nucleus. 
Table [6] lists photometric measurements of the 3 most 
prominent optical knots. The individual knot luminosi- 
ties oi Lb ~ 3x lO 8 -L0are substantially high er than other 
luminous star clusters in Mrk 266 (see £ 3.6 Fig. 27 1, and 
their orientation in an arc embedded within a triangular 
bow-like structure suggests they may be energized by an 
AGN -driven outflow in a situation similar to M51 (Cecil 



1988 1 rather than by unusually massive star clusters. 

comparison of the B-band emission with pub- 
lished high-resolution 18 cm radio continuum observa- 



tions (iThean et al.||2001| is shown in Figure [IT] The 18 
cm MERLIN map exhibits limb-brightening at the loca- 



map 

tion of the optical arc and knots 0'.'4 (240 pc) west of the 
nucleus, and a secondary radio continuum peak is located 
0'.'6 (360 pc) NE of the nucleus. An axis at position angle 
56° intersects the radio/optical nucleus, the radio emis- 
sion 0f!6 NE of the nucleus, and B-band knots situated 
/ .'4 and Of! 9 SW of the nucleus. Bi-conical AGN outflows 
with aligned radio continuum and optical emission-line 
radiation with dimensions of <1 kpc are commonly ob- 
served in Seyfert galaxies and in some LINERS, espe- 
cially those that are nearly face-on such as NGC 1068 
and M51. (See, for example, the review of Veilleux et 



al. (2005).) The optical morphology of the JNE nuclear 
environment is remarkably similar to the radiative bow 
shock observed 5" (230 pc) south of the nucleus of M 51, 



which is also a LINER (|Cecil][1988 Ford et al.|1985|. 



The combined data suggest that the AGIN m Mrk 266 
NE is ejecting radio-emitting plasma along an axis with 
PA m 56° , which is also aligned with a bi-conical ioniza- 
tion cone that extends to a radius of 1.2 kpc from the 
nucleus. The 6 cm radio power (4.2 x 10 22 W Hz -1 ) and 
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Figure 17. Comparison of new 0.44 (HST ACS) im agery with 
published 18 cm radio continuum (MERLIN) emission jThean et| 
|al,|200l" Fig. 1) in the circumnuclear region of Mrk 266 INHi. lop 
(a.): HST ACS B-band image in grayscale; the same B-band data 
are contoured (blue) and overlaid here and in the following two 
panels. Center (b): 18 cm radio continuum interferometry with 
natural weighting and spatial resolution 0'.'5 x 0'.'3 HPBW (thick 
black contours) overlaid with ACS B-band (thin blue contours). 
Bottom (c): 18 cm radio continuum interferometry with uniform 
weighting and spatial resolution 0'.'2 x O'.'l HPBW (thick black 
contours) overlaid with ACS B-band (thi n blue contours). The 
plus signs in panels (b) and (c) are from ( |Thean et al.|[2001| see 
their Fig. 1) and they mark the optical nuclear positions, wnere 
the marker extent indicates the 2cr positional uncertainty. These 
positions are in good agreement with the astrometry we derived 
for the ACS images. The scale bar indicates 1", and the field of 
view is 6'.'0 x 4"3 for all three panels. 

spectral index (a 20 , 6i2cm = 0.61 ± 0.04) of Mrk 266 NE 
are consistent with op tically thin synchrotron emission 



(Mazzarella et al. 1988 1, and these parameters are s imilar 



et al. 1996). The radio jet appears to be running into 
dense material concentrated 074 (240 pc) SW from the 
nucleus, but it is able to es cape to a greater distance in 
the opposite direction (Fig. 17). This hypothesis is con- 
siste nt wi th the velocity field measured for the molecular 
gas (j |3.7[ ) which shows a rotating disk aligned roughly or- 
thogonal to this putative ionizat ion cone. The B-I color 
map from the HST data (Figure f26j) reveals that dust is 
distributed primarily along the NS direction, which again 
is roughly orthogonal to the AGN jet/outflow axis. Re- 
gions inside a radius of ^2" from the nucleus of Mrk 266 
NE (Fig. 26) are as blue as the Northern Loop, which 
is known to be d ominated by bright [O III] A5007 and 
Ha line emission flHutchings et aL |1988 Ishigaki et al 



2000 



28 



Davies et al.| ( |2000 ) provided independent evi- 



dence for shock-heated gas m this region in the form of 
the radial profile of the 1-0 S(l) H 2 2.12 /im emission line. 
Furthermore, the B r7 emission-line image presented by 



Davies et al. (2000) shows a bimodal structure oriented 
at .FA « 50 centered on the NE nucleus, providing ad- 
ditional support for our hypothesis. While most of this 
evidence is indirect, its combined weight favors our in- 
terpretation that the extranuclear emission around Mrk 
266 NE is dominated by an AGN-powered outflow that 
is generating spatially extend shocks. 

Similar events may be occurring in the SW galaxy. The 
HST images show numerous compact structures over a 
1 kpc size region that share the same NS orientation 
as the prominent stellar bar seen in the H-band im- 
age. It is highly likely that thes e compact structures 
are yo ung stellar clusters (see ^3.6). Observations at 1.6 
GHz (iThean et al. 1120011) show sub-arsccond structure 



aligned IN W-SF with respect to the SW nucleus, simila r 



to previous 5 Ghz measurements (Mazzarella et al 



1988) 



Comparing the new HST images with the extended ra- 
dio emission maps suggests the circumnuclear emission 
in this case is dominated by star formation. An AGN 
outflow in the SW galaxy may be present but is difficult 
to detect due to its relatively high inclination (sa 42°). 
It is very difficult to disentangle starburst-driven winds 
from AGN-driven winds in galaxie s containing b oth a 



Seyfe rt 2 and a powerful starburst |Veilleux et ak] (e.g. 
2005 1 . Although evidence for an AGN-driven outflow in 



the circumnuclear region of Mrk 266 SW is weak, there 
is ample evidence for a starburst-driven superwind on a 
larger scale based on the filamentary Ha and [O III] line 
emission that comprise the Northern Loop, and from the 
fainter emission-line filaments labeled "eastern arm-like 



region " and "west knot" in Figure 2 of Ishigaki et al. 
(2000), all of which appear to connect or align with the 
SW nucleus. Properties of the kpc -sca le outflow revealed 
by the new data are discussed in |3.4| 

3.2.5. Hard X-Ray Properties of the Nuclei 

Basic prope rties of the Chandra data for Mrk 266 were 
presented by IBrassington et all (120071) . We have per- 



to other luminous AGNs in the local universe (e.g., Rush 



28 The [0 III] A5007 line emission around Mrk 266 NE jHutchingsl 
|et al.|1988| l is elongated along the sam e po sition angle as the HS 1 
B-band image contoured here in Figure [l7| This provides evidence 
that the HST B-band image is dominated by emission-line gas in 
this region. The passband of the F435W filter omits the [O III] 
lines, but it includes t he [O III A3 727 line observed in the spectrum 
of Mrk 266 NE l |Kim et al.|1995| and also at large rad ii within the 
ionization cones ol many TTONs (e.g. |Tadhunter|2008} . 
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Mazzarella et al. 



HR map(S:1.S-4ksV H:4-7keV) 



Table 6 

Photometry of Knots Embedded in the Arc 0'.'4 west of the NE 
Nucleus 



DEC (deg) 



R.A 


Dec 


m F435W 


H1FS14W 


hh:mm:ss.sss 


dd:mm:ss.ss 


mag 


mag 


(1) 


(2) 


(3) 


(4) 


13:38:17.75 


48:16:40.9 


19.94 ±0.03 


19.02 ± 0.05 


13:38:17.75 


48:16:41.1 


20.09 ±0.03 


18.89 ± 0.04 


13:38:17.75 


48:16:41.3 


20.27 ±0.03 


19.06 ± 0.04 



Note. Columns (1) and (2): J2000 R.A. and Dec, with 

la uncertainty in the absolute astrometry of 0'.'5. Columns (3) 
and (4): Apparent magnitudes and uncertainties measured in the 
HST/ACS F435W (B) and F814 (I) bands, respectively. 

formed additional processing to reveal previously un- 
known structural de tails and physical properties. As 
shown in Figure 18 i, the SW (Seyfert 2) nucleus has 
a much softer X-ray spectrum and lower luminosity than 
the NE (LINER or composite AGN/Starburst) nucleus. 
The luminosity and hard X-ray spectrum of the NE nu- 
cleus provide a clear signature of an obscured AGN. 
The X-ray spectrum obtained of the integrated system 
with XMM (Fig. [6]) shows strong Fe Ka emission at 
6.4 keV. This high-ionizat on feature has been o bserved 

in the 



in the ULIRGs A rp 220 (Iwasawa et al 
LIRG NGC 3690 (jBallo et al ||2004|, 



2005], 
and more recently 
in the integrated (stacked) spectrum of 30 LIRGs in the 
GOA LS sample that are rel atively faint in the hard X-ray 
band ( jlwasawa et~al]|2009| . Modeling the global XMM 
spectrum of Mrk 266 with a thermal (soft, extended) 
component due to star formation and an absorbed power- 
law (hard, compact) component from the AGNs results 
in an absorption column of Nh = 6.5 x 10 22 cm~ 2 , which 
is close to the lower limit s deri ved from the 9.6 /im sili- 
cate absorption feature (j ]3.2.1 ). 

Unlike XMM, Chandra was able to resolve the nuclei 
and has detected two separate sources of (2-7 keV) X-ray 
emission, with the bulk of the combined emission ema- 
nating from the NE component (Fig. |3j). The physical 
origin of the bright Fe Ka line observed in the XMM spec- 
trum of the integrated system is less clear. The Chandra 
spectrum of the total emission from Mrk 266 (Fig. JHJ) 
does not show the Fe Ka line. However, the large aper- 
ture (33" radius) used to produce this spectrum intro- 
duces substantial noise in the background in the higher 
energy bands where there is very little extended emis- 
sion. Detailed inspection of the Chandra data within 
small apertures centered on the nuclei provides a strong 
clue that the Fe Ka line originates primarily from the 
SW nucleu£p*| indicative of a reflection-dominated spec- 
trum of a heavily obscured AGN. Since optical spectra 
reveal the presence of Seyfert 2 activity in the SW nu- 
cleus, the weakness of the hard X-ray emission can be 
understood as a consequence of suppression of X-rays by 
Compton thick absorption. Heavy obscuration towards 
the SW nucleus as inferred from the Spitzer IRS spec- 

29 The 6-7 keV (Fe K) image contains 10 counts from the NE 
nucleus and 4 counts from the SW nucleus, while the neighboring 
3-6 keV band shows a much larger contrast, with 79 counts from 
the NE nucleus and only 1 count from the SW nucleus. In other 
words, the spectrum of the SW nucleus has almost no detected flux 
at energies > 4 keV, but it shows a sudden rise at the Fe K band 
(Fig. pj. This suggests that the hard band spectrum of the SW 
nucleus is dominated by flux from the Fe Ka line. 
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Figure 18. Maps of the hardness of the Chandra X-ray emis- 
sion, (a): The hardness ratio between 1.5 and 7 keV, defined 
as HR hard = (F 4 _7 kcV — F 1 .5_4 kc v)/(F 4 _7 kc v ± F li5 _4 ke v). 
The hard X-ray emission from the SW nucleus is too faint to 
measure this flux ratio, with the exception of the peak pixel 
which has HRi lar( j = —0.4 (blue). (b): The hardness ratio 
between 0.4 and 1.5 keV, defined as HR so f t = (Fo.s_i.5 ke v — 
Fo.4-o.8kev)/( F 0.8-l.5keV± F 0.4-0.8kov); this illustrates the tem- 
perature variation over the diffuse emission. HRhard and HR so f t 
are plotted using colors as indicated in their respective legends. 
The soft band X-ray data are overlaid as contours drawn with lin- 
ear intervals between 1.4% and 12% of the peak brightness (in the 
NE nucleus) to emphasize the extended, diffuse emission. 
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trum and the large column density of obscu ring material 
estimated from the CO (1-0) observations ({3.7) support 
this interpretation of the X-ray data. 

If the 6-7 keV photons in the Chandra spectrum of 
the SW nucleus are due to Fe Ka line emission, as our 
analysis suggests, the flux of 1.2j^g x 10~ 6 ph cm~ 2 s _1 
accounts for the majority of the Fe Ka line flux detected 
in the XMM-Newton spectrum. This means that the 
NE nucleus accounts for most of the escaping hard X- 
ray continuum, but it emits little Fe K emission. The 
equivalent width of Fe K expected from a source ab- 
sorbed by Nh ~ (6 

the NE nucleus (Table [2]), is at most 0.1 keV (e.g.|Awaki 



of extensive shocks in the region ({ 3.2.2 ) that may be 
triggered by outflow from the AGN (T 3.2.4 ) . 



et al.|1991||Krolik et alTp994||Ghiselrini, Haard t k Matt 
1994p , while the equivalent width measured in the XMM- 
Newton spectrum is ~ 0.5 keV. This mismatch can be 
readily explained if the bulk of the Fe Ka line emission 
is from the SW nucleus. The Fe Ka line lum inosity is 
erg s 



3 x 10 40 



1 . Using the relation derived in |Iwasawa| 
et al. (12005b for the Type 2 AGN in the HyLIRG IRAS 



F15307+3252, the intrinsic 2-10 keV luminosity of the 
obscured AGN in Mrk 266 SW is - 1 x 10 42 f" 1 erg s" 1 , 
where f is the (unknown) visible fraction of the reflecting 
matter in the line of sight. The observed F X /F[ Q m] ratio 
for the SW nucleus also indicates a highly absorb ed X - 



ray source with Nh > 10 cm . (See Appendix {5.3 ) 



Finally, it should be emphasized that the Fe Ka line 
detections are marginal. However, we gain confidence 
through independent 2a detections for the global system 
with XMM and for the SW galaxy with Chandra, com- 
bined with line flux values that agree within the uncer- 
tainties in the calibrations and S/N (Table |2|. A detec- 
tion of Fe has also been reported a t the 90% confidenc e 
level with ASCA GIS observations ( |Risaliti et aL]|2000| ). 
Confirmation requires more sensitive observations. 

3.2.6. The LINER in Mrk 266 NE 

The relative strength of the PAH features, silicate ab- 
sorption, and continuum shape indicate Mrk 266 NE 
is typical of an "IR LINER" or "Transition" between 
IR-luminous LINERs and "Type 1" LINERs. Although 
the mid-IR data for Mrk 266 NE are only sufficient to 
compute 3 of the 6 parameters in diagnostic diagrams 



(Fig. 2 of Sturm et al. (20061), the relatively small 
6.2/im/11.2^m PAH flux ratio (1.1) is on the low end 
of the distribution observed in "IR LINERs" (where the 
transition to Type 1 occurs), whereas the continuum flux 
ratios f„(15/mi)/f„(6/mi) = 16 and f l/ (30/im)/f zy (6/xm) = 
54 place Mrk 266 NE on the extreme end of the distribu- 
tion of "IR LINERs" with very red, warm dust continua. 
Mid-IR spectral data alone do not permit a clear dis- 
tinction between the vario us kinds of energy sources tha t 
can power a LINER (e.g., iKewley et al.||2006| [Ho||2008j ). 
However, there are at least tour indicators that Mrk 266 
NE is energized primarily by an AGN, in addition to 
shock excitation: (1) small PAH equivalent widths and 
a deficiency of 6.2 and 7.7 urn PAH emission relative to 
11.3 /im PAH emission ({3.2.1); (2) morphological ev- 
idence for a radiative bow shock within an ionization 



cone aligned with radio plasma in an outflow ({ 3.2.4); 
(3) spect ral pr operties of the luminous, hard X-ray point 
source (^TgJ Sj); and (4) a high H 2 (l - 0) S(l) to Br 7 
flux ratio (Imanishi et al. 2009). Indeed, the strength 
and spatial extent of H 2 emission indicate the presence 



3.3. Between the Colliding Galaxies 
3.3.1. A Search for Counterparts to the Radio Emission 

The bright radio continuum sou rce located between the 
nuclei of Mrk 266 (Figs. [8i and 191 was interpreted as 
a source of enhanced synchrotron emission induced by 
sho cking of the ISM at the i nterface of the merging galax- 
ies ( Mazzarella et al. 1988 ) . Direct imaging of this area 
with H ST ACS and NlCMOS, Spitzer IRAC and MIPS 
( §2.6.1 ), and spectral mapping with Spitzer IRS ([2.7 1 
reveales no obvious source at optical continuum, near- 
IR, or mid-IR wavelengths. This strongly rules out any 
third (dust-obscured) galaxy nucleus or a luminous off- 
nuclear star-forming region as obser ved in other interact- 
ing systems su ch as the Antenna e (Mirabel et al. 19981 
and II Zw 96 ( |Inami et al. 2010). Additonally, there is 
no detectable emission from this region in the near-UV 
or far-UV (Fig. [7] and [8h) , and there is no detection of 
warm molecular gas in the co ntin uum- free H 2 S(3) 9.67 
/im emission-line image (Fig. 10). Although a bridge of 
CO (1-0) emission has be en detected between the g alax- 
ies ( |Imanishi et al. 2009), it is shown below ({3.7 1 that 
this material is located to the north of the central radio 
continuum emission. Next we turn attention to the only 
other waveband where a source of radiation has been 
detected that is spatially correlated with the radio con- 
tinuum emission between the galaxies. 

3.3.2. X-Rays Between the Nuclei 

Substantial soft X-ray emission has been detected be- 
tween the nuclei of Mrk 266 by Chandra. Analysis by 
Brassington et al. (2007) indicates that the X-ray emit- 
ting gas in this region has an intrinsic luminosity of 
9.8(±0.5)xl0 40 erg s _1 and has a temperature of 1.07 
keV, approximately twice the temperature of much more 
diffuse emission surrounding the two galaxies (omitting 
the nuclei and northern source). Brassington et al. sug- 
gested that this X-ray emission is likely associated with 
the collision between the two galaxies, a s proposed previ- 
ously for the radio continuum emission ( Mazzarella et al. 
19881. The present work extends these results using re- 



processing of the Chandra data to enhance the definition 
of structure in the diffuse X-ray emission and to further 
characterize its physical properties. 

To begin with, Figures [3] and [19] illustrate the re- 
markable correspondence between the 0.4-7.0 keV X-ray 
emission and the 20 c m r adio continuum emission. As 
shown in Figures [8] and 19 the higher resolution 6 cm ra- 
dio continuum emission(0'.'3x0'.'4 FWHM), resolves the 
20 cm source into a flattened structure with two sub- 
components. These components are coincident with the 
brightest knot of soft (0.4-4.0 keV) X-ray emission lo- 
cated between the galaxies. The contours of the full-band 
(0.4-7.0 keV) X-ray image in Figure 18 further delineates 
the structure of the diffuse e miss ion between the nuclei 
and to the south-west. Figure [l8p presents the hardness 
of the soft X-ray spectrum (HR so f t ), which indicates tem- 
perature variations within the extended emission. 

Figure [20] provides a detailed view of the enhanced 
soft X-ray emission between the two galaxies. The X-ray 
spectra in the three circled areas show that the SE region 



20 




Figure 19. Mrk 266 Chandra 0.4 - 7 keV (full band) X-ray image 
in pseudocolor, VLA 20 cm radio continuum in white contours, and 
VLA 6 cm radio continuum in black cont ours. T he X-ray data are 
based on image processing described in \ 2.3.11 of this study, and 
the radio data are from |Mazzarella et ah] 



The brightest 20 

cm contours have been omitted so as not to obscure the 6 cm and 
X-ray emission. (See also Fig. The field of view is 45" x 45" 



and the coordinate grid is in J2TJO0 equatorial coordinates 

has a softer spectrum indicating a lower temperature (kT 
0.35 keV) than the other two regions (kT w 0.8 keV). 
These differences are also reflected in the HR so f t map, 
where the SE region has HR so f t = —0.16 and the region 
at the radio continuum peak and 4" NE of that location 
have HR so ft ranging from 0.2 to 0.4. Therefore, the X- 
ray brightness certainly increases at the radio peak, but 
the temperature of the X-ray emitting gas at that loca- 
tion is similar to the region 4" (2.4 kpc) to the NW. A 
caveat in this analysis is that when dust obscuration is 
also varying spatially, HR so f t may not translate directly 
to real variations in temperature. 

Table [7] presents physical parameters estimated from 
spectral analysis of the soft X-ray emission illustrated in 
Figure [20] The most straightforward interpretation is 
that theX-ray emission is the result of gas being shock- 
heated to T ~ 10 7 . This physical process would also 
give rise to the observed enhanced synchrotr on radiation 



1 



(Mazzarella et al. 1988). Davies et al 



S(lj emission from 



(2000) interpreted 
he central region 



the lack ot H2 

of Mrk 266 as evidence for very fast shocks with relative 
velocities of « 1000 km s _1 produced at the interface of 
superwinds emerging from the two galaxies. Such fast 
shocks would generate the observed non-thermal radio 
continuum emission and also heat the gas to > 10 6 K, 
far above the dissociation temperature of H2 . This would 
explain the lack of detected H2 emission lines between 
the galaxies in the near-infrared observations of Davies 
et al. and in the mid-in frare d IRS spectroscopy presented 
in the current study (<j277j). This physical scenario also 
predictes higher temperatures in the shocked region com- 
pared with the diffuse superwind. This is exactly what 
our refined analysis of the Chandra data reveals. 

The estimated cooling time in this region is only 4 Myr. 
This is very short compared to the dynamical time scale 



■ w 




Figure 20. A close-up view of the Chandra soft-band (0.4 - 7 
keV) X-ray emission in the region between the nuclei and to the 
SW. The three circled regions centered on the radio continuum 
peak (70 counts), 4" to the NW (70 counts), and 4" SE (39 counts), 
are used to analyze the hardness of the X-ray emission to deduce 
differences in the gas temperature. 



Table 7 

Spectral Analysis of X-rays Between the Nuclei 



Parameter 


Value 




Units 


(1) 


(2) 




(3) 


0.5-2 keV Flux 


2.7x10" 14 


erg 


— 1 — 9 

3 cm 


Temperature (kT) 


72+ - 11 

u -' -0.09 


keV 




Luminosity a 


1.4xl0 41 


erg 


3" 1 


Volume 


8.1xl0 65 


cm 3 




Electron Density 


l.lxlO" 2 


cm" 


3 


Thermal Energy 


1.8xl0 55 


erg 




Gas Mass 


1.2xl0 7 


Mq 




Cooling Time 


4xl0 6 


yr 





Note. — Column (1): Physical parameter computed from anal- 
ysis of the soft band (0.5-2.0 keV) Chandra emission in the region 
between the colliding disks of Mrk 266. Column (2): Value. Col- 
umn (3): Units. The computations assume an ang ular- scale dis- 
tance of 129 Mpc and the scale of 0.59 kpc/" (see j ]1.3| , and the 
emitting region is approximated by a tube 4" in diameter and 13" 
in length at a position angle of —45°. 

a The bolometric luminosity of this thermal gas, as deduced via 
extrapolation from fitting the data at energies of 0.01-10 keV. 



of the merger (a few hundred million years), and there- 
fore it is likely that we are observing this emission during 
a short-lived phase of the interaction. This would explain 
why so few interacting systems have been found with ra- 
dio continuum and X-ray emission enhanced between the 
galaxies to the same extent as in Mrk 266. However, as 
shown below, there is evidence that gas in this region 
may be reple nish ed from the powerful starburst-driven 
superwinds ((3.4), or from gravitational infall of mate- 
rial toward the center of mass of the system (|3.7[). 



3.4. The Superwind and Northern Loop 

An outstanding characteristic of Mrk 266 is the exten- 
sive filam entary nebula of io nized hydrogen ~30 kpc in 
diameter ( Armus et al.||1990 ) contained within an X-ray 
nebula ~1UU kpc in extent. Optical emission-line diag- 
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nostics and kinematics indicate this is one of the most 
energetic examples of an ou tflowing, starburst-driven su- 
perwind ( Wang et al.|1997 ). Mrk 266 has the third high- 
est Ha luminosity m a sample of 40 local in frared-bright 
galaxies stud i ed by Armus et al.| ( |1990[ ). |Kollatschny| 
|fc Kowatsch| (1998) presented deep observations with 
KUSAT HRl that revealed luminous soft X-ray emission 
spatially correlated with a loop of Ha and [O III] A5007 
line emission ~14" long by 7" wide extending ^10-34" 
(6-20 kpc) to the north of the system (herein called the 
Northern Loop), leading these authors to interpret the 
X-ray emission in this region as hot post-shock gas in a 
high- veloc ity "jet" (most likely emanating from the SW 



nucleus). Kollatschny & Kowatsch (19981 argued also 



that the X-ray emission is too luminous and the opening 
angle of the X-ray/ [O III] emitting region is too small to 
be explai ned by a superwind. A third possibility, sug- 



gested by |Brassington et aL (2007), is that the northern 
structure is a tidal tail containing in situ star formation 
that is responsible for the soft X-ray emission and Ha 
emission. In the subsections below we exploit the new 
multiwavclcngth data to make the case that the North- 
ern Loop is the most luminous component of a superwind 
energized by star formation and one (or both) AGNs. 




Figure 21. Chandra 0.4 - 7 keV (full band) X-ray data contrast- 
ing the high surface-brightness emission (highest resolution) with 
the large-scale diffuse emission. Left (a): Processed for maximum 
spatial resolution and smoothed with a 2-pixel Gaussian kernel 
(same as Fig. |3l bottom left). Right (b): The same data smoothed 
with a 6-pixel Ijaussian kernel to highlight diffuse emission. 



3.4.1. Spatially Correlated Soft X-rays and Ha Filaments 



Our reprocessing of the Chandra data (Fig. 21 1 re- 
veals that the X-ray emission detected in the directnigh- 
resolution data is a high-surface-brightness ridge embed- 
ded in a longer spur that extends an additional ^20" (12 
kpc) toward the NE. A map of the hardness of the soft 
X-ray spectrum (Fig. 18 1 shows the northern region has 
HRsoft ~ —0.2, indicating a relatively cool temperature 
of kT« 0.3 keV. If this is a region of in situ star for- 
mation, one would expect to find a high density of star 
clusters; however, the nearest detected optical star clus- 
ters are more than 5" (3 kpc) SW from the brightest X- 
r ay an d radio continuum emission in the Northern Loop 
(j |3.6[ ). Therefore, the properties of this feature must be 
ascribed to something other than starlight. Figure [22] 
shows that the diffuse soft X-ray emission (0.4 - 1.0 keV) 
corresponds remarkably well with observed Ha filaments 



(a) 





(b) 




(d) 




Figure 22. Soft X-rays compared to other wavelengths: (a): 
Chandra soft (0.4-1 keV) X-ray image re-binned with 2"x 2" pixels 
(grayscale). The contours are the same data smoothed with a 3- 
pixel Gaussian kernel to emphasize the faint, diffuse emission; the 
same contours are also overlaid in the other three panels, (b): 
HST/ ACS B-band (0.44 /im, F435W) image (grayscale) with the 
0.4 - 1 keV X-ray data overlaid as magenta contours a nd the VLA 
20 cm radio continuum image (Mazzarclla et al. 1988) overlaid as 
orange contours, (c): Ha image psn igaki et al. 2U00J with the 
same contours as panel (b). (d): 24 p,m Spitzer Mil's image with 
the same contours as panels (b) and (c). Each field is 1'.2 x 1'.2 
and the scale bar indicates 15". 

and knots that extend out to ~18 kpc (30") to the north 
and to the east of the galaxies. In the northern region, 
the soft X-rays have highest surface brightness near the 
outer edge of the Northern Loop detected in Ha and in 
the B band (including [O II] A 3727 line emission), which 
is also co-spatial with 20 cm radio continuum emission. 
Close spatial correlation between hot X-ray emitting gas 
and warm Ha-emitting gas is a signature of winds in ob- 
j ects ranging from ULIRGs to dwarf starburst galaxies 
( |Grimes et aL|[2?J05] ). 



3.4.2. Structure in the Northern Loop at ^100 pc Scales 

The filamentary morphology of the Northern Loop re- 
vealed by the HST imagery is evidence for an expending 
shell or "bubble" that is fragmenting due to Rayleigh- 
Taylor instabilities in a starburst- or AGN-driven su- 
perwind. Three dimensional hydrodynam ic simulations 
of ga lactic winds predict this phenomena (Cooper et al. 
20Q8|> . Elongated structures « 0.5 - 1" (300-600 pc) in 
length protruding inward from the western, outer edge 
of the Northern Loop are highly suggestive of Rayleigh- 
Taylo r fingers or bow shocks predicted by sim ulations 



(e.g., |Heckman et al.||1990| |Cooper et al.||2008| ). These 
protrusions appear to have a ra dial distribution cen- 
tered on the NE nucleus (Fig. 23), indicating this AGN 



may be a primary ener gizing source for th is co mpo 



nent of the supe rwind. Wang et al. (1997) and Ishi 



gaki et al. (20001 have shown that filaments throughout 



Mrk 266, including the Northern Loop, have line ratios 
(e.g., [O IIIJA5007/H/3) indicating photoionization from 
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Figure 23. HST/ ACS B-band (0.44 (im, F435W) image in 
grayscale with a field of view and intensity mapping chosen to 
emphasize the filamentary structure in the Northern Loop of Mrk 
266. The arrows point out elongated structures 0.5-1" (300-600 
pc) in length on the western side of the loop that appear to be 
aligned radially with the NE nucleus. 

an AGN dominates over ionization from massive stars, 
and shocks that could explain the LINER-like nebular 
spectrum have insufficient energy to account for the lu- 
minosity of the extended nebula. However, the luminous 
starburst in the SW galaxy is also likely to be a major 
contributor to the over all energetics of the superwind. 

The morphological, kinematic and ionization proper- 
ties of the prominent Northe rn Loop can be interpr eted 
with the superwind model of Heckman et al. (19901. In 
this scenario, the star-burst driven gaseous outflow has 
inflated a "bubble" that is highly elongated along the di- 
rection of least resistance (maximum pressure gradient), 
and the closed loop of optical line emission indicates the 
superwind in Mrk 266 is in a relatively young radiative 
cooling phase that occurs before the wind expands fur- 
ther during a "blow-out" phase. This again is consistent 
with the idea that we are observing a relatively short 
lived and dynamic phenomena. 

3.4.3. Dust in the Wind 

The detection of 24um emission co-spatial with the 
Northern Loop (Fig. 22 i) is strong evidence for thermal 
continuum emission from dust grains entrained in the su- 
perwind. By estimating the dust temperature and grain 
emissivity profile, the mass of dust can be de rived follow- 
ing th e single temperature thermal model of Hildcbrand 
(1983), which is given by 



M dust = [(4/3) a p /Q V ][S V D 2 /B(T,v)l 



(2) 



where a is the average grain size (0.1 /an), p is the den- 
sity of grain material, Q u is the grain emissivity charac- 
terized as Q250/™(250/!to/A)^, S v is the flux density of 



the thermal emission at frequency is, D is the distance, 
and B(T, v) is the Planck funct ion. Assum i ng /3 = 1 
and other grain properties from Hildebrand ( 1983 ) re- 

=3~ 



sults in [(4/3) a p/Qn^m] = 9.6xlO~ J g cm" 2 at 24 /im. 
In more convenient units, the above equation becomes 



M, 



dust 



= 0.016 S„(24pm) D 2 [e 



(599.49/T dust ) 



- 1] M(7 



(3) 



where S v {2Apm) is in Jy, D is in Mpc, T is in K, and 
Mdust is in Mq. The 24 /jm flux density measured for 
the Northern Loop (17 mJy, or 2% of the total flux) 
corresponds to a dust mass Md(24um) ss 5 x 1O 7 M0. 



Using the calibration of Alton et al. ( 1999 ) based on the 
known case of M82, the total tar-infrared luminosity ratio 
of L ir (Mrk 266)/L ir (M82) = 18 predicts a dust mass of 
- 2 x 10 7-8 M in the entire wind of Mrk 266. This 
would imply that the 24 pm emission from the Northern 
Loop contains at least 25% (and possibly the majority) 
of the total dust mass within the superwind of Mrk 266. 

However, this result depends on the assumption that 
most of the observed flux is due to dust emission alone 
and that we have made a reasonable determination of 
Tdust ■ One possible source of alternate emission is the [O 
IV] 25.9 pm emission line. We estimate this contribution 
to be less than 10%P^1 Since there is no detection in 
any other mid-IR or sub-mm passband available to help 
constrain Tdust, a characteristic temperature of 37 K is 
assumed based on modeli ng of sub-mm obs ervations of 
the dusty outflow in M82 ( |Alton et alT]|1999| ). 

The overall physical interpretation presented above 
adds strong suppo rt for the superwind interpretation of 
Wang et al. (1997) in which stellar winds and supernovae 
are driving a powerful outflow of gas and dust. The 
Northern Loop is likely the cross section of the most lumi- 
nous "bubble" of expanding ionized gas in the wind. The 
fragmented morphology and radial orientation of protru- 
sions in this region suggest it is strongly influenced by ra- 
diation or a gaseous outflow energized by the NE galaxy, 
and the detection of co-spatial 24 pm emission provides 
strong evidence for dust mixed with the ionized gas in 
the superwind. The new observations are difficult to rec- 
oncile with a jet (ejection from a n AGN) as suggested by 
Kollatschny & Kowatsch (1998) or with a (stella r) tidal 
tail as hypothesized by Hrassin gTon et al.| ( 2007[ ). 



3.5. SED Analysis 
3.5.1. Decomposition of Dust Emission Components 

Figure [24] presents results of a multicomponent model 
decomposition of the IRS spectrum (Fig. 11) combined 
with the broad-band SED (Table [l2j Fig [9} of global 
measurements for the Mrk 266 system using the CAFE 
software package ( Marshall et al. |[200 7p T j A very good 
fit to the SED is provided by a model consisting of PAH 
emission, atomic and molecular line emission, plus cold 

30 This estimate is based on the assumption that the mixture 
of ionized gas and thermal dust in the superwind of Mrk 266 may 
resem ble a scaled up version of a planetary nebula such as NGC 
2346 ( |Su et al.||2004| . The published IRS spectrum of Mrk 266 
shows a bngnt [u IV] 25.9 /im emission line, however it does not 
permit separation of emission from the SW nucleus and the diffuse 
emission in the Northern Loop. 

31 CAFE could not be used to constrain the fits for the indi- 
vidual galaxies because the IRS spectral map could n ot re solve the 
galaxies over the 14-38/im range of the LL module (92.7b. 
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Figure 24. Res ults of multicompone nt fitting using the CAFE 
software package (Marshall ct al. 2007) to model the SED of the 
global emission from Mrk 266. Solid black circles are broad-band 
continuum measurements, and open black squares are IRS spec- 
tral measurements. The solid orange curve represents the sum of 
the various fitted components identified in the legend. The cor- 
responding fit parameters and uncertainties are given in Table [8] 



(26 K), cool (72 K), warm (235 K) and hot (1500 K) 
thermal dust components heated by the combination of 
an unobscured photospheric interstellar radiation field 
(ISRF), an obscured starburst (SB) component with a 
strength constrained by the observed PAH emission, and 
an AGN (accretion disk) component with a strength con- 
strained by the fitted hot dust emission. Uncertainties 
in the derived dust temperatures are ±1 K. 

Table |S] lists the results of the fit. The derived 
extinction-corrected PAH flux ratios appear to be fairly 
"normal," located between the curves for ion ized and 
neutral PAHs in Figure 16 of|Draine fc Li|p00T|). The pa- 



rameter L source is the sum of the derived intrinsic source 
luminosities Ljsrf + Lsb + Lagn, and Ldust is the to- 
tal dust luminosity. The primary results are the relative 
fractions of L source , the fractions of Ldust, and the dust 
masses characterized as arising from cold, co ol, warm and 
hot therm al dust components. As noted in |Marshall et| 
al. ( 2007 ), the four dust components used to model galax- 



les are typically heated to characteristic temperatures of 
Tcoid ~ 35K, T coo i « 80K, T warm « 200K, and T ho t ~ 
1400K, where the actual temperatures (noted above and 
m Table [8| for Mrk 266) are determined from fitting of 
the radiation field energy density of the corresponding 
radiation sources. Although the model suggests that a 
hot dust component accounts for «10% of the total dust 
luminosity, this is uncertain because the fitting of this 
component is weighted heavily by the mid-IR wavelength 
region where emission from hot dust competes with stel- 
lar photospheric emission. Reconciliation of this rela- 
tively low AGN contribution with spectral diagnostics 
that suggest a ~50% AGN contribution in the individ- 



ual galaxies (f 3.2.1 1 may be due primarily to an aperture 
effect. The latter measurements were made in apertures 
that isolate the nuclear regions, and they are thus more 
sensitive to AGN signatures than the larger apertures 
used to construct the SED for the global system which 
encompass more spatially extended starburst emission. 

It is illustrat ive to compare Mrk 2 66 with the galax- 
ies modeled by Marshall et all (120071). Overall, the SED 



Parameter 


ISRF+SB+AGN 


(1) 


(2) 


Twarm(9.7/^m) a 


2.7 i 0.1 


PA IT £ 9/,m /PATT 7 7nm 


n 98 -1- O 0*3 


PATH" 11 ^i/m/PAH 7 7//m 

L Jr\lL X _L . OLlil.il I ± iTll 1 ■ 1 LL III 


n 98 -I- 09 


LlSRF /Lgource 


0.19 


/Lsource 


0.71 


LAGN /Lsource 


0.10 


Lsource (Lq) 


2.4(±0.2) x 10 11 


Lcoid (26 K)/L dust 


0.57 ±0.09 


L CO oi (72 K)/L dust 


0.20 ± 0.03 


Lwarm (235 K)/L dus t 


0.12 ±0.01 


L hot (1500 K)/L dust 


0.11 ± 0.01 


Ldust (L©) 


2.2(±0.2) x 10 11 


M cold (26 K) (M ) 


1.5(±0.4) x 10 8 


M coo i (72 K) (M Q ) 


1.5(±0.6) x 10 5 


M warm (235 K) (M ) 


7.0(±0.2) x 10 2 


M hot (1500 K) (M ) 


6.0 x IO -2 



Note. — Parameters resulting from SED decomposition using 
the sum of ISRF + starburst + A GN radiation sources, as described 
in the text and plotted in Figure [24] 

a These values are likely closer to the intrinsic physical attenuation 
than the much lower value of T warm (9.7/im) — 0.7 ± 0.1 derived from 
spline fitting that estimates only the apparent strength of the silicate 
absorption feature |3.2. 1| 

The L warm /Ld us t values arc corrected for the attenuation by 1 - 
e -r„,m(9.7H/ Tmrm {9 .7^m). 

of Mrk 266 can be modeled as a starburst dominated 
LIRG, with weak AGN indicators from the thermal dust 
or PAH emission, as is also the case in NGC 6240 and 
NGC 2623. While the cold dust component in Mrk 266 
has essentially the same temperature as derived for NGC 
6240 (29 K), the cool (72 K) and warm (235 K) dust com- 
ponents of Mrk 266 are significantly warmer than those 
in NGC 6240 (61 K and 193 K). Together the warm and 
hot components comprise 23% of the total dust luminos- 
ity of Mrk 266. For context, the sum of the (apparent) 
warm and hot dust fractions are 22% in NGC 6240, 77% 
in Mrk 463, and 82% in the QSO PG 0804 + 761. Within 
the uncertainties, Mrk 266 has the same cold dust mass 
as NGC 6240. Compared to the other two dual-AGN 
systems which have had SED decomposition using this 
technique, the Tdust distribution in Mrk 266 is quite sim- 
ilar to NGC 6240, and both of these systems have much 
smaller fractions of warm+hot dust than Mrk 463. The 
Mrk 266 SED fitting results also illustrate that only 700 
Mq of dust at T > 200 K can generate the same luminos- 
ity (4.4 x 10 10 L Q ) as 1.5 x 10 5 M Q of dust at T 72 K. 

3.5.2. Flux Ratios 

The physical nature of AGNs is often revealed via flux 
ratios. W e begin by determ ining the "q" parameter as 
defined by |Yun et al.| ( pOOlj ) , 

W^F^ry), (4) 



qFia- 



[W 



[W 

14/ 



where FIR [W m" 2 ] = 1.26 x 10~ 14 (2.58 • S 60/im [Jy] + 
Sioo/im[Jy]), and the flux densities are those measured 
by IRAS for the global system. In addition, the qpiR 
parameter can be defined as 



q 2 4=log[(S 2 4 M m)/(Sl.4GHz)], 



(5) 

where q24 is the 24 /xm flux density mea sure d by Sp itzer 
MIPS. The data in APPENDIX Tables [ill and fT3l give 
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qFm = 2.0 and q 2 4 = 0.95 for the global system, and 
q 24 = 1.77 (1.14) for Mrk 266 SW (NE). Mrk 266 is 
typical of galaxies in the IRAS 2 Jansky Survey, 98% of 
which fall in the range 1.64 < (qFm) < 3.02 (population 
mean (qFm) = 2.34). A useful comparison for the q 2 4 
values is the sample mean (qFm) = 0.83 ± 0.31 (la dis- 
persion) for over 770 galaxies in t he Spitzer First Loo k 
Survey with redshifts less than 0.4 ( Marleau et al.|2007 ). 
The global system and Mrk 266 NE are within la of this 
sample mean, and although q 2 4 for Mrk 266 SW is large 
compared to most LIRGs, it is within the 3er dispersion 
of the FLS population. Therefore, Mrk 266 shows no 
significant deviations from the infrared-radio correlation 
followed by star-forming galaxies and radio-quiet AGNs. 

The ratio of infrared to ultraviolet luminosities pro- 
vides a diagnostic of the amount of reprocessed versus 
escaping radiation from young, massive stars. Mrk 266 
has an observed ratio of Lj r /Luv = 20, where Luv is 
the sum of the NUV and FUV luminosities measured by 
GALEX, and it has a ratio L uv /(L ir + L uv ) = 0.047 
(5%). With a correction for the finding that Lj r is pow- 
e red ro ughly equally by the AGNs and star formation 
( §3.2.1| ), about 10% of the total star formation in Mrk 
266 is observable directly as escaping UV photons, and 
90% is UV radiation reprocessed into thermal infrared 
emission. Since the star formation rate (SFR) scales di- 
rectly with the starburst component of Luv+Lir, the sys- 
tem has SFRuv/SFRxotai = 10%, which is on the high 
end of the distribution for the GOALS sampl e ranging 
0.2% - 18% with mean 3% ( |Howell et aIj[20lo| . 

Lastly we examine the ratio of nam X-rays to in- 
frared flux. Among 13 local ULIRGs studied by Armus| 
et al. ( |2007 ), those with AGN-dominated mid-IR spec- 
tra have a ratio of hard X-rays to infrared luminosity 
Lnx/Lir > 10~ 3 , whereas starburst-dominated ULIRGs 
have Lnx/Lir < 10~ 4 . The XMM-Newton measurement 
for the global Mrk 266 system gives Lnx/Lir = 3.5x 10~ 4 , 
and the higher resolution Chandra measurements give 
Lnx/Lir = 5.0 x 10~ 5 and L H x/L ir = 1.0 x 10~ 3 for the 
SW and NE nuclei, respectively. The component galax- 
ies of Mrk 266 are therefore at extreme ends of the range 
observed to date in the local (U)LIRG population: The 
value of Lnx/Lir for Mrk 266 NE is 40% larger than NGC 
6240, which itself is similar to ratios fo und in optical/UV - 
selected Seyfert galaxies and quasars ( Ptak et al.||2003 1. 
Mrk 266 SW, on the other hand, has one of the smallest 
Lnx/Lir ratios among l ocal (U)LIRGs, sim ilar to Arp 220 
and IRAS 22491-1808 ( |Armus et al.|2007[). This is consis- 
tent with new evidence presented in §3.2.5| that the X-ray 
source in Mrk 266 SW is Compton-thick. The fact that 
the luminosity ratio of the SW and NE galaxies in hard 
X-rays (L H x(SW)/L H x(NE) = 0.16) is close to the recip- 
rocal of the ratio in the infrared (L ir (SW)/Lir(NE) = 3.3) 
and in cold molecular gas (L C o(SW)/L C o(NE) = 4.9; 
[3.7) is a strong clue that the factor of « 20 disparity 
in the Lnx/Lir ratios for these galaxies is mainly due to 
differences in the column density of material absorbing 
UV/X-ray photons and re-emitting the energy in the far- 
infrared, rather than an order of magnitude difference in 
the intrinsic hard X-ray luminosities of the AGNs. 

This same conclusion about high column density ma- 
terial as the responsible agent for the variance of flux ra- 
tios can be arrived at by simply comparing the two nuclei. 



The data presented in |2.6.2| show that the luminosity ra- 
tio of the two galaxies, TTMrk 266 SW)/L(Mrk 266 NE), 
is ss 0.7 in the far-UV and soft X-ray bands, transitions 
through 1.0 in the near-infrared (1.2 /im), increases to a 
maximum of 5.7 at 8 /jm, and drops to 3.5 at 70 /mi (the 
SED peak) and 3.3 in the total infrared emission (L; r ). 
The L(Mrk 266 SW)/L(Mrk 266 NE) ratio drops further 
to 2.3 at 850 ^m, and to 1.4 at 20 cm. Thus, although 
Mrk 266 SW clearly dominates this merger in the mid-IR 
and far-IR, the two galaxies have similar luminosities in 
the near-IR, UV, soft X-rays, and radio continuum. It is 
only at the highest energies with available data for the 
separate galaxies where the situation is reversed: Mrk 
266 NE is 6.4 times more luminous in the hard (2-7 keV) 
X-r ay ba nd than its companion Mrk 266 SW. As shown 
§3.2.5 this is more likely due to a higher column den- 
~ bsorbing material along our line of sight to the 



m 

sity of a 



SW nucleus (compared to the NE nucleus) than to an 
intrinsic difference in their power sources. 

3.6. Star Clusters 
3.6.1. Cluster Detection and (B - I) Color Map 

The HST/ ACS images of Mrk 266 reveal numer- 
ous very compact, marginally resolved sources scattered 
throughout the system. Their over all properties indi- 
cate they are star clusters (SCs) and associations sim- 
ilar to those observed in other galaxies undergoing ex- 
tensive bursts of star formation. The detecti on method 
foll owed the procedure described in detail by |Modica et| 
aD ( [20121 ) for SCs in the late stage LIRG merger IC 883. 
First, masks outlining the gala xy were applied to the im- 
ages. Then Source Extractor (Bertin & Arnouts 19961 
was used to fit and subtract the underlying galaxy and 
perform point source detection in the B and I band im- 
ages. The extracted objects were processed with IDL 
routines to perform photometry and apply selection cri- 
teria that accepted only reliable sources with a FWHM 
of 1.7-4 pixels, have S/N > 5, and are detect ed in both 
bands; details are given in |Modica et al.| ( [2012[ ). This pro- 
cedure resulted in the detection ot 120 point sources with 
21 mag < B < 27 mag and 19.5 mag < I < 25.7 mag 
(uncorrected for extinction); their coordinates and pho 
tometric pro perties are listed in the Appendix (Table 

Figure 25 shows the direct HST B-band image (a) 
the image after processing with Source Extractor to fit 
and subtract the local background and extract SCs (b). 
The density of SCs is highest in three exte nded regions 
identified with green ellipses in Figure 25 d: (1) within 
and surrounding the disk of the SW galaxy; (2) along 
a band extending ~5-15" (3-9 kpc) north-west from the 
NE galaxy, passing near the base of the Northern Loop; 
and (3) in a group centered ^8" (5 kpc) north-east of the 
NE nucleus. There are 22 SCs within a projected radius 
of 5" (3 kpc) around the nucleus of Mrk 266 SW, and 
only 6 SCs within the same radius around Mrk 266 NE 
(black circles in Fig. 25 d) . The ratio of surface densities 
of detected clusters in the SW and NE galaxies (3.7) is 
nearly th e sam e as the ratio of their infrared luminosities 
(3.5, see [3.5 1, and it is also similar to t he ratio of their 
cold molecular gas masses (4.9, see [3.7 1. This situation 
is consistent with the a pparen t Hubble types and 1.6 
/im bulge-to-disk ratios (j |3.1.1| ) of the two galaxies (i.e., 
SBb disk galaxies contain more gas to fuel star formation 
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Figure 25. Star clusters (SCs) extracted by Source Extractor identified with small circles overlaid on the 0.44 /im (B, F435W) HST 
image before background subtraction (a) and after background subtraction using a median-masking technique (b). Panel (b) also identifies 
notable regions discussed in the text: three areas where most of th e S Cs are concentrated (large green ellipses), each of which contains a 
smaller concentration of very blue SCs (blue ellipses; see also Fig. [2611; regions with radius 5" centered on the nuclei of the two galaxies 



(black circles); and a dearth of SCs bisecting the system at PAs 

than SO/a galaxies). Despite a nearly 4-to-l ratio in the 
number of clusters detected in the central 3 kpc of the 
SW and NE galaxies, across the entire system there is a 
nearly equal number of clusters on both sides of a notable 
gap about 2" (1.2 kpc) wide running between the galaxies 
at PA w -47° (red band in Fig. |25Jd). 

A color map revealing the variation of (B - I) magni- 
tudes throughout Mrk 266 is shown in Figure [261 with an 
overlay of symbols coded to indicate the (B - I) colors of 
the SCs (uncorrected for extinction). Dust lanes present 
in the direct images (Figs. [1} [13] and [14) are more clearly 
delineated in this color map as filamentary regions with 
(B - 1) > 2 mag. In the SW galaxy, the reddest dust lanes 
are on the east and west sides of the nucleus, and aligned 
with the major axis of the underlying starlight detected 
in the 1.6 /jm image (Fig. 3); they are also aligned wit h 
the apparent rotation axis of the molecular gas (see { 3.7 ) . 
In the NE galaxy, the dust lane has an "S" shape that 
crosses the nucleus at PA« —10°. Higher than aver- 
age extinction explains the lack of detected clusters in 
these circumnuclear dust lanes. However, in the region 
between the galaxies the color map shows bluer (B - I) 
values, suggesting the apparent dearth of clusters along 
this gap (red band in Fig. |25|b) may be intrinsic 



(red band). 



Some of the bluest SCs (blue circles in Fig. 26) are 
loc ated in regions identified with blue ellipses in Figure 
[25p , including groupings centered ~7" and ~13" north- 
west of the SW nucleus (inside regions 1 and 2 noted 
above) , and a notable chain «5" in length located north- 
east of the NE nucleus (inside region 3 noted above). 
The detected clusters are a direct probe of star forma- 
tion occurring in regions with relatively low extinction. 
However, the fact that the total luminosity of the un- 
obscured clusters (Mb = —16.25 mag, or Lb = 5 x 10 8 
Lq) is only 0.1% (1.6%) of the total infrared (B-band) 
luminosity indicates that these clusters provide a minor 
contribution to the total energy output of the system. 

The very blue colors of the Northern Loop and fila- 




Figure 26. A grayscale F435W - F814W (B - I) color map of 
Mrk 266, with star cluster positions color-coded by their individual 
(B - I) colors in ranges as follows: [F435W] - [F814W] (B - I) = 
- 0.5 mag: blue circles; 0.5 - 1 mag: green squares; 1-1.5 mag: 
yellow triangles; 1.5-2.0 mag: orange diamonds; (B-I) > 2.0 mag: 
red stars. The scale bar on the right indicates the (B-I) colors of 
the pixels as observed, uncorrected for extinction (reddening). 

ments throughout the system are likely a signature of 
strong nebular [O II] A3727I and [O III] A4364A line 



emission in the B band (see [3.4). The NE nucleus, the 
arc 0'.'4 to the west, and the bi-conic structure in which 
they are embedded are comparably blue, with (B - I) ss 
-0.2 mag. Similarity in color with filaments throughout 
the system that are visible in [O III] emission-line imag- 
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Figure 27. Mn vs. B — I color-m agnitude diagram for SCs detected in Mrk 266. Evolutionary tracks computed using population synthesis 
models jBruzual &: Charlot||2003 I assuming an instantaneous starburst and Solar metallicity for cluster masses of 10 4 (red), 10 5 (green) 
and 10 6 (blue) Mq are shown. Marks on the 10 6 Mq track indicate the age of the starburst. The data are plotted with no extinction 
correction (a), and with a mean extinction correction of Av = 1.2 mag estimated from various measurements as described in the text (b). 
The vector represents 1 mag extinction in the V band. Only SCs with uncertainties in B — I less than 0.2 mag are shown. 

E(B — V) = 0.5(±0.2) mag observed fo r a large sam- 
ple of local, luminous starburst galaxies Ca lzetti et al.| 
(1994). To estimate the total extinction at particular 
wavelengths, we assume Ry = 3.1 as found in the gen- 
eral ISM of our Galaxy. However, this can be considered 
a minimum, because regions containing young star clus- 
ters in Mrk 266 (and other LIRGs) could have a signif- 
icantly larger value of Ry ~ 6 as measured in Galactic 
star-forming regions s uch as the Orio n Nebula and giant 
molecular clouds (e.g., |McCa"lT||2004[ ). 
Application of a mean extinction co rrec tion of Ay = 

I^F 3 ! shifts the 



ing (Ishigaki et al. 2000) lends support to the hypothesis 
presented in [ 3.2.4[ that the circumnuclear region of Mrk 
266 NE is an AGN ionization cone. 

3.6.2. Star Cluster Color-Magnitude Diagram 

Figure [27] is a B — I vs. Mb color-magnitude diagram 
with evolutionary tra cks constructed using pop ulation 
synthesis models from Bruzual & Chariot ( 2003[ ). Most 
of the SCs have observed (apparent) colors of 0.5 mag < 
£B - I) < 2.0 mag and -9 mag < M B < -13 mag (Fig. 
27 1). Although there is an age/mass degeneracy in the 
uncorrected measurements, at first glance most of the 
SCs appear to have ages of a few hundred Myr to 1 Gyr. 
However, given the orientation of the extinction vector 
in Figure [27j which is nearly parallel to the age tracks in 
the population synthesis models, it is clear that the SCs 
are very likely to be significantly younger and/or more 
massive than they appear when uncorrected for extinc- 
tion and reddening due to dust. There are other indica- 
tors for high but variable dust extinctio n: ( a) dust lanes 
evident in the color (B - I) map (Fig. 26); (b) the ob- 



1.2 mag to the observed colors (Fig 
points to an intrinsic mean color of (B 
and to a mean luminosity of Mb 



servation that ~90% of the total luminosit y is g enerated 
by dust emitting in the thermal infrared (E 3.5 ) ; and (c) 
evidence for dust entrained in the su perwmd as far as 
~25"(15 kpc) from the center ( g3.4.3[ >. 

Due to this high variability we can only apply a mean 
extinction estimate using the following indicators: (a) 
The emission- line Balmer decrement (Ha/H/3) measured 
from long-slit optical spectroscopy resulted in E(B-V) 
values of 0.72 mag and 0.22 mag fo r Mrk 266 SW and NE, 



respectively ( Veilleux et al. 1995 



(b) Long-slit spec- 
troscopy along three position angles found E(B — V) = 
0.2-0.5 througho ut the inner -20" (10 kpc) of Mrk 266 
( Wang et al.|1997 ). These measurements indicate a char- 
acteristic extinction of E(B — V) = 0.4 (Ay = 1.2 mag) 
in Mrk 266, which is similar to the sample mean of 

32 The identifications of the NE and SW nucl ei of Mrk 266 were 
mis takenly swappe d in the tables and figures of |Kim et al.| JT995II 
and |Vei!leux et al.| (1995) . 



0.0 mag 
-12 mag (Fig. [27j>). 
Comparison of the mea n extinction-corrected data to the 
model curves (Fig. [27jj) leads to the following results for 
the detected cluster population on average: (a) most of 
the SCs have ages between 1 and 50 Myr; (b) although 
the current data cannot remove the age/mass degeneracy 
for the reddest clusters (B - I > 0.5 mag), about half of 
the objects are very blue (B - 1 < 0.2 mag), very young (1 
- 5 Myr), and are consistent with masses of — 10 5 M ; 
(c) the 11 most luminous SCs have Mb < —13.5 mag 
and are located near the (blue) track for 10 6 Mq clusters 
with ages of only — 1-10 Myr; (d) the two very luminous 
objects with Mb ps -15.8 mag are located 0'.'75 E (the red 
cluster with B - I « 1.1) and Of! 91 W (the blue cluster 
with B - I w -0.6) of the SW nucleus, respectively (see 
Fig. 26 1; and finally (e) if a mean extinction of Ay = 1.8 
mag (Ry w 6) is more applicable than Ay = 1.2 mag in 
regions containing the blue SCs, as argued above, than 
most of the SCs are on average only ~1 Myr old and have 
masses spanning — 10 5 to 10 6 Mq. These results indicate 
the detected SCs arc massive and likely relatively young. 
However, they represent only the "tip of the iceberg" in 
terms of the extensive starburst which is mostly obscured 



33 The Galactic foreground e xtinction in the dir ection of Mrk 
266 is only E(B-V) = 0.013 mag jSchlegel et al.|1998| ; with caveats 
documented in NED, this corresponds to Ay = 0.U4 mag. 
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a = 2.01 ± 0.13 




a = 2.03 ± 0.12 



Absolute Magnitude 

Figure 28. Luminosity functions of the extracted star clusters in the F43 5W (P>) band and F81 4W (I) band before (blue) and after (black) 
correction for foreground contamination and completeness as described bylModica et al.N2012| and Vavilkin et al. (2012, in preparation). 
The solid line represents a linear fit (via \ 2 minimization) to the corrected (black) histogram bins over the ranges —9.75 < Mb < —12.75 
and —10.75 < Mi < —13.25. Alpha indicates the slope of the SC luminosity function in each bandpass over these ranges of Mb and Mi, 
which have not been corrected for the estimated mean extinction (see text). 



by dust in optical and UV observations. 

3.6.3. Star Cluster Luminosity Functions 



Figure 28 is a plot of the luminosity function (LF) de- 
rived for the SCs detected in the B and I band images. 
The index (a) of a power law fit to the LF, after correc- 
tion of the cluster counts for foreground conta mination 
and complet eness using techniques described in Modica 



et al. 



is a 



(|2012h and Vavilkin et al. (2012, in preparation 
TU± 0.1 in both the B and I bands. 
It is revealing to compare the power-law slope of the 
star cluster luminosity function (LF) in Mrk 266 with 
the star cluster LFs observed in other merger systems. 
In the well-studied Antennae Galaxy, the cluster LF has 
a v = -2.13 ± 0.07 tiWhitmore et ah |2010| ). In the early- 
stage mergi ng LIRG system 1 1 Zw 96, an = — 1.1 ± 0.3 



was found ( |Inami et al. 2010). For clusters in two well- 
known ULlRGs that appear to be in intermediate and 
late stages of the merger proc e ss, th e following slopes 

1.6±0.1 
For SCs 



(1998): a B 



were measured by Surace et al 
for Mrk 463, and a B = -1.8 ± ij\2Tor" Mrk 231 
in the late-stage (si ngle nucleus) merger system IC 883, 
Modica et aL] ( 2012| ) derived LF indices (a B = -2.2 ±0.2 
and ai = —2.0 ± 0.2) that are remarkably similar to the 
values found in Mrk 266. In another late-stage LIRG 
merger remnant, NGC 34 (Mrk 938), the star cluste r 
LF has a v = -1.75 ± 0.1 (iSchweizer & Seitzer| 120071). 



(Both IC 883 and NGC 34 are in the GOALS sample, 
with L ir = 5.4 x 10 11 and 2.5 x 10 11 L & , respectively.) 
It is also noteworthy that the LF slopes in Mrk 266 and 
these other (U)LIRGs are similar to the LFs of young 
star clusters in the irregular galaxi es SMC and LMC , 
which have mean a = -2.0 ± 0.2 ( |Gieles et al]l2006[ ). 
The specific frequency of young star clusters is defined 
as T N = Nd x 10(°- 4 ( m b+i5)) ; w here N d is the number 
of clusters above the specified luminosity limit normal- 
ized by the B-band lum inosity of the host galaxy (e.g., 
Larsen &: Richtler] [l999). Correcting for completeness in 
the cluster detection, Mrk 266 has Tn = 0.91, which is 
similar to completeness-corrected Tn values observed in 
the cluster-rich LIRG sample of Vavilkin et al. (2012, 
in preparation). We conclude that a and Tn for star 
clusters in Mrk 266 are similar to other LIRGs. 



Various disruption mechanisms have been described in 
the literature to model the luminosity, mass and age dis- 
tributions of star clus ters in interacting galaxies (e.g. 



Whitmore et al. 
the most 



2007j). Given that Mrk 266 has one of 
ul superwinds known in the local uni- 



verse (j |3.4[ ), this raises the question of whether star clus- 
ters are being disrupted by this outflow/feedback process 
to a larger degree than in systems such as the Antennae 
with much less powerful superwinds and lower star for- 
mation rates. The dearth of detected SCs in the inter- 
action zone of Mrk 266 noted above, despite relatively 
low exti nctio n and molecular gas spanning part of this 
region ( |3.7[ ), and despite the presence of SCs in other 
regions at equally large distances from the galaxy cen- 
ters, may be a sign of cluster disruption or suppression 
of star formation in this narrow band between the two 
galaxies. A possible explanation would be fast shocks 
and extensive turbulence at the interface of the colliding 
galaxies or colliding super win ds expanding from the two 
galaxies, as discussed in S3.3 However, there is no sig- 



nature of cluster disruption on a global scale in Mrk 266, 
as might be revealed through an a or Tn value that is 
significantly different than the SC populations observed 
in (U)LIRGs with less powerful superwinds. 



3.7. Molecular Gas Properties 



3.7.1. 



The Integrated Emission from CO, NCH, and HCO 

Interferometric observations of Mrk 266 in HCN(l-O), 
HCO+(1-0) and CO(1-0) wi th the Nobeyama Mil lime- 



ter Array were presented by Imanishi et al 
this section these data are used to explore t 



(j2009j). In 

le distribu- 



tion of molecular gas with respect to various features in 
the system, and to derive the physical conditions of the 
cold molecular gas. Figure [29] presents contours of the 
integrated emission from these three molecular transi- 
tions overlaid on images ranging from X-rays through ra- 
dio continuum. There are intriguing differences between 
the spatial distributions of these molecular species. Al- 
though the SW galaxy contains about 5 times more flux 
in CO (1-0) than its companion, the NE galaxy is clearly 
detected in all three molecular transitions. The bulk of 
the integrated CO (1-0) emission in the system is aligned 
with the major axis of the SW galaxy, with a centroid 
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(a) 0.4-0.7 keV (Chandra) 




(d) 1.6 fan (2MASS H) 




(g) S(3) H2 9.7 fim (Spitzer IRS) 




(b) U (XMM-OM) 




(e) 1.6 Mm (HST NICMOS) 




(h) 11.3 /im PAH (Spitzer IRS) 




(c) 0.44 ftm (HST ACS) 




(f) 8 fan (Spitzer IRAC) 




(i) 20 cm (VLA) 

A 




in the HST B-band im age, w hich we interpret to be the 
AGN of Mrk 266 NE (P^I). 



Figure 29. The integrated emission of Mrk 266 in CO(l-O) (red contours), HCN(l-O) (cyan contours), and HCO + (1-0) (yellow contours) 
with a spatial resolution of 4'.'3 x 3" 4 (PA= —31.5°) from the Nobeyama Millimeter Array jlmanishi et al,||2009[ l are shown superposed 
on grayscale representations of the following images: (a) Chandra 0.4-7 keV X-ray image, (b) AMM-UM U band, (c) HST ACS 0.44 fim, 
(d) 1.6 jim (2MASS H), (e) 1.6 (im (HST NICMOS), (f) 8.0 jtm (IRAC), (g) H 2 S(3) 9.7 jtm (IRS), (h) 11.3 fim PAH (IRS), and (i) 20 
cm (VLA). The 2MASS 1.6 /tm (H) image is included to show the accuracy of the spatial registration of the 1.6 /im NICMOS image. All 
panels display a field of 20"x 20". The scale bar is 5". 

near the X-ray/infrared nucleus. CO (1-0) is also de- 
tected in a bridge connect ing th e two galaxies, which is 
discussed further below (j ]3.7.2 1. The HCN gas in the 
SW galaxy shows double structure with peaks located 
sal. 2 kpc (2") to the NE and SW of the region where 
the CO (1-0) and HCO + (1-0) are most heavily concen- 
trated at the location of the X-ray /IR/radio nucleus. If 
this morphology is due to limb-brightening in a rotating 
disk, the HCN disk is tilted ~10° with respect to the axis 
of the CO (1-0) and 1.6 /jm emission (Fig 



2!) 



. Most 

of the CO emission in the NE galaxy is offsetTEiy ~1 kpc 
(1'.'5) to the SE of the X-ray/optical nucleus. The CO (1- 
0) in this region is spatially coincident with the peak of 
the HCO + (1-0), whereas the dense gas traced by HCN 
is concentrated very close to the nucleus about 1 kpc to 
the NW. Figure [29^ shows that the centroid of the HCN 
emission is aligned within 0'.'2 of the eastern point source 



Applying the formulation of |Evans et al. ( 2005 1 to com- 
pute L' co , and assuming a CO luminosity to H2 gas mass 
conversion factor a — 1 M© ( K km s -1 pc 2 ) -1 suggested 
by observations of ULIRGs ( |Downes fc Solomon 1998), 
results in estimates of the warm molecular gas mass for 
the two galaxies: 3.4 x 10 9 and 7.0 x 10 s M© for Mrk 266 
SW and NE, respectively. The sum of these components 
corresponds to a total mass of M(H2) = 4.1xl0 9 M©. 
However, the maps also include an extended component, 
and the integrate d CO(l-O) line flux o f the total system 
is 180 Jy km s _1 ( Imanishi et al.|2009 |. The correspond- 
ing total molecular hydrogen gas mass of Mrk 266 from 
the Nobeyama observations is thus M(Ha) = 7.0xl0 9 
Mq. Within the uncertainties, this is consistent with 
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(a) 7850 km/s (112.25 GHz) 



o 




5" 

< > 



(d) 8207 km/s (112.12 GHz) 




(g) 8317 km/s (112.08 GHz) 




(b) 8097 km/s (112.16 GHz) 




(e) 8235 km/s (112.11 GHz) 




(h) 8427 km/s (112.04 GHz) 

m 




(c) 8152 km/s (112.14 GHz) 




(f) 8290 km/s (112.09 GHz) 




(i) 8455 km/s (112.03 GHz) 




I 



Figure 30. The spatial dis tribution of the CO (1 -0) emission as a function of radial velocity (frequency) using a channel map from the 
Nobeyama Millimeter Array jlmanishi et al.||2009| . The CO data are contoured in red and superposed on a grayscale representation of 
the HST 1.6 fim image to illustrate the kinematics of the CO with respect to the old stars in the galaxies. For comparison, the integrated 
emission from HCN (1-0) is plotted (cyan contours) in panel (h), and the integrated emission from HCO + (1-0) is shown (yellow contours) 
in panel (i). The mean radial velocity (and frequency) is shown in each panel. All panels display a 24"x 24" field, and the scale bar is 5". 

206 Jy km s -1 meas ured with single-dish observations 
(Sanders et al. 1986) and a corresponding estimate of 
M(H 2 ) = 2.5 x 10 y M when adjusted to a = 1. If the 
CO to H 2 conversion factor appropriate for Mrk 266 is 
similar to GMC s in our Milky Way (a = 4 ± 1, |Rad- 



Table 9 

Molecular Gas Properties of Mrk 266 



ford et al.|1991 ) , the actual gas masses ma y be ~4 times 
larger than the estimates computed here. ( Sanders et al. 
dl986| > derived M(H 2 ) = 1.4 x 10 10 M assuming a = 4.) 
Assuming a weighted mean dust temperature of 5 K in 
the CO-emitting regions based on SED fitting (j |3.5.1 1 



and using the obs erved CO (1-0) lin e widths, applying 



the formulation of Evans et ah] ( 2005 ) enables estimation 
of the radius of the unresolved (JU-emitting regions, the 
dynamical masses of the galaxies, as well as the column 
densities. These results are presented in Table [9j 

3.7.2. Spatial Distribution and Kinematics of CO (1-0) 



Parameter 


SW 


NE 


Total 


(1) 


(2) 


(3) 


(4) 


CO (1-0) flux [Jy km s" 1 ] a 


88 


18 


180 


CO (1-0) FWHM [km s- 1 ]* 


500 


250 




V CQ [K km s" 1 pc 2 ] 


3.4xl0 9 


7.0xl0 8 


7.0xl0 9 


M(H 2 )] [M ] 


3.4xl0 9 


7.0xl0 8 


7.0xl0 9 


R(CO) [pc] 


208 


134 




^Dynamical 


1.2xl0 10 


1.9xl0 9 




M(H 2 )/M Dynamical 


0.28 


0.37 




Column Density [cm 2 ] 


1.6xl0 24 


7.9xl0 23 





Figure 30 shows the spatial distribution of CO (1-0) in 
velocity channels ranging from 7850 to 8455 km s _1 , in 



Note. — Column (1): Observed or derived parameter. Columns 
(2)-(4): For the SW galaxy, NE galaxy, and the total system, 
respectively — the value of the parameter in Column (1). 
a The CO (1-0) line fluxes and line widths are from [Tmanishi et al.| 
(120091; all other parameters are derived here for the first time. 
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judiciously chosen steps varying from 30 to 250 k m s 1 . 
produced from t he FITS channel map provided by |Iman- 



ishi et al.| ( |2009[) . The CO (1-0) contours are superposed 
on the i/i J'JMCMOS 1.6 /im image for comparison with 
the highest resolution tracer available for the (old) stellar 
mass in the galaxies, and to minimize confusion from dust 
obscuration and young star-forming regions that domi- 
nate the imagery at shorter wavelengths. 

The CO (1-0) emission in Mrk 266 SW shows a strong 
velocity gradient ranging from 7850 km s _1 on the SE 
side to 8235 km s _1 on the NW side of the disk. At pro- 
gressively higher velocities, the centroid of the CO (1-0) 
emission is located at larger radii, extending as far as «5" 
(3 kpc) NW from the SW nucleus. The molecular gas in 
this off- nuclear region has a radial velocity ps 450 km s _1 
greater than the SW nucleus, and this gas coincides spa- 
tially with the most luminous region of off-nuclear star 
formation that is prominent in optical and UV images 
(see Fig. [8]). Although the density of SCs is not signifi- 
cantly higher than in the surrounding areas, t his r egion 
contains many of the most luminous SCs (Fig. 25). 

The low Leo /Lir ratio in Mrk 266 NE is consistent with 
its apparent near-IR Hubble type of SO/a (pec), suggest- 
ing that its pre-merger gas content was relatively low 
and incapable of sustaining a high SFR. Figure [30] shows 
that a large quantity of molecular gas detected over the 
velocity range 8150 — 8320 km s _1 is bridging between 
the galaxies. The highest concentration of CO on the 
NE side of the bridge is centered sal kpc (1'.'5) to the SE 
of the nucleus of Mrk 266 NE. As shown in Figure [29} , 
the CO (1-0) in this off-nuclear region coincides with the 
peak of the HCO+ (1-0) emission, whereas the HCN (1-0) 
emission is concentrated at the NE nucleus. This result 
is consistent with research that has called into question 
the use of HCN as a tracer of dense gas fueling star for- 
mation, suggesting that in some cases HCN is energized 
by an X-ray D issociation Region (XDR) associated with 
an AGN (e.g., |Gracia-Carpio et al.||2008j . 

The distribution of CO (1-0) in Mrk 266 suggests the 
system is currently in a short-lived phase in which ~40% 
of the cold molecular gas is falling toward the center of 
mass of the system, or transferring between the galax- 
ies along a tidal bridge as predicted by merger simula- 



tions (e.g., Cox et al. 2006). Mrk 266 appears to be in 



a similar evolutionary stage as three other LIRGs in the 
GOALS sample that have a large fraction of their total 
CO emissio n located bet ween the nuclei: VV 114 (Yun 
eTal} [1994] ) NGC 6090 (|Bryant fc Scoville||l999"l ), anTJ 
NGC 6240 (Engel et al |[2oT0| ). However, unlike VY 114 



where the 1.4 GHz radio continuum emission has a sim- 
ilar spatial distribution as the CO, the enhanced radio 
continuum emission between the nuclei of Mrk 266 peaks 
«2"(1.2 kpc) south of the CO bridge (Fig. [29}) . Another 
difference is that whereas the bul k of the JiCO + in VV 
114 is located between the nuclei ( Imanishi et al.j[2007 ), 
in Mrk 266 the HCO + is still bound to the galaxies. 

3.8. Wider Implications 

The proceeding sections have shown that Mrk 266 ex- 
hibits the simultaneous occurrence of a variety of phe- 
nomena including dual AGNs, an outflow containing a 
likely radiative bow shock in the narrow-line region of the 
NE nucleus, a starburst with a population of young star 



clusters, a galactic scale "blow-out" expelling dust in a 
superwind, shocked gas between the galaxies, and molec- 
ular gas falling toward the center of mass or transferring 
between the galaxies. In this closing section, Mrk 266 is 
compared with other local major mergers and (U)LIRGs, 
and we attempt to place the system into a tentative 
sequence with a small number of other confirmed dual 
AGNs (with kpc scale nuclear separations) that may rep- 
resent stages in the formation of binary AGNs (accreting 
SMBHs with pc-scale orbital separations). 

3.8.1. Emission Between Interacting Galaxies 

Interacting systems exhibiting substantial radio con- 
tinuum, X-ray, or molecular gas emission between the 
galaxies are rare, but a few have been studied in detail. 
By comparing these systems with Mrk 266 we may be 
able to gain further insight into this phenomenon. 

The Taffy Galaxy (VV 254): This interacting pair 
has «50% of its total 1.5 GHz flux emerging between the 
companions, which has been interpreted as synchrotron 
radiation from cosmic ray electrons trapped in magnetic 
field lines and stret ched like taffy as t he galaxies pass 
through each other (|Condon et al.|1993 ) . Imaging in 450 



and 850 /im continuum emission and in multiple tran- 
sitions of CO indicate the bridge between the galaxies 
consists of interstellar dust and gas extrude d from the 
galax y disks via direct cloud-cloud collisions (Zhu et al. 
2007[) . Mid- infrared observations have detected substan- 
tial PAH emission in the bridge ( jJarrett et al.|[T999| . 

The regions between the components of Tatty and Mrk 
266 both contain large fractions of their total radio con- 
tinuum and CO emission, yet they contain little or no 
evidence for luminous star formation. However, unlike 
Taffy, there is no detected infrared enhancement between 
the galaxies of Mrk 266. In addition, while the radio con- 
tinuum and CO emission are roughly spatially aligned 
and uniformly distributed between the Taffy galaxies, in 
Mrk 266 the central radio emission is much more cen- 
trally concentrated and the bridge of CO (1-0) emission 
is spatially decoupled from the radio continuum and X- 
ray emission. The absence of a discrete mid-IR or far-IR 
source between the galaxies of Mrk 266 suggests shocks 
may have destroyed most of the dust grains in this re- 
gion. The collision in the Taffy system was estimated 
to occur ~ 2 x 10 7 yr prior to its current configuration 
( Condon et al.|1993 ). Its global star formation rate is not 
significantly fiighcr than in non-interacting spiral galax- 



ies (Zhu et al. 



2007), and no evidence for an AGN has 



been published tor either galaxy. 

The Antennae (NGC 4038/4039): This galaxy 
pair is well known for having —30% of its total b olomet- 
ric lumin osity, -40% of its total CO (1-0) flux (|Zhu et 



| al.|2003 ) , and a high co ncentration of young star clusters 
( |VV hit more et ah 1999) located in an off-nuclear region 
between the colliding disks. We have shown that in Mrk 
266 the molecular gas bridging the galaxies is decoupled 
from the central radio/X-ray source, and there is no evi- 
dence for star formation between the galaxies. Compared 
to Mrk 266, the Antennae exhibit: (a) —5 times lower lu- 
minosity in the far-infrared and in CO (1-0); (b) a cooler 
effective dust temperature; (c) a large region of diffuse 
X-rays, but a weaker superwind; (d) a pair of long tidal 
tails rather than diffuse, asymmetric tidal debris; and (e) 
low-luminosity LINER and starburst (H II) nuclei. 
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Early simulations suggested we are viewing the Anten- 
nae galaxies so on after their first encounter (e.g 
ski et al.|1996 



star burst in t 



Dubin 

but a recent model that accounts tor the 



le overlap region constrains the system to 
—40 Myr after the second encounter and only —50 My r 
before the galaxy centers fully merge ( Karl et al.| [2010). 
The close projected separation and small relative veloc- 
ity of the nuclei imply Mrk 266 is in a similar stage of 
the merger process as the Antennae. The lack of well de- 
fined tidal tails in Mrk 266 may be due to any number of 
conditions, i ncludi ng a relatively small disk mass in the 
NE galaxy (j pTO I; a retrograde encounter or an orbital 
geometry (e.g., a smaller periga lactic distance) that p ro- 
duces more diffuse tidal debris ( jDubinski et al.|1996 ); or 
a large ratio of dark matter halo mass to visible (bulge 
plus disk) mass that can cause tails produced early in the 
encounter to fall back onto the galaxies quickly, result- 
ing in late-stage tid al debris resembling loops and shells 
(Mihos et al. 1998). It is noteworthy that the extent of 
the diffuse optical emission in Mrk 266 (« 103 kpc) is 
similar to the end-to-end length of the tails in the An- 
tennae. A detailed simulation of kinematic data for Mrk 
266 is needed to constrain its true dynamical age. 

NGC 6240: Mrk 266 appears to be in a similar evo- 
lutionary stage as the well-known system NGC 6240, 
in which X-ray obs ervations by Chandra have also re- 
vealed dual AGNs (Komossa et al. 2003) and a large 
fract ion of the molecular gas is located between the nu- 
clei dTacconi et al. 1999). Although the total infrared 



luminosity is about halt that of NGC 6240, Mrk 266 has 
a much warmer dust temperature (f60/j m /fioo/im = 0.72) 
compared with NGC 6240 (f60Mm/fioo/xm = 0.26). Like- 
wise, decomposition of the infrared SEDs indicates warm 
(cool) dust components of 225 K (73 K) in M rk 266 com- 
pared to 193 K (61 K) in NGC 6240 (pXTj). The 0.4-10 



keV X-ray luminosity of Mrk 266, 9.5 x 10 ergs s , 
is 57% that of NGC 6240, and within the uncertainties 
of the X-ray flux absorption corrections their intrinsic 
L x /Lj r ratios are similar (5 — 7 x 10~ 4 ). The projected 
separation of the nuclei in NGC 6240 is «1 kpc (1'.'8). 

Stephan's Quintet: This is a well studied compact 
group containing a large intergalactic shock region with 
substantial radio continuum, X-rays, infrared and molec- 
ular gas emission. Unlike Mrk 266, the main shock region 
of Stephan's Quintet exhibits broad (870 km s _1 ) rota- 
tional lines of warm H2 with an order of magnitud e higher 
surface brightn ess than the soft X-ray emission ( | Apple-) 
ton et aL] [2006). This dichotomy may be explained by a 
difference in the shock conditions in the two systems. De- 
tailed modeling suggests the shock velo city in Stephan's 
Q uintet is < 300 km s " 1 ( |Cluver et al]|2010[ ). As noted 
by Davies et al. (2000), if the intergalactic shock front in 
Mrk 266 is produced by the collision of superwinds ex- 
panding outward from the two galaxies at — 300 km s _1 , 
adding this AV — 600 km s -1 to the collisional velocity 
of the two galaxies can produce a net shock velocity ap- 
proaching ~ 1000 km s _1 . This would heat the gas fa r 
above the dissociation temperature of H 2 ( |Lim 2003), 
which is consistent with a temperature of — 10 y K de- 
duced from modeling of the soft X-ray data ( |3.3.2|) . An- 
other difference between these systems is that although 
-40% of the CO (1-0) in Mrk 266 is located in a bridge 
between the galaxies where little or no star formation is 



taking place, in Stephan's Quintet CO (1-0) is detected 
in regions where star formation is occurring as traced by 
24 /jm dust emission and UV continuum emission from 
young, massive stars ( Cluver et al.||2010 ). 

In summary, the radio continuum and X-ray emission 
between the galaxies in Mrk 266 is more similar to the 
large-scale shock detected in Stephan's Quintet than the 
off-nuclear star formation observed in the Antennae. The 
lack of bright H2 line emission and infrared emission from 
dust grains in the central region of Mrk 266 may be due 
to the presence of a very fast shock that effectively de- 
stroys most of the dust grains, sublimates the molecular 
gas, and inhibits star formation. The ridge of — 10 7 
K gas between the galaxies of Mrk 266, which is nearly 
— 10 times hotter than in other regions of the diffuse X- 
ray nebula and has a correspondingly short cooling time, 
is consistent with this scenario. Since the region with 
enhanced soft X-rays and radio continuum emission be- 
tween the galaxi es i s not well aligned with the molecular 
gas bridge (Fig. 29 ) and the physical processes generat- 
ing these features are likely quite different (shocks gener- 
ated via superwinds and gas transfer via gravitation tidal 
forces), we suspect that these two short-lived phases are 
not directly linked. However, their simultaneous appear- 
ance is not a mere coincidence, as both phenomena are 
causally related to the ongoing galaxy merger. 

3.8.2. Will Mrk 266 become a ULIRG? 

It is well established that most local LIRGs and effec- 
tively all ULIRGs involv e interacting and merging gas- 
rich spiral galaxies (e.g., Sanders & Mirabel 1996), and 
the ULIRG phenomenon occurs predominantly in late- 
stage mergers when dust is heated to high temperatures 
by being transported by gravitational torques (with the 
rest of the ISM) close to the central power sour ces (e.g., 
Mazzarella et al.|1991b||Mihos fc Bothun||1998[ ). Clearly 
ULIRGs must pass through a L1KG phase as the star 
formation rate and/or AGN fueling rate increases the 
dust heating to produce a luminosity of L; r > 10 12 L®. 
At least two lines of evidence - one based upon merger 
simulations and one based on molecular gas properties - 
suggest that Mrk 266 will transform into a ULIRG. 

Numerical simulations show the timescale for major 
mergers (from initial approach through final coalescence 
of the nuclei) averages ~ 2 Gyr, but spans 1 — 8 Gyr, de- 
pending on the circularity of the orbit, the v irial radius 
and circular velocity (e.g., Jiang et al. 2008). Compar 
ison of the observed pr operties of Mrk 266 with hydro 



dynamical sim ulations ( Barnes fc Hernquist 1996 



fc Hernquist|l996| |Hop 



tins et al. 



Cox et a 



Mihos 



suggests we are viewing this system about 50% to 70% of 
the way through the merger process. The presence of an 
asymmetric halo of low surface-brightness optical emis- 
sion spanning over 100 kpc, a relatively small projected 
nuclear separation of 6 kpc, and a low relative veloc 
ity of th e nuclei of AV = 135 km s _1 (Khachikian et 



aL]|2006 ) indicate the two galaxies are likely undergoing 



their second or third encounter, with each passage damp- 
ing the relative velocities of the nuclei, strong dynami- 
cal friction causing the orbit to decay, and tidal t orques 
drivin g strong inflows of gas. As pointed out by Mihos 
(1999), the timescale for violent relaxation is > 1 Gyr 
in the outer regions (low V c ), but only — 100 My in the 
inner regions (high V c ). Thus, depending on whether 
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the nuclei are on the verge of free-fall (which can oc- 
cur as rapidly as « \ of a rotation period) or whether 
they will orbit one more time before fully merging, the 
estimated time to coalescence is about 50 to 250 Myr. 
(We are referring here to the merger of the stell ar sys - 



tems, not the SMBHs. The latter is discussed in \ 3.8.3 
Based on Monte Carlo simulations of merger models to 
match the observed distribution of nuclear separations in 



ULIRGs ([Murphy el al.| | 1996[ ), Mihos (1999| concluded 
that most ULIRGs are In the final 20% of the merger 
process. Thus, Mrk 266 appears to be in an interesting 
evolutionary stage when the nuclei are about to enter the 
final phase of coalescence characteristic of ULIRGs. 

The total cold molecular gas mass of Mrk 266 is sim- 
ilar to that o f ULIRGs observed locally and at high 
redshifts (e.g., Solomon & Vanden Bout 2005). Given 
that «40% of the total CO (1-0) emission is m a bridge 
between the galaxies, there is a large reservoir of cold 
molecular gas available to form more stars or to fuel the 
AGNs as the galaxies and nuclei fully coalesce. With 
L co = 6.7 x 10 9 L and L ir /L' co = 50.4, Mrk 266 
falls directly on the well known correlation between 



log(L j r /L^ ri ) and log(L; r ) ( |Solomon fc Vanden Bout 



2005 1 . W ithin the evolutionary scenario of Sanders et al 



qi988a|pp , we expect Mrk 266 will evolve into a ULIRG 



initially with a relatively "cool" global far-infrared dust 
temperature similar to Arp 220, follo wed by a "warm' 



ULIR G phase exemplified by Mrk 463 ( Mazzarella et al. 
1991a Bianchi et al.|[2008 ). Given the new evidence tor 



d ust en tramment in the extensive superwind of Mrk 266 
( ^3.4.3 1 , this system illustrates that the "blow-out phase" 
when large quantities of dust are expelled to reveal pre- 
viously obscured optical/UV AGNs can begin during an 
intermediate stage of the merger before the nuclear coa- 
lescence phase when ULIRGs are most likely to occur. 

ULIRGs have a much lower space de nsity than LIRGs 
locall y and at intermediate redshifts ( |Le Floc'h et al." 
2005|. In addition, the mean luminosity of ULIRGs 



in the local universe is lower than ULIRGs and sub- 



millimeter galaxies at high redshifts (e.g., Solomon & 
Vanden Bout|20 05 ) . The resolution of these observations 



may well lie in the pre-merger gas content of the partic- 
ipant galaxies. The higher that gas content, the larger 
the burst of star formation and subsequent dust heat- 
ing. In the case of Mrk 266, the NE component contains 
w 5 times less molecular gas than its companion; this is 
consistent with its apparent morphological type (SO/a) 
and correspondingly high bulge-to-disk ratio found at 1.6 
imp.. Using the well-e stablished correlation of Leo with 
L; r (Young et al. 1986 ), a merging pair involving two late- 
type spirals with molecular gas masses similar to Mrk 
266 SW (SBb) would have Leo a f ew times higher than 
is observed in the actual Mrk 266 (SBb + S0/a) pair, re- 
sulting in a t otal Lj r at the high end of the lo cal ULIRG 



distribution (Solomon & Vanden Bout 2005) 



In the hierarchical model, a galax y's Hub 



)ble type is a 

byproduct of its merger history (e.g., Schweizer 20001. In 



this scenario, although merging played an important role 
in galaxy evolution over all cosmic time, high-z galax- 
ies and proto-galaxies were very different than those ob- 
served at intermediate and low redhifts, and prior to 
z ~ 1 the Hubble sequence was not yet in place. However, 
from z ~ 1 to the present, because many late- type spiral 



galaxies merged and transformed into early-type spirals, 
SO, or elliptical galaxies, the fractio n of late-type spi - 
rals has decreased substantially (e.g., Oesch et al.|2010 ). 
Locally at z ~ 0, this makes pairings of gas-rich, late- 
type galaxies less likely than pairings of one late-type 
with one early-type galaxy (or of two early-type galax- 
ies). Therefore, at low redshifts a smaller fraction of 
major mergers will have sufficient gas to fuel star forma- 
tion or AGN accretion at the level required to reach the 
ULIRG threshold during early to intermediate phases of 
the merger. As argued above, however, this would not 
prevent a ULIRG from developing later if one of the com- 
panions has > 7 x 10 9 Mq of molecular gas available to 
be driven toward the center of mass of the system to fuel 
two massive AGNs as they inspiral toward final coales- 
cence. Mrk 266 is apparently in this physical situation. 

3.8.3. Evolution from Dual to Binary AGNs 

In this final section we turn attention to the signifi- 
cance of Mrk 266 in context with other confirmed dual 
AGN systems in an attempt to understand the low de- 
tection rate of binary AGNs with sub-parsec scale or- 
bits found in recent surveys. We also propose a putative 
evolutionary sequence involving a quantitative parame- 
ter available for a small number of major merger systems 
with confirmed dual AGNs and candidate binary AGNs. 

Current t heory predicts the fo llowing evolutionary se- 
quence (e.g., Colpi &: Dotti|2009 and references therein): 
(a) pairs of massive galaxies separated by many times 
the diameter of their luminous disks and bulges will 
merge in ~1 Gyr (a few dynamical time-scales) due to 
dynamical friction between their massive dark matter 
halos; (b) the two pre-existing SMBHs in their cores 
(with M S MBH/Mbuigo ~ 0.1%) will form dual AGNs as 
tidal torques drive the ISM inward to fuel accretion at a 
higher rate than existed prior to the encounter; during 
this stage the nuclei are merging from a separation of 
> 1 kpc to ~10 pc; (c) in many cases conditions may 
be conducive to forcing the black holes pairs to inspiral 
over a period of ~ 100 Myr to a separation of < 10 pc, 
at which point they enter Keplerian orbits and may be 
detectable as binary AGNs; (e) the SMBHs will even- 
tually either completely coalesce or else recoil, resulting 
in the most powerful sources of gravitational waves pre- 
dicted to exist. Considering the relatively small number 
of dual AGNs confirmed to date, it is tempting to con- 
clude that it is rare for both SMBHs to have accretion 
rates high enough to power two luminous AGNs simul- 
taneously. However, given that the column density of 
gas and dust in (U) LIRGs is high enough to make the 
central ~ 100 pc optically thick in the mid -infrared, and 
ev en in t he far-infrared in many cases (e.g.,|Thompson et 



al. 



2005), many more dual AGNs will likely be detected 



h observations in other spectral regions. Indeed, X- 
ray surveys sug gest ~60% of optically classifie d LINERs 



are true AGNs (Gonzalez-Martin et al. 



2006), and high 



resolution radio continuum imaging finds a high fraction 
(74%) of AGNs amon g optically classifie d LINERs in the 
1 Jy ULIRG sample ( Nagar et aL"fl2003[ ) . 

Along with Mrk 266 and 1NGC 6240, the following lo- 
cal (U)LIRGs have also been confirmed to contain dual 
AGNs based on detection of two compact, luminous 
sources of h ard X-rays: Mrk 1 71 (Arp 299 = NGC 3690 



+ IC 694) ( |Ballo et al.||2004| ) and Mrk 463 ( |Bianchi et 
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al. 2008). Mrk 171 is similar to NGC 6240 in terms 



of infrared luminosity (Lj r = 8.5 x 10 11 L©), and pro- 



However, with 
(uncorrected for 



jected nuclear separation (0.9 kpc, 3'.'c 
L x (2 - 10 keV) = 2.7 x 10 41 ergs s" 
absorption), Mrk 171 has L x /L ir = 8.4 x 10~ 5 , which is 
an order of magnitude smaller than Mrk 266 and NGC 
6240. With L x (2 - 10 keV) = 3.8 x 10 42 ergs s" 1 , Mrk 
463 has a hard X-ray luminosity 9 times greater than 
Mrk 266 and 16 times greater than Mrk 171. Mrk 463 
has L x /Lj r = 1.9 x 10 -3 , which is —2-5 times larger than 
NGC 6240 and Mrk 266, and -23 times larger than Mrk 
171. A remarkable dual AGN is 3C 75, in which both nu- 



clei are emi tting powerful double-sided radio jets (Owen 
this system consists of a pair of elliptical 



et al. 



1985) 



galaxies with projected separation 6.4 kpc (14"). 

There is growing evidence that the current small 
number of confirmed dual AGNs will increase, primar- 
ily with observations of hard X-rays. For example, 
in a local sample of hard X-ray selected AGNs, 10% 
have been found to reside in dual sy stems with sepa- 
rations < 100 kpc ( |Koss et al^|2012| P1 Dual AGNs 
and binary AGN candidates are also Deing found in 
surveys probing to higher redshifts. The z = 0.36 
system COSMOS J100043. 15+020637.2 has been con- 
firmed to contain tw o luminous AGNs wit h a separa- 
tion of 2.5 kpc (0'.'5 ) dComerford et al.|[2009| ). In CXOC 
J100043. 1+020637, [Ovano et al.| ( |2010f have interpreted 
the presence of optical morphology indicating a major 
merger, a large velocity offset, and an X-ray iron line 
with an inverted P-Cygni profile as a recoiling SMBH or 
a slingshot effect in a triple system. The published SED 
indicates L x /L ir 0.018, which is orders of magnitude 
larger than in the local (U)LIRGs. Four more dual AGNs 
with kpc-scale separations and 0.13 < z < 0.21 have al so 
been discovered recently in the SDSS ( Liu et aLp OlO). 

What appears to be the first bona fide binary AGIN was 
discovered in 4C +37.11, which contains two compact ra- 
dio sources with a projected separation of only 7.5 mas 
(7 pc). Based on the radio source compact ness, variabil- 
ity, motio n in the jets, and spectral shapes, [Rodriguez et| 
al. ( 2006 ) concluded that both radio cores contain AGNs 
orbiting with a rotational period of 1.5 x 10 yr within 
a single elliptical galaxy host. It is notable that the sum 
of the SMBH masses in 4C +37.11, - 1.5 x 10 8 M , is 
similar to Mrk 266; this implies that LIRGs with bulge 
masses similar to the companions in Mrk 266 may end up 
as remnants similar to 40 +37.11. The soft X-ray lumi- 
nosity of 4C +37.11 is « 170 times larger than Mrk 266. 
(Hard X-ray observations of 4C +37.11 are not available.) 

Very few additional binary AGNs have been confirmed, 
although there are a few candidates. Observations of 
double-peaked emission lines in the radio-quiet QSO 
SDSS J153636. 22+044127.0 (z = 0.388) have been inter- 
preted as a candidate binary SMBH with masses 8 x 10 8 
and 2 x 10 7 inferred to be sep a rated by only 0.02 mas (0.1 
pc) ( |Boroson k LauCT|[2009j) . IWrobel k Laor| ( [2009] ) re- 
ported detection ot two 8.5 GHz radio continuum sources 
separated by / .'97 (5.1 kpc) within the aperture of the 

34 However, it should be noted that many pairs with such wide 
separations may not be gravitationally bound, and a hard X-ray 
selected sample is naturally biased toward a much higher AGN frac- 
tion than an infrared-selected sample. The number of confirmed 
dual AGNs is currently less than 3% of the GOALS sample. 



original spectroscopic observations, suggesting this is an- 
other kpc-scale dual AGN system. However, further ob- 
servations show the second radio source is an elliptical 
galaxy that i s not the origin of the r ed- or blue-shifted 
Balmer lines ( Lauer k Boroson 2009 ) . Another possibil- 
ity is that SDS« J 153636. 22+U44I27.0 is similar to other 
"double-peak" AGNs such as 3C 390.3 and Arp 102B, 
which were also promoted as candidate binary AGNs but 
later found to be more consistent with a rapidly rotating 
accre tion disk within a single broad-line region (|Gaskell| 
2010| . Additional observations are needed to distinguish 



between these alternate hypotheses. 

Table [10] lists basic parameters for the small number 
of confirmed dual AGNs and putative binary AGNs 
discussed above. Due to the high degree of uncer- 
tainty in age dating mergers with factors such as 
nuclear separation, the morphology of tidal debris, 
and varying techniques applied in the simulations of 
specific systems, we propose a merger sequence based 
on the quantitative parameter L x /Lj r , where L x is 
the 2-10 keV hard X-ray from XMM where available, 
otherwise the 2-7 keV band from Chandra. Accord- 
ing to the evolutionary paradigm whereby ULIRGs 
transform into U V/optical QSOs during t he merger of 



gas-rich galaxies (Sanders et al 
to which the luminous optical/ U V /A-ray 
from the accretion disk is obscured should 



1988a|bJ, the degree 
emission 
diminish 



as superwinds driven by star formation and/or AGNs 
diffuse or expel the dust from the system. To the 
extent that increasing log(L x /L; r ), given in parenthe- 
sis below, is a reliable gauge of the degree to which 
the embedded AGN-powered UV/X-ray sources have 
emerged from their obscuring cocoons, a presumptive se- 
quence is: Mrk 171 (-4.1) -» Mrk 266 (-3.5) -> 
NGC 6240 (-3.1) Mrk 463 (-2.8) -> 

COSMOS J100043.15 + 020637.2 (-1.7)FJ 

Despite the growing number of Kpc-scale dual 
AGNs, the number of sub-parsec binary AGNs in the 
class of 4C +37.11 remains very small. Currently 
SDSS J153636. 22+044127.0 is the only candid ate out of 
« 17, 500 AG Ns with z < 0.70 in the SDSS ( |Lauer k 
Boroson 2009 ) . Similarly, a recent study of VLB A data 



tor over 3100 radio- luminous AGNs fou nd only one bi 
nary, the previously known 4C +37.11 ( Burke-Spolaor 
2011|). Since major mergers provide a natural process 



to form binary SMBHs, t he appa r ent d earth of binary 
AGNs needs explanation. |Gaskell| ( |2010[ ) has noted two 
possibilities: either the lifetime of binary SMBHs is much 
shorter than predicted (i.e., they inspiral very rapidly 
rather than remaining in a stalled orbit), or fueling of the 
accretion disk(s) is somehow suspended. The former ex- 
planation seems unlikely because the very presence of the 
accretion disks forms a significant buffer against rapid 
inspiral and effectively "softens" the SMBH merger pro- 
cess. A possible physical mechanism for the latter sugges- 
tion may be that AGN jets have the potential to impart 
enough energy into the inner parsecs of a galaxy to im - 
pede gas inflow to the accretion disk ( [Nagar et al.|2 005 ) . 
We propose another physically viable alternative: If dual 
AGNs are confirmed in many more (U) LIRGs in inter- 

35 SDSS J153636.22+044127.0 is omitted from the proposed se- 
quence because its nature as a true binary AGN is controversial, 
and the required IR and X-ray measurements are not yet available. 
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Table 10 

Dual to Binary AGN Evolution Based on L x (Hard)/L; r 



Object 


z 


Dl 


Sep 


Sep 


&V L ' 
(km/s) 


log( -^4 '- ) 


Notes 


Refs 






(Mpc) 


(") 


(kpc) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Mrk 171 (NGC 3690 + IC 694) 


0.01041 


51.1 


3.8 


0.9 


11.88 


-4.1 




1 


Mrk 266 (NGC 5256) 


0.02786 


129 


10 


6.0 


11.53 


-3.5 




2 


NGC 6240 


0.02448 


116 


1.8 


0.9 


11.93 


-3.1 




3 


Mrk 463 


0.05036 


233 


3.8 


3.8 


11.77 


-2.8 


Young radio jet? 


1 


3C 75 


0.02315 


97.9 


14 


6,1 






Twin radio jets 


5 


COSMOS J100043. 15+020637.2 


0.36060 


1944 


0.5 


2.5 


11.59 


-1.7 


Recoiling SMBH? 


6 


4C +37.11 


0.05500 


242 


0.0073 


0.008 






Binary AGN 


7 


SDSS J153636.22+044127.0 


0.38930 


2133 


1.0? 


10~ 4 or 5.1 






Binary or dual AGN? 


8,9 



Note. — Basic parameters for confirmed dual AGNs and putative binary AGNs. Column (1): Object name. Column (2): Heliocentric 
redshift from NED (see citations therein). Column (3): Luminosity distance computed by NED. Column (4): Projected nuclear separation 
in arcseconds. Column (5): Projected separation in kpc. Column (6): Log of the infrared luminosity in Solar units. Column (7): Log of the 
ratio of observed hard X-ray luminosity (2-10 keV from XMM, otherwise 2-7 keV from Chandra, uncorrected for absorption) to infrared 
luminosity. Column (8): Notes. Column (9): References for the redshi ft in column (2), t he hard X-ray flux use d to compute the ratio in 
co lumn (7), or the Not es in column (8); the co des a re as follows: 1 =|Ballo et al.| (j2004|t; 2 = this article; 3 =|Komossa et al.|([2003} ; 4 
= |Bianchi et al.| l|2008 [); 5 = |Owen et ah] | |1985| ; 6 = |Civano et al.| j201Tj{ ; 7 = |Kodriguez et al~|j2006| ; 8 = |Boroson fc Lauer| ^M^ Vl = 
Wrobel & Laorl 12009b . Measurements taken from the literature have been adjusted to H = 70 km s _1 Mpc 1 , il M = 0.27, Hy = 0.73. For 
COSMOS J10U043.15+020637.2, L ir was estimated from the Spitzer IRAC 8 ^m and MIPS 24 /im measurements in |Civano et al.| | |2010"| l 
utilizing equation^ of the current study. 



mediate stages of the merger process with kpc scale sep- 
arations (i.e., if most of the LINERs in the GOALS and 
IRAS 1 Jy ULIRG samples are confirmed to be ionized 
by a pow er-law source as sug gested by the radio obser- 
vations of Nag ar et al.| ( |2003| ) , this raises the possibility 
that there will be an insufficient fuel supply remaining to 
power luminous binary AGNs late in the merger process. 

The lifetime of nuclear starbursts and AGNs is in- 
versely proportional to the efficiency of gas depletion 
and ranges betw een ~ 10 7 for ULIRGs and ~ 10 8 yr 
for LIRGs (e.g., [Solomon fc Vanden Bout|[2005l ). The 
lifetime of powerful radio-lou d AGIMs has also been es- 



2008J). This is 
re final coales- 



timated to be only ~ 10 7 yr (Bird et al 
similar to the dynamical time scale of t 
cence of the nuclei with high circular velocity at the core 
of a merger, but only 1-10% of the ~ 10 9 Gyr time scale 
for dynamical relaxation of the stars in outer orbits of 
the remnant galaxy (e.g., |Mihos 1999). Therefore, it is 
possible that starbursts and dual AGIN s may consume or 
expel the bulk of available gas relatively early in many 
major mergers, leaving insufficient fuel to power QSOs 
that can be detected spectroscopically through the stage 
when the SMBHs inspiral to form a true binary and then 
coalesce or recoil. This hypothesis counters simulations 
that find AGN a ccretion rates peak du ring the SMBH 
coalescence (e.g., Johansson et al.p009 ). 



4. SUMMARY AND CONCLUSIONS 

We have presented the results of a multi-wavelength 
campaign on the luminous infrared galaxy (LIRG) sys- 
tem Markarian 266 that has yielded the following conclu- 
sions regarding the morphology, nuclei, dust, molecular 
gas, outflows, star clusters, and radiation between the 
nuclei, as well as implications for a better understanding 
of the origin of ULIRGs and binary AGNs. 

1. HST images show that the "double- nucleus galaxy" 
observed in low-resolution images is a strongly interact- 
ing pair; we classify the SW and NE galaxies as SBb 
(pec) and SO/a (pec), respectively. Deeper imaging re- 
veals asymmetric tidal debris spanning wlOO kpc. At H 
band, the companions are more luminous than L* galax- 



ies, and their estimated stellar masses are 6.1 x 10 10 Mq 
(SW) and 4.4 x 10 10 M Q (NE). The stellar bulge lu- 
minosities imply remarkably similar nuclear black hole 
masses of 2 x 10 s M . 

2. The NE (LINER) nucleus is emitting 6.4 times more 
hard X-ray (2-7 keV) flux than the SW (Seyfert 2) nu- 
cleus, yet the luminosity of the NE galaxy is only w20% 
of the SW galaxy at 24 /im and in CO (1-0) molecular 
gas emission. The Chandra spectra indicate that the SW 
nucleus is likely the primary source of a bright Fe Ka line 
detected in the XMM spectrum of the total system, con- 
sistent with a reflection-dominated X-ray spectrum of a 
heavily obscured AGN. The F(X-ray)/F([0 III]) flux ra- 
tios also indicate the SW nucleus is much more heavily 
obscured than the NE nucleus. 

3. Soft (0.4-2 keV) X-ray emission extends 15 kpc to 
the north of the system, between the nuclei, and 3 kpc 
west of the SW nucleus. The smoothed X-ray structure 
corresponds remarkably well with filaments in Ha and 
[O III] images. A ridge of X-ray emission between the 
galaxies has an X-ray spectrum consistent with T ~ 10 7 
K shock-heated gas, strengthening the idea that the cor- 
responding radio emission is shock-induced synchrotron 
radiation. The lack of detected 9.7 /im H2 S(3) line emis- 
sion supports the suggestion that a very fast shock at the 
interface of superwinds expanding from the two galaxies 
might generate the observed non-thermal radio contin- 
uum emission and heat the gas above the dissociation 
temperature of H2. The derived cooling time of the X- 
ray emitting gas between the nuclei is only «4 Myr. 

4. The HST optical images reveal a circumnuclear arc 
containing three knots 240 pc west of the NE nucleus. 
The arc is embedded within B-band emission with a fila- 
mentary, bi-conic morphology extending over a radius of 
1.2 kpc. This region appears to be dominated by nebu- 
lar emission from the narrow-line region of a double-sided 
AGN ionization cone. Radio continuum emission at 1.6 
GHz peaks at the optical/near-IR nucleus and has com- 
ponents in alignment with the nucleus and the circumnu- 
clear B-band knots (PA=56°). A plausible explanation 
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is that the optical knots are clouds entrained in a shock 
front produced by an AGN-powered collimated plasma 
outflow (jet). We liken this region to the radiative bow 
shock 230 pc south of the LINER in M51. 

5. Mid-IR spectral diagnostics suggest the bolomet- 
ric luminosity is powered by roughly an equal mixture 
of radiation from AGNs and star formation, with sub- 
stantial scatter among the various methods. Newly con- 
structed SEDs give infrared luminosities of L; r = (2.3, 
0.7, 3.4) xlO 11 L Q for the SW, NE, and total system 
respectively. Decomposition of the SED of the total sys- 
tem indicates that cold (26 K), cool (72 K), warm (235 
K) and hot (^1500 K) thermal dust components con- 
tribute approximately 57%, 20%, 12%, and 11% of the 
bolometric luminosity, respectively. The total cold dust 
mass estimate is 1.5(±0.4) x 10 8 M Q . 

6. HST imaging resolves the Northern Loop extending 
6-20 kpc from the galaxies into a fragmented morphol- 
ogy suggestive of Rayleigh- Taylor instabilities. A few 
structures 300-600 pc in length are radially aligned with 
the NE nucleus, suggestive of bow shocks predicted by 
hydrodynamic simulations of galactic superwinds. De- 
tection of 24 /xm emission in the Northern Loop implies 
that at least ~ 2xlO 7 M0 of dust is being swept out of 
the system by the superwind during a "blow-out" phase 
that is well underway prior to the final galaxy merger. 

7. Approximately 120 star clusters detected with HST 
are concentrated in the SW galaxy and at the base of 
the Northern Loop; most have estimated ages less than 
50 Myr and masses of <~ 10 5 M . The ratios of cluster 
surface densities, Leo and Lir within 3 kpc of each nu- 
cleus are similar (i.e., L(SW)/L(NW) w 4- 5). The star 
cluster luminosity function is similar to what has been 
found in other LIRGs, and the unobscured clusters con- 
tribute little to powering the total infrared luminosity. 

8. We conclude that Mrk 266 NE is an AGN based 
on the presence of: (1) an obscured, hard X-ray point 
source; (2) a radiative bow shock aligned with ioniza- 
tion cones and apparent radio plasma outflow; (3) PAH 
emission with small equivalent widths and a deficiency 
of 6.2 and 7.7 ^m flux relative to 11.3 (im flux; and (4) 
a H 2 (l - 0) S(l) to Br7 flux ratio of 3.5. The spatially 
extended H2 emission demonstrates that LINERs can ex- 
hibit strong signatures of an AGN and shock excitation. 

9. The bulk of the CO (1-0) emission in the system is 
aligned with the major axis of the SW galaxy. HCN (1-0) 
is aligned with H 2 S(3) 9.7 fim emission in the nucleus 
of Mrk 266 NE. However, the CO (1-0) and HCO+ (1-0) 
emission peaks are located ~1 kpc SE of the NE nucleus, 
suggesting that either most of the cold molecular gas 
has already been stripped from the NE galaxy, or gas is 
being transferred from the SW galaxy to the NE galaxy. 
Approximately 40% of the total CO (1-0) emission is 
bridging the galaxies, likely in-falling toward the center of 
mass or transferring between the galaxies. In this regard, 
Mrk 266 appears to be in an evolutionary stage similar 
to VV 114, NGC 6090, and NGC 6240. 

10. Two lines of evidence suggest Mrk 266 may evolve 
into a ULIRG: First, in the context of merger simula- 
tions, the galaxies lie within an asymmetric, low surface- 
brightness halo of tidal debris spanning «100 kpc, with 
a projected separation of 6 kpc and a relative velocity of 
only 135 km s _1 , indicating they are likely undergoing 



their second or third encounter with only 50 to 250 Myr 
remaining until they merge via tidal dissipation. There- 
fore, Mrk 266 appears to be in a short-lived stage when 
the nuclei are about to enter the final phase of coales- 
cence characteristic of ULIRGs. Second, the total cold 
molecular gas mass of w 7 x 10 9 M is similar to local 
and high-redshift ULIRGs. Since «40% of the total CO 
(1-0) emission is located between the galaxies, this reser- 
voir is available to form more stars and to fuel the AGNs 
as the stellar systems and nuclei inexorably coalesce. 

11. We propose that Mrk 266 belongs to an evolution- 
ary sequence in which dual AGNs with kpc separations 
represent precursors to putative binary AGNs with pc 
scale orbital radii. In this scenario, where Mrk 266 is in 
an intermediate phase between Mrk 171 and NGC 6240, 
the global L x /L ir ratio increases by over four orders of 
magnitude as obscuring material is expelled by outflows 
to gradually expose previously obscured AGNs. 

12. Since major mergers provide a natural process to 
form SMBH/AGN pairs, the scarcity of confirmed and 
candidate binary QSOs in large spectroscopic and VLBI 
surveys is unexpected. While two possibilities are raised 
in the literature, that either the SMBHs inspiral very 
rapidly, or fueling of the accretion disks is quenched dur- 
ing the binary phase, we propose a third hypothesis. 
Since the gas-depletion lifetime of ULIRGs and LIRGs 
is 10-100 times shorter than the time-scale for creation 
of a binary SMBH (>1 Gyr), it is possible that, in most 
instances, the gas will be consumed by star formation 
and accretion during a dual AGN phase long before the 
SMBHs inspiral to sub-pc separation. 
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Table 11 

Elliptical Apertures for Regions of Interest 



Table 12 

Photometric Data for Mrk 266 



Total System 



Region 


RA, DEC (J2000) 


a" 


b" 


PA° 


(1) 


(2) 


(3) 


(4) 


(5) 


SW 


(13:38:17.28,4-48:16:32.1) 


5.0 


3.0 


120 


NE 


(13:38:17.79,4-48:16:41.2) 


5.0 


3.0 


120 


Between Galaxies 


(13:38:17.51,448:16:36.8) 


6.0 


2.0 


115 


Northern Loop 


(13:38:17.54,448:16:56.0) 


17.0 


7.5 


10 


Total System 


(13:38:17.66,448:16:35.8) 


40.0 


40.0 





Note. — Parameters for elliptical ape rtu res used for photo- 
metric measurements, as illustrated in Fig. |31| 

Facilities: CXO (ACIS-S), XMM-Newton (EPIC pn, 
OM), GALEX (FUV, NUV), HST(ACS, NICMOS), 
UH:2.2m (), Spitzer(IRAC, MIPS, IRS) 

Objects: Mrk 266, Mrk 266 SW, Mrk 266 NE 

5. APPENDIX 

5.1. Photometric Data for Mrk 266 

Figure [31] displays regions parameterized in Table [TT] 
that were used to pe rfor m aperture photometry tabu- 
lated in Tables p] and fl3l as plotted in Figure |9|(jj2 e^- 




jicnskpc) 




Figure 31. The regions used to perform aperture photometry of 
the nuclei, central region, and Northern Loop superposed on the 24 
fim Spitzer MIPS image (upper left), the HST ACS B band image 
(upper right), the Chandra full-band X-ray image (lower left), and 
the 20 cm radio continuum image (lower right). All four images 
are shown at the same scale and orientation. 



5.2. Clarification of Regions Covered in Published 
Spitzer IRS Data for Mrk 266 

As noted in j ]2.7[ mid-IR spectroscopy of Mrk 266 ob- 
tained in staring mode with the Spitzer IRS is compli- 
cated due to aperture affects. To facilitate interpretation 
of the IRS spectroscopy published to date, the available 
observations are summarized here, and Table [14| provides 
clarification of which emission features have been mea- 
sured for each component of the system. 

In the low spectral resolution (R 65-130) observations 
(AOR #3755264; GTO program P00014; PI J. Houck), 
the Short-Low (SL, 5.2-14.5 //m) aperture of the IRS 



Bandpass 


v (Hz) 


U (Jy) 


*/„ (Jy) 


Ref. 


(i) 








{■>) 




(A\ 
\ L ) 


(5) 


2-10 keV (ASCA) 


1 


45E+18 


2 


48E-08 


5 


OE-09 


2 


2.0-7.0 keV (XMM) 


1 


09E+18 


2 


11E-08 


6 


3E-09 


1 


0.7-7 keV (ASCA) 


9 


31E+17 


3 


65E-08 


7 


3E-09 


2 


0.5-2.0 keV (XMM) 


3 


02E+17 


5 


96E-08 


1 


2E-08 


1 


GALEX FUV 


1 


96E+15 


7 


94E-04 


3 


1E-05 


1 


XMM-OM UVW2 


1 


41E+15 


2 


06E-03 


8 


0E-04 


1 


GALEX NUV 


1 


32E+15 


1 


63E-03 


3 


4E-05 


1 


XMM-OM UVW1 


1 


03E+15 


2 


12E-03 


2 


6E-04 


1 


XMM-OM U 


8 


71E+14 


3 


87E-03 


1 


4E-03 


1 


HST ACS F435W (B) 


6 


89E+14 


8 


92E-03 


4 


1E-05 


1 


V (Johnson) 


5 


42E+14 


1 


82E-02 


8 


6E-04 


3 


LW2 (ISOCAM) 


4 


44E+13 


1 


14E-01 


1 


4E-02 


4 


R (Johnson) 


4 


28E+14 


2 


86E-02 


8 


0E-04 


3 


HST ACS F814W (I) 


3 


68E+14 


2 


97E-02 


5 


5E-05 


1 


2MASS J (1.2 ,um) 


2 


43E+14 


5 


44E-02 


6 


5E-04 


1 


2MASS H (1.6 (im) 


1 


80E+14 


6 


42E-02 


8 


8E-04 


1 


2MASS Ks (2.2 //m) 


1 


39E+14 


6 


02E-02 


8 


3E-04 


1 


Spitzer 3.6 (im 


8 


44E+13 


3 


82E-02 


5 


7E-03 


1 


Spitzer 4.5 ^tm 


6 


67E+13 


3 


05E-02 


4 


6E-03 


1 


Spitzer 5.8 (im 


5 


23E+13 


6 


80E-02 


1 


OE-02 


1 


Spitzer 8.0 /im 


3 


81E+13 


1 


79E-01 


2 


7E-02 


1 


IRAS 12 fira 


2 


50E+13 


3 


20E-01 


1 


8E-02 


5 


Spitzer 24 /im 


1 


27E+13 


9 


55E-01 


1 


9E-01 


1 


IRAS 25 fim 


1 


20E+13 


1 


07E+00 


2 


8E-02 


5 


IRAS 60 urn 


5 


00E+12 


7 


25E+00 


3 


3E-02 


5 


Spitzer 70 /im 


4 


20E+12 


8 


95E+00 


1 


8E+00 


1 


IRAS 100 fira 


3 


00E+12 


1 


01E+01 


1 


4E-01 


5 


Spitzer 160 /^m 


1 


92E+12 


6 


41E+00 


1 


3E+00 


1 


170 Atm (ISO) 


1 


76E+12 


8 


73E+00 


2 


6E+00 


6 


SCUBA 850 ixm 


3 


53E+11 


8 


20E-02 


1 


7E-02 


7 


20.0 GHz (OVRO) 


2 


00E+10 


1 


49E-02 


1 


3E-03 


8 


4.85 GHz 


4 


85E+09 


4 


70E-02 


7 


1E-03 


9 


1.46 GHz (VLA) 


1 


46E+09 


1 


07E-01 


6 


2E-03 


8 


Texas 365 MHz 


3 


65E+08 


3 


53E-01 


4 


3E-02 


10 


151 MHz (6C) 


1 


52E+08 


5 


90E-01 


4 


OE-02 


11 


Note. — Reference codes 


in Column (5): l=New measurements 


in this article; 2=|Ueda et al. 1 12005); 


3= 


= Huchra (1977); 4= 


Ramos 


Almeida et ah] (2007); 


5= 


Sanders et al. (2003): 


6 


=|5tickel et al. 


(2004); 7=L)uimeeral. 


(200U); 8=Kdelson (1987) 


9 


=|necke~ et al. 


(19yi); 10= Douglas et al. 


T996t; TT 


=|riales et al. 


(1990). 





contained Mrk 266 NE (off center) , and missed both Mrk 
266 SW and the central radio/X-ray source; the Long- 
Low (LL, 14-38 /mi) slit contained Mrk 266 SW and the 
radio /X-ray source between the nuclei, but it missed Mrk 
266 NE. Neither aperture had the orientation required 
to sample the Norther n Loop. The low-re solution IRS 
spectrum displayed in Bra ndl et al.| ( |2006| ) is therefore 
a mixture of emission from Mrk 266 NE in the 5-14 /tm 
region with emission from Mrk 266 SW (and the radio/X- 
ray source between the nuclei) in the 14-38 /im region. 

AOR #3755264 also contained hi gh-resolution obser 
vations (spectral resolution ^600). iDudik et al. (2007 



Fig. 8) plotted the SH and LH apertures of these obser 
vations as an overlay on a 20 cm image, although they did 
not present the resulting spectra; their diagram clearly 
illustrates the aperture affects that complicate the in- 
terpretation of IRS spectra for this system. The SH slit 
contained only the radio source between the nuclei (miss- 
ing both nuclei) , while the LH slit contained 266 SW and 
the central radio source but it missed Mrk 266 NE. Ob- 
servations with the SH (9.9 - 19.6 /im) module were also 
obtained on 2005 January 10 in Spitzer program P03237 
(PI. E. Sturm, AOR #10510592); these data cover the 
SW nucleus, but not the NE nucleus. The available high- 
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Table 13 

Photometric Data for Mrk 266 



Components 



Bandpass v (Hz) /„ (Jy) or (Jy) Ref. 

(1) (2) (3) (4) ' (5) 



Table 14 

Regions Covered in Published IRS Measurements 



Feature 
(1) 



Flux (W m- 2 ) 
(2) 



EQW (fan) Ref. 
(3) (4) 



Between the Nuclei 



Mrk 266 SW 


2 0-7 f) keV 






2 


W ill \J o 


5 


1E- 


-10 




f) 5-2 keV 


3 


02E+17 


7 




\ 


5E- 


-09 




YUM OM TIVW2 




41E+15 




68E-04 


I 


1E- 


-05 




XMM-OM UVW1 


1 


03E+15 


2 


72E-04 


5 


4E- 


-05 


1 


"\7" TV T TV T /~\ TV J T T 

XMM-OM U 


8 


71E+14 


5 


12E-04 


2 


2E- 


-05 




HST ACS F435W (B) 


6 


89E+14 


1 


31E-03 


1 


2E- 


-05 




HST ACS F814W (I) 


3 


68E+14 


5 


03E-03 


2 


3E- 


-05 




2 MASS J (1.2 (um) 


2 


43E+14 


1 


04E-02 


2 


5E- 


-04 




NICMOS H (1.6 (im) 


1 


87E+14 


1 


45E-02 


1 


7E- 


-04 




2 MASS H (1.6 /tm) 


1 


80E+14 


1 


46E-02 


3 


3E- 


-04 




2 MASS Ks (2.2 /tni) 


1 


39E+14 


1 


53E-02 


3 


1E- 


-04 




Spitzcr 3.6 /xm 


8 


44E+13 


1 


50E-02 


2 


3E- 


-03 




Spitzer 4.5 /im 


6 


67E+13 


1 


47E-02 


2 


2E- 


-03 




Spitzer 5.8 /an 


5 


23E+13 


4 


59E-02 


6 


9E- 


-03 




Spitzcr 8.0 /im 


3 


81E+13 


1 


28E-01 


1 


9E- 


-02 




Spitzcr 24 /im 


1 


27E+13 


5 


91E-01 


1 


2E- 


-01 




Spitzer 70 (an 8 


4 


20E+12 


6 


01E+00 


1 


2E+00 




Spitzer 160 /tm a 


1 


92E+12 


4 


30E+00 


8 


6E- 


-01 




2 cm 


1 


50E+10 


4 


90E-03 


7 


0E- 


-04 


12 


6 cm 


4 


89E+09 


1 


01E-02 


5 


0E- 


-04 


12 


20 cm 


1 


47E+09 


3 


24E-02 


1 


4E- 


-03 


12 


Mrk 266 NE 


2.0-7.0 keV 


1 


09E+18 


1 


29E-08 


3 


9E- 


-09 




0.5-2.0 keV 


3 


02E+17 


9 


93E-09 


2 


0E- 


-09 




XMM-OM UVW2 


1 


41E+15 


2 


45E-04 


2 


0E- 


-05 




XMM-OM UVW1 


1 


03E+15 


3 


63E-04 


9 


1E- 


-05 




XMM-OM U 


8 


71E+14 


6 


49E-04 


3 


0E- 


-05 




HST ACS F435W (B) 


6 


89E+14 


1 


59E-03 


1 


5E- 


-05 




HST ACS F814W (I) 


3 


68E+14 


6 


02E-03 


2 


8E- 


-05 




2 MASS J (1.2 ^m) 


2 


43E+14 


9 


62E-03 


2 


4E- 


-04 




NICMOS H (1.6 /im) 


1 


87E+14 


1 


22E-02 


1 


6E- 


-04 




2 MASS H (1.6 (im) 


1 


80E+14 


1 


22E-02 


3 


0E- 


-04 




2 MASS Ks (2.2 /im) 


1 


39E+14 


1 


09E-02 


2 


6E- 


-04 




Spitzer 3.6 /tm 


8 


44E+13 


7 


35E-03 


1 


1E- 


-03 




Spitzer 4.5 /tm 


6 


67E+13 


5 


90E-03 


8 


9E- 


-04 




Spitzer 5.8 /tm 


5 


23E+13 


9 


57E-03 


1 


4E- 


-03 




Spitzcr 8.0 /tm 


3 


81E+13 


2 


23E-02 


3 


3E- 


-03 




Spitzcr 24 /tm 


1 


27E+13 


3 


28E-01 


6 


6E- 


-02 




Spitzer 70 /im a 


4 


20E+12 


1 


74E+00 


7 


0E- 


-01 




Spitzer 160 (im a 


1 


92E+12 


1 


25E+00 


5 


0E- 


-01 




2 cm 


1 


50E+10 


4 


10E-03 


6 


0E- 


-04 


12 


6 cm 


4 


89E+09 


1 


22E-02 


4 


0E- 


-04 


12 


20 cm 


1 


47E+09 


2 


35E-02 


1 


0E- 


-03 


12 



Mrk 266 SW 



2.0-7.0 keV 


1 


09E+18 


9 


19E- 


-10 


2 


8E- 


-10 


1 


0.5-2.0 keV 


3 


02E+17 


8 


60E- 


-09 


1 


7E- 


-09 


1 


2 cm 


1 


50E+10 


2 


70E- 


-03 


5 


0E- 


-04 


12 


6 cm 


4 


89E+09 


6 


00E- 


-03 


2 


7E- 


-03 


12 


20 cm 


1 


47E+09 


2 


15E- 


-02 


9 


0E- 


-04 


12 


Northern Loop 


0.5-2.0 keV 


3 


02E+17 


4 


63E- 


-09 


9 


3E- 


-10 




XMM-OM UVW2 


1 


41E+15 


1 


94E- 


-04 


1 


7E- 


-06 




XMM-OM UVW1 


1 


03E+15 


2 


10E- 


-04 


4 


8E- 


-06 




XMM-OM U 


8 


71E+14 


4 


28E- 


-04 


3 


0E- 


-05 




HST ACS F435W (B) 


6 


89E+14 


8 


84E- 


-04 


1 


6E- 


-05 




HST ACS F814W (I) 


3 


68E+14 


2 


50E- 


-03 


2 


1E- 


-05 




Spitzcr 24 /tm 


1 


27E+13 


3 


62E- 


-02 


7 


2E- 


-03 




20 cm 


1 


47E+09 


2 


60E- 


-03 


5 


0E- 


-05 


12 



Note. — The reference codes in Column (5) are as follows : l=New 
measurements in this article; 12=Mazzare lla et al.| ( [T988^ . 

a The 70 and 160 /tm flux densities for the SW and NE galaxies 
were estimated using the 8 and 24 /tm measurements as described 
in the text. 



[S IV] 10.51/im (SH) 


9.00(±0.38) 


X 


io- 


21 


B09 


PAH 10.6/tm (SH) 


0.25(±0.05) 


X 


lo- 


20 


B09 


PAH 10.7/mi (SH) 


0.12(±0.03) 


X 


ur 


20 


B09 


PAH 11.0/tm (SH) 


1.09(±0.05) 


X 


io- 


20 


B09 


PAH 11.3/tm (SH) 


14.90(±0.74) 


X 


10" 


20 


B09 


PAH 12.0/tm (SH) 


0.73(±0.10) 


x 


10" 


20 


B09 


H 2 (S2) 12.28/tm (SH) 


3.72(±0.32) 


X 


io- 


21 


B09 


HI (7 - 6) 12.37/im (SH) 


0.69(±0.05) 


X 


io- 


-21 


B09 


PAH 12.74/im (SH) 


7.13(±0.36) 


X 


io- 


20 


B09 


[Nell] 12.81/im (SH) 


57.04(±1.64) 


X 


io- 


21 


B09 


PAH 13.5/tm (SH) 


0.13(±0.06) 


X 


io- 


21 


B09 


PAH 14.2/tm (SH) 


0.57(±0.06) 


X 


io- 


20 


B09 


[Ne V] 14.32/tm (SH) 


7.96(±0.17) 


X 


io- 


21 


B09 


[CI II] 14.38/tm (SH) 


0.46(±0.07) 


X 


io- 


21 


B09 


[Ne III] 15.55/Ltm (SH) 


27.95(±0.76) 


X 


io- 


21 


B09 


PAH 16.45/tm (SH) 


0.84(±0.08) 


X 


10" 


20 


B09 


PAH 17.0/tm (SH) 


1.07(±0.14) 


X 


io- 


20 


B09 


H 2 (51) 17.03/tm (SH) 


8.40(±0.97) 


X 


io- 


-21 


B09 


[P III] 17.89/tm (SH) 


0.67(±0.15) 


X 


io- 


21 


B09 


[Fe II] 17.95/tm (SH) 


0.48(±0.14) 


X 


io- 


21 


B09 


[S III] 18.71/tm (SH) 


24.33(±0.69) 


X 


io- 


21 


B09 



Mrk 266 SW + Source Between Nucleus 



[Ne III] 15.55/im (LL) a 








-19 


PAH 17.0/tm (LL) 


0.51 


X 


io- 


[S III] 18.71/tm (LL) a 








-21 


[Fe III] 22.93/tm (LH) 


2.15(±0.29) 


X 


io- 


[Ne V] 24.31/tm (LH) b 


11.12(±0.59) 


X 


10" 


-21 


[O IV] 25.89/tm (LH) 


52.94(±5.66) 


X 


10" 


-21 


[O IV] 25.89/tm (LL) a 








-21 


[Fe II] 25.98/tm (LH) 


4.35(±0.70) 


X 


io- 


S 2 (50) 28.22/tm (LH) 


< 1.60 


X 


io- 


-21 


[5 III] 33.48/tm (LH) 


50.92(±2.91) 


X 


io- 


-21 


[5 ///] 33.48/tm (LL) a 








-21 


[Si II] 34.81/im (LH) 


87.02(±4.44) 


X 


io- 



403 



B06 
B06 
B06 
B09 
B09 
B09 
B06 
B09 
B09 
B09 
B06 
B09 



Mrk 266 NE 



PAH 6.2/im (SL) 


0.59 x 10- 


19 





619 


B06 


PAH 7.7/im (SL) 


0.92 x 10- 


1!) 





467 


B06 


PAH 8.6/im (SL) 









126 


B06 


H 2 5(3) 9.665/tm (SL) a 
[5 IV] 10.51//m (SH) a 










B06 










B06 


H 2 5(2) 12.279/tm (SL) a 










B06 


[Nell] 12.81/im (SL) a 




19 






B06 


PAH 11.3/tm (SL) 


0.47 x 10- 





422 


B06 


PAH 14.2/im (SL) 


0.08 x 10- 


1!) 





065 


B06 



Note. — Column (1): Spectral line identification followed in 
parentheses by the IRS module used: SL = 5.2-14.5 /tm, 3'.'7x57", 
R ~ 60-127; LL = 14.0-38.0 /im, ltf'5xl68'.'3, R ~ 57-126; SH = 
9.9-19.6 /im, 4'.'7xll'.'3, R ~ 600; LH = 18.7-37.3 /im, ll?lx22'.'3, 
R ~ 600. Column (2): Flux in W cm" 2 . Co lumn (3): Equivalen t 
width in ura. Column (4): Refere nce: B06 = |Brandl et al.| (2006); 
B09 = |Bernard-Salas et al.|{2009) . 

a These features are vis ible in the low-res spectrum published 
by Branch et al. (20061, but fluxes were not published. 
b TEis" 



is consistent with a 
1.19(±0.06) x 10~ 20 W cm- 2 



[Ne V] 24.32 /im flux of 
measured from the same LH 



aperture by Dudik et al. (2007 Table 2) 



resolution spectra were published in an atlas of Spitzer 
spectra for a sample of starburst galaxies ( Bernard-Salas 
etalp059 |. 
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Table 15 — Continued 



5.3. Evidence That Fe Ka Originates in Mrk 266 SW 

Another clue to the origin of the Fe Ka emission in 
Mrk 266 comes from analysis of the X-ray data in com- 
bination with optical spectroscopy. The ratio of X-ray 
(2-10 keV) continuum flux to (dereddened) [O III]A5007 
emission-line flux has been shown to be anti-correlated 
with Nh and with the equivalent width of the Fe Ka line, 
thus providing an in dicator for Compton-thick AGNs 
(Bassani et al. 1999). Usingthe F x (hard) values mea- 
sured from Chandra (Tabl e [2]) and the de-redde ned [O 
III] A5007 line fluxes from |Veilleux et aL] ([19951) ( 7 - 8 x 



1(T 14 and 8.5 x 1(T 15 erg IF 1 " for the SW and NE nu- 
clei), the F x /F[ in] ratios are 0.28 and 16 for the SW 
and NE nuclei, respectively. A ratio of F x /F[q hi] < 0.1 
indicates a highly absorbed source (Nh > 10 24 cm -2 ), 
whereas F X /F[ Q m ] >> 1 is a signature of a relatively 
unobscured (Type 1) AGN ( |Bassani et aLl|1999|) . This 
result adds compelling evidence that the Fe Ka line de- 
tected in the XMM spectrum of Mrk 266 originates from 
the SW (Seyfert 2) nucleus. Although Mrk 266 SW is 
a strong candidate for containing a Compton-thick X- 
ray source, observations at higher energies and greater 
sensitivity are required for confirmation. 



5.4. Star Clusters Photometry ({3.6) 



Table 15 

Star Cluster Photometry 



SC# S R.A. 


8 Dec 


B 


a 


M B 


I 


(7 


B - I 


(J 


ss.ss 


mm ss.s 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



001 


17 


19 


16 


32 


1 


21 


45 





03 


-14 


16 


21 


11 





05 


0.31 





06 


002 


17 


35 


16 


31 


9 


21 


52 





04 


-14 


09 


19 


48 





04 


2.01 





06 


003 


17 


14 


16 


33 


9 


22 


50 





01 


-13 


11 


21 


25 





08 


1.22 





09 


004 


17 


35 


16 


30 


9 


22 


91 





13 


-12 


70 


21 


20 





08 


1.68 





15 


005 


18 


10 


16 


13 


8 


23 


07 





01 


-12 


54 


22 


17 





03 


0.57 





03 


006 


17 


39 


16 


30 


9 


23 


08 





07 


-12 


53 


21 


21 





05 


1.84 





09 


007 


17 


16 


16 


36 


8 


23 


14 





05 


-12 


47 


22 


34 





09 


0.77 





10 


008 


17 


17 


16 


36 


5 


23 


24 





06 


-12 


37 


22 


58 





08 


0.63 





10 


009 


17 


72 


16 


34 


1 


23 


53 





02 


-12 


08 


22 


83 





03 


0.67 





01 


010 


17 


02 


16 


36 


6 


23 


54 





09 


-12 


07 


22 


49 





13 


1.02 





16 


011 


17 


08 


16 


31 


7 


23 


63 





01 


-11 


98 


20 


91 





02 


2.69 





04 


012 


16 


96 


16 


32 





23 


84 





09 


-11 


77 


23 


17 





20 


0.64 





22 


013 


16 


15 


16 


39 


9 


23 


97 





01 


-11 


64 


22 


65 





01 


1.29 





01 


014 


16 


79 


16 


14 


4 


24 


00 





03 


-11 


61 


23 


24 





04 


0.73 





05 


015 


16 


83 


16 


39 


2 


24 


09 





10 


-11 


52 


23 


87 





16 


0.19 





19 


016 


17 


68 


16 


38 


4 


24 


20 





04 


-11 


41 


23 


37 





16 


0.80 





16 


017 


17 


45 


16 


28 


8 


24 


23 





13 


-11 


38 


22 


76 





14 


1.44 





19 


018 


18 


64 


16 


44 


5 


24 


28 





06 


-11 


33 


23 


80 





07 


0.45 





09 


019 


17 


06 


16 


31 


9 


24 


30 





15 


-11 


31 


23 


23 





22 


1.04 





27 


020 


16 


80 


16 


36 


2 


24 


40 





07 


-11 


21 


23 


45 





13 


0.92 





15 


021 


16 


91 


16 


45 


6 


24 


54 





03 


-11 


07 


23 


82 





04 


0.69 





05 


022 


16 


71 


16 


10 


7 


24 


60 





03 


-11 


01 


23 


74 





04 


0.83 





05 


023 


18 


16 


16 


12 


3 


24 


71 





07 


-10 


90 


25 


32 





81 


-0.64 





81 


024 


17 


22 


16 


35 


2 


24 


72 





08 


-10 


89 


22 


91 





12 


1.78 





14 


025 


16 


62 


16 


45 


5 


24 


73 





01 


-10 


88 


23 


74 





04 


0.96 





06 


026 


17 


27 


16 


27 


1 


24 


74 





11 


-10 


87 


23 


43 





16 


1.28 





19 


027 


17 


29 


16 


28 


4 


24 


76 





26 


-10 


85 


23 


47 





32 


1.26 





41 


028 


17 


57 


16 


42 


1 


24 


83 





19 


-10 


78 


23 


83 





19 


0.97 





27 


029 


17 


19 


16 


30 


1 


24 


85 





12 


-10 


76 


25 


81 


1 


05 


-0.99 


1 


06 


030 


16 


32 


16 


11 





24 


86 





02 


-10 


75 


22 


23 





01 


2.60 





02 


031 


16 


77 


16 


15 


2 


24 


91 





01 


-10 


70 


24 


22 





06 


0.66 





07 


032 


16 


97 


16 


44 


1 


24 


92 





07 


-10 


69 


23 


74 





05 


1.15 





09 


033 


17 


09 


16 


37 


9 


24 


92 





15 


-10 


69 


23 


32 





16 


1.57 





22 


034 


16 


53 


16 


31 


3 


21 


94 





03 


-10 


67 


23 


81 





01 


1.10 





05 


035 


17 


23 


16 


40 


6 


24 


95 





08 


-10 


66 


24 


73 





22 


0.19 





23 



SC# <5 R.A. 


<5 Dec 


B 




a 


M B 


I 




(7 


B - I 


ex 




ss.ss 


mm ss.s 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


036 


18 


48 


16 


45 


2 


24 


97 





05 


-10.64 


24 


43 





07 


0.51 


0.09 


037 


18 


45 


16 


36 


6 


25 


00 





08 


-10.61 


24 


33 





21 


0.64 


0.22 


038 


18 


69 


16 


11 


3 


25 


00 





07 


-10.61 


24 


62 





13 


0.35 


0.15 


039 


17 


94 


16 


31 


3 


25 


05 





06 


-10.56 


24 


07 





07 


0.95 


0.09 


040 


18 


38 


16 


45 


9 


25 


12 





04 


-10.49 


24 


27 





08 


0.82 


0.09 


041 


16 


76 


16 


45 


6 


25 


13 





07 


-10.48 


24 


64 





11 


0.46 


0.13 


042 


18 


17 


17 


04 


5 


25 


14 





03 


-10.47 


24 


17 





03 


0.94 


0.04 


043 


18 


73 


16 


44 


9 


25 


16 





11 


-10.45 


24 


10 





11 


1.03 


0.16 


044 


16 


97 


16 


4.3 





25 


19 





07 


-10.42 


24 


11 





08 


1.05 


0.11 


045 


16 


70 


16 


37 


2 


25 


20 





13 


i n a 1 
-1U.41 


24 


57 





25 


n fin 


U.zo 


046 


16 


92 


16 


42 


2 


25 


24 





06 


-10.37 


24 


46 





07 


0.75 


0.09 


047 


17 


21 


16 


24 


6 


25 


28 





04 


-10.33 


24 


63 





06 


0.62 


0.07 


048 


16 


77 


16 


38 


7 


25 


32 





14 


-10.29 


24 


72 





31 


0.57 


0.34 


049 


17 


16 


16 


28 


1 


25 


34 





28 


-10.27 


24 


01 





34 


1.30 


0.44 


050 


17 


22 


16 


26 





25 


35 





07 


-10.26 


24 


30 





08 


1.02 


0.11 


051 


16 


42 


16 


30 


9 


25 


36 





05 


-10.25 


23 


23 





03 


2.10 


0.06 


052 


18 


32 


16 


43 


7 


25 


36 





09 


-10.25 


24 


48 





22 


0.85 


0.24 


053 


17 


20 


16 


13 


2 


25 


36 





06 


-10.25 


24 


24 





10 


1.09 


0.12 


054 


17 


00 


16 


34 


6 


25 


38 





22 


-10.23 


24 


57 





58 


0.78 


0.62 


055 


18 


57 


16 


51 


2 


25 


39 





04 


-10.22 


24 


32 





03 


1.04 


0.05 


056 


18 


56 


16 


44 


9 


25 


42 





07 


-10.19 


24 


91 





10 


0.48 


0.12 


057 


16 


77 


16 


32 


1 


25 


48 





27 


-10.13 


24 


40 





33 


1.05 


0.43 


058 


18 


66 


16 


45 


6 


25 


48 





13 


-10.13 


24 


41 





11 


1.04 


0.17 


059 


17 


13 


16 


14 


7 


25 


50 





06 


-10.11 


24 


80 





10 


0.67 


0.12 


060 


15 


85 


16 


42 


5 


25 


50 





04 


-10.11 


24 


53 





04 


0.94 


0.06 


061 


18 


77 


16 


12 


3 


25 


51 





08 


-10.10 


24 


97 





12 


0.51 


0.14 


062 


16 


20 


16 


33 


6 


25 


53 





04 


-10.08 


24 


37 





03 


1.13 


0.05 


063 


18 


41 


16 


37 


7 


25 


53 





09 


-10.08 


24 


57 





20 


0.93 


0.22 


064 


16 


72 


16 


35 


3 


25 


53 





18 


-10.08 


23 


99 





19 


1.51 


0.26 


065 


18 


36 


16 


45 


1 


25 


54 





07 


-10.07 


25 


06 





16 


0.45 


0.17 


066 


17 


00 


16 


43 





25 


59 





10 


-10.02 


24 


76 





16 


0.80 


0.19 


067 


17 


63 


16 


31 


7 


25 


60 





16 


-10.01 


24 


06 





19 


1.51 


0.25 


068 


19 


08 


16 


13 


2 


25 


62 





04 


-9.99 


24 


03 





02 


1.56 


0.04 


069 


16 


63 


16 


46 


2 


25 


64 





11 


-9.97 


24 


75 





13 


0.86 


0.17 


070 


17 


33 


16 


10 


6 


25 


65 





18 


-9.96 


24 


68 





17 


0.94 


0.25 


071 


18 


73 


16 


43 


8 


25 


65 





22 


-9.96 


25 


84 





51 


-0.22 


0.56 


072 


16 


64 


16 


38 


9 


25 


66 





09 


-9.95 


24 


82 





14 


0.81 


0.17 


073 


16 


09 


16 


31 


2 


25 


68 





04 


-9.93 


24 


51 





04 


1.14 


0.06 


074 


17 


94 


16 


37 





25 


70 





10 


-9.91 


24 


86 





30 


0.81 


0.32 


075 


18 


06 


16 


35 


1 


25 


70 





12 


-9.91 


25 


18 





29 


0.49 


0.31 


076 


15 


58 


16 


41 





25 


72 





05 


-9.89 


24 


68 





04 


1.01 


0.06 


077 


20 


51 


16 


18 


3 


25 


77 





05 


-9.84 


24 


57 





04 


1.17 


0.06 


078 


17 


70 


16 


31 


7 


25 


77 





15 


-9.84 


25 


10 





33 


0.64 


0.36 


079 


18 


29 


16 


13 


1 


25 


80 





12 


-9.81 


24 


49 





22 


1.28 


0.25 


080 


18 


74 


16 


43 


6 


25 


82 





20 


-9.79 


26 


29 





62 


-0.50 


0.65 


081 


18 


50 


16 


18 


8 


25 


92 





08 


-9.69 


25 


25 





08 


0.64 


0.11 


082 


16 


72 


16 


15 


3 


25 


95 





22 


-y.Du 


25 


45 





27 


C\ A"7 


U.oo 


083 


16 


75 


16 


37 


9 


26 


02 





25 


-9.59 


24 


63 





19 


1.36 


0.31 


084 


16 


50 


16 


16 


1 


26 


06 





07 


-9.55 


24 


71 





06 


1.32 


0.09 


085 


17 


31 


16 


39 


9 


26 


08 





20 


-9.53 


25 


29 





19 


0.76 


0.28 


086 


16 


56 


16 


45 


3 


26 


08 





12 


-9.53 


25 


45 





20 


0.60 


0.23 


087 


16 


86 


16 


45 


9 


26 


12 





13 


-9.49 


25 


54 





23 


0.55 


0.26 


088 


18 


54 


16 


45 


3 


26 


13 





14 


-9.48 


25 


81 





28 


0.29 


0.31 


089 


16 


09 


16 


46 


2 


26 


17 





09 


-9.44 


24 


96 





07 


1.18 


0.11 


090 


18 


01 


16 


28 


7 


26 


19 





11 


-9.42 


24 


78 





09 


1.38 


0.14 


091 


16 


85 


16 


38 


1 


26 


19 





10 


-9.42 


25 


87 





63 


0.29 


0.75 


092 


16 


11 


16 


18 


2 


26 


20 





08 


-9.41 


25 


42 





11 


0.75 


0.14 


093 


19 


64 


16 


13 


9 


26 


23 





07 


-9.38 


25 


11 





05 


1.09 


0.09 


094 


17 


86 


16 


32 


2 


26 


24 





13 


-9.37 


24 


23 





08 


1.98 


0.15 


095 


16 


14 


16 


12 


6 


26 


27 





09 


-9.34 


25 


22 





09 


1.02 


0.13 


096 


17 


20 


16 


26 


1 


26 


29 





16 


-9.32 


26 


50 





79 


-0.24 


0.81 


097 


19 


10 


16 


36 





26 


30 





09 


-9.31 


25 


37 





09 


0.90 


0.13 


098 


16 


90 


16 


27 





26 


30 





18 


-9.31 


24 


45 





12 


1.82 


0.22 


099 


16 


73 


16 


35 


7 


26 


31 





11 


-9.30 


24 


08 





21 


2.20 


0.46 


100 


18 


68 


16 


50 


1 


26 


43 





09 


-9.18 


25 


79 





12 


0.61 


0.15 


101 


17 


50 


16 


16 


2 


26 


16 





09 


-9.15 


25 


47 





08 


0.96 


0.12 


102 


16 


45 


16 


28 


4 


26 


50 





11 


-9.11 


26 


66 





39 


-0.19 


0.41 


103 


15 


87 


16 


39 


3 


26 


51 





09 


-9.10 


25 


41 





09 


1.07 


0.13 


101 


11 


99 


16 


28 


8 


26 


53 





09 


-9.08 


25 


11 





07 


1.36 


0.11 


105 


18 


08 


16 


17 


9 


26 


55 





10 


-9.06 


25 


66 





10 


0.86 


0.14 


106 


20 


16 


16 


24 


9 


26 


56 





10 


-9.05 


25 


16 





06 


1.37 


0.12 


107 


18 


69 


16 


45 


9 


26 


56 





36 


-9.05 


25 


34 





27 


1.19 


0.45 


108 


16 


41 


16 


35 





26 


58 





13 


-9.03 


25 


16 





13 


1.39 


0.18 


109 


17 


49 


16 


11 


3 


26 


63 





09 


-8.98 


25 


69 





10 


0.91 


0.13 


110 


15 


71 


16 


21 


7 


26 


63 





10 


-8.98 


25 


49 





09 


1.11 


0.13 



Dual AGNs, Outflows, Shocks, and Young Star Clusters in Mrk 266 
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Table 15 — Continued 



SC# 5 R.A. 


S Dec 


j3 


a 


iliB 


j 


a 


B — I 


a 




ss.ss 


111111 BB.O 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


(1) 


(2) 




(A\ 


(K\ 
V>> 




V 1 ) 


\°) 


(Q~\ 


(, iu ; 


111 


15.77 


1 fi AA ^ 


ZD. i D 


n 1 9 
u. 1Z 


-o.oo 


OA R9 


u.uo 


9 11 




112 


17.22 


16 44.0 


26.77 


0.21 


-8.84 


25.77 


0.41 


0.97 


0.46 


113 


16.42 


16 46.4 


26.85 


0.16 


-8.76 


25.87 


0.19 


0.95 


0.25 


114 


16.04 


16 21.8 


26.89 


0.12 


-8.72 


26.14 


0.16 


0.72 


0.20 


115 


16.78 


16 57.1 


27.07 


0.18 


-8.54 


25.73 


0.12 


1.31 


0.22 


116 


18.37 


16 40.5 


27.12 


0.88 


-8.49 


24.25 


0.28 


2.84 


0.92 


117 


15.93 


16 50.7 


27.13 


0.17 


-8.48 


25.98 


0.13 


1.12 


0.21 


118 


15.67 


16 48.2 


27.16 


0.13 


-8.45 


25.59 


0.09 


1.54 


0.16 


119 


17.56 


16 56.9 


27.23 


0.22 


-8.38 


25.73 


0.13 


1.47 


0.26 



Note. — Column (1): Star cluster (candidate) number. Cols. 
(2) and (3): J2000 S R.A. and 8 Dec of the object. To save space, 
offsets in seconds of time in R.A. and offsets in arcminutes and arc- 
seconds in Dec are listed. Full coordinates are obtained by adding 
these values to R.A.=13h38m and Dec=48d (e.g., 13h38ml7.19s; 
48dl6m32.4s). The la uncertainty in the coordinates is ±0'/3. 
The objects can be identified as "Mrk 266 SC###" or as "GOALS 
Jhhmmssss+ddmmsss." For example: Mrk 266 SC001 = GOALS 
J13381719+4816324. Col. (4): apparent mag in the HST ACS 
F435W (B) filter. Col. (5): la uncertainty in the B mag. Col. 
(6): absolute B mag. Col. (7): apparent mag in the HST ACS 
F814W (I) filter. Col. (8): la uncertainty in the I mag. Col. (9): 
B - I color mag. Col. (10): la uncertainty in the B - I color mag. 
Absolute B mag and B-I colors have been corrected for fore ground 
Galactic extinction using values from |Schlegel et al.|jl"9 98) as pro- 
vided by NED: A B = 0.056 mag; At = U7U25 mag' The se dat a 
have not been corrected for extinction within Mrk 266 (see j]3.6.2|l. 
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